Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,346,6750TG42.9% 132.8 / 13.4 77T218P (ACC→CCC) fumBanaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base T (7/37);  new base G (33/0);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 5.75e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAT  >  NC_000913/4346527‑4346825
                                                                                                                                                     |                                                                                                                                                      
aaCGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTAc                                                                                                                                                      <  2:622632/151‑1 (MQ=255)
  cGCCTGACCCTGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGc                                                                                                                                                    <  1:182179/151‑1 (MQ=255)
   gCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGTGCAGGACCGGGGGGAc                                                                                                                                                      >  2:188538/1‑148 (MQ=255)
   gCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTAACCAGGGGACg                                                                                                                                                     >  2:502188/1‑149 (MQ=255)
       gACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCAtt                                                                                                                                                 <  1:856862/149‑1 (MQ=255)
          cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTtctc                                                                                                                                            <  1:1564/151‑1 (MQ=255)
          cATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTtctc                                                                                                                                            <  2:920309/151‑1 (MQ=255)
           aTGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTcg                                                                                                                                           <  1:172646/151‑1 (MQ=255)
            tGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGGGCCGCAACGGGGGGACGCATTTtct                                                                                                                                             >  1:183689/1‑148 (MQ=255)
                cgttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTTTcgacg                                                                                                                                        >  2:411328/1‑149 (MQ=255)
                  ttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGTGCAGTACCGAGGGGACGCATTTTCTCGACGAg                                                                                                                                      >  2:327197/1‑149 (MQ=255)
                  ttcCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAgg                                                                                                                                     >  1:372728/1‑150 (MQ=255)
                    cccTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTTGTTTCCGAAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGa                                                                                                                                    <  2:4514/149‑1 (MQ=255)
                      cTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCCGGACCGGGGGGACGCATTTTTTTGACGAGGAAGtt                                                                                                                                >  2:645110/1‑151 (MQ=255)
                          cGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTc                                                                                                                            <  1:970902/151‑1 (MQ=255)
                              ggCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGtt                                                                                                                        <  2:188296/151‑1 (MQ=255)
                               gCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGttt                                                                                                                       <  1:815073/151‑1 (MQ=255)
                                     cATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCggg                                                                                                                 <  2:105444/151‑1 (MQ=255)
                                      aTCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGGTCAGGC‑GCCATACCGGGGGGGCGCAATTTCTCTACGAGGGAGTTTTTTAGTTTgatgg                                                                                                                  >  2:920595/1‑145 (MQ=255)
                                          tAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑GAGGCGGCCGTACCGCGGGGACGCCTTTTTTTGACGAGGGAGGTTTTTCGGTTGTCggggggc                                                                                                             >  1:222731/1‑146 (MQ=255)
                                           aaTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTca                                                                                                            <  1:743978/150‑1 (MQ=255)
                                            aTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTcagc                                                                                                          <  1:155495/151‑1 (MQ=255)
                                            aTAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCCGCAGTACCGAGGGGACGCCTTTTTTCGGCGAGGAAGGTTTTCAGTTTGCCGGGGGGcagc                                                                                                          >  2:918544/1‑151 (MQ=255)
                                             tAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGGGGCGTACCGAGGGGAAGCCATTTCTCGCCGGGGGAGGTTTTTAGTTTGCCGGGAGTgag                                                                                                           >  2:351889/1‑149 (MQ=255)
                                                  aGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGGAGGTTTTTAGTTTGCCGGGAGGCAGCAggg                                                                                                      >  2:651386/1‑149 (MQ=255)
                                                    cgcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCttt                                                                                                  <  1:210094/151‑1 (MQ=255)
                                                     gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTAACGAGGGGACGCCTTTTTTCGGCGAGGGAGGTTTTTAGTTTTCCGGGGGTCAGCAGGGCttt                                                                                                  >  1:875908/1‑150 (MQ=255)
                                                     gcTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGCAACGAGGGGACGCATTTTTTCGCCGAGGGAGGTTTTCCGTTTGGCGGGGGTCCACAggggctt                                                                                                  >  1:161145/1‑146 (MQ=255)
                                                       ttGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CCAGCTGCAGTCCCGACGGGACGCTATTTCTTCCCCCGGAAGGTTTGCAGGTTGTCGGGGGGCCGCAGGGGttttg                                                                                                >  1:6971/1‑148 (MQ=255)
                                                        tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGGACCGAGGGGACGCATTTTCTTGACGAGGGAGGGTTTTAGTTTGTCCGGAGGCAGCAGGGCttttggt                                                                                              >  1:293046/1‑147 (MQ=255)
                                                        tGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGGACCCCGGGGACGCATTTTCTTTACGGGGAAGTTTTTTAATTTGCCGGGAGGCCGCCGGGGTTTGGtt                                                                                              >  2:765356/1‑151 (MQ=255)
                                                          cTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCCGCAGTACCGCGGGGACGCATTTTCTCTACGAGGGAGGTTTTCAGTTTGCCGGGGGTCCGCAGGGCCTTgggt                                                                                              >  1:37225/1‑147 (MQ=255)
                                                              cTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGTGCCGGACCGGGGGGACGCCTTTTTTCGCCGAGGAAGGTTTTCCGTTTGCCGGGGGTCAGCAGGGGTTTGGgtt                                                                                             >  1:271100/1‑146 (MQ=255)
                                                                tgccGGCTTCCAGGTTTGTTTCCGCAGACGTTCAACCAATCATAAACGCGATAATGTACGGCGTG‑CATTCTGCAGTACCGAGGGTACTCATTTTCTCGACGAGGAAGTTATACAGTTTGCCGGGAGTCTGCAGGGCTTTGGTTTCCTgg                                                                                        <  1:645913/146‑1 (MQ=255)
                                                                    ggTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑GAGGCTGCAGTACCGAGGGTAAGCATTTTCTCGACGGGGGAGGTTTTTAGGTTGTCGGGGGGCAGCAGGGCTTTGGTTTTCTggggg                                                                                     >  2:960679/1‑145 (MQ=255)
                                                                     gTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGAt                                                                                 <  2:68727/151‑1 (MQ=255)
                                                                         tCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACg                                                                              <  1:503443/150‑1 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCTCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGt                                                                             <  1:80645/150‑1 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGt                                                                             <  1:502106/150‑1 (MQ=255)
                                                                          cAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGt                                                                             <  2:899451/150‑1 (MQ=255)
                                                                               tggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTttgt                                                                        <  2:531025/150‑1 (MQ=255)
                                                                                ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑GAGGCTGCAGTCACGCGGGGACCCATTTTCTCGACGGGGAAGTTTTTTAAGTTTTCGGGGGGCAACAGGGGTTTGGTTTCCTTGGAGAGGTACGGTttggt                                                                       >  2:302669/1‑148 (MQ=255)
                                                                                ggttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTCTCGACGAGGAAGTTTTTTCGTTTGCCGGGGGTCAGCAGGGCTTTGGGTTCCTGGTAGAGATACGTTttgtt                                                                       >  2:225988/1‑150 (MQ=255)
                                                                                 gttTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCATTTTTTCCACGAGGAAGGTTTTTCGTTTGCCGGGGGTCCGCCGGGCCTTTGTTTCCCGGGAGAGATAAGtttttt                                                                        >  2:751683/1‑146 (MQ=255)
                                                                                   ttCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGGTGCAGGAACGCGGGGACGCATTTTTTCGACGGGGGGGGTTTTCAGTTTTGCGGGGGTCAGCAGGGGGTTGGTTTCCTGGTAAAGGTACGGTttgtttgg                                                                    >  1:714418/1‑147 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCag                                                                  <  2:937008/150‑1 (MQ=255)
                                                                                     ccGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCag                                                                  <  2:204773/150‑1 (MQ=255)
                                                                                        cAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGGACGCCTTTTTTCGACGAGGAAGTTTTTTAGTTTGCCGGGGGTCAGCAGGGCTTTGGGTTTCTGGTAGAGATGCGTTTTTTTGGCAGAgcc                                                              >  2:993957/1‑151 (MQ=255)
                                                                                            cGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTAACCAGGGGACGCATTTTCTCGACGAGGGAGGTTTTCAGTTTGCCGGGAGTCAGCCGGGGTTTGGGTTTCTGGTAGAGATACGTTTTGTTGGCAGAgcc                                                              >  2:77079/1‑147 (MQ=255)
                                                                                              tACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCCAGGGGACGCAATTTTTCGACGAGGGAGTTTTTTTGGTTTGCGGGGGGCAGCCGGGGTTTTGGTTTCTGGGAGAGATACGGTTTTTTGGGAGAgccgccgcc                                                        >  1:759458/1‑151 (MQ=255)
                                                                                               aCCGCTAATCACAAACGAGATATGGTACGGGGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGGGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCGGGTAGAGATACGTTTTGTTGGCAGAgccgccgc                                                         <  1:405888/149‑1 (MQ=255)
                                                                                               aCCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTAACCAGGGGACGCATTTTCTCGACGAGGGAGGTTTTTAGTTTTCCGGGGGTCAGCCGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGGAGAgccgccgcc                                                        >  2:365840/1‑150 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGG‑CAGGCCTTAGTACCGAGGGTACGCCTTTTTTCCACGCGGAAGTTTTTTTATTTTGCGGGAGACATTCGGGCGTTGGTTTCCTGGTAGAGGTATGTTTTTTTTGGAGAGCCGCCGCCttt                                                     >  2:922815/1‑149 (MQ=255)
                                                                                                   ccAATCACAAACGCGATATGGTACGGCGGG‑CAAGGTGGAGGAACGAGGGGACGCCATTTTTTGACGAGGGAGGTTTTTAGGTTGCCGGGGGGTAGCCGGGCTTTGGTTTTCTGGGAGAGATAAGGTTTGTTGGGGGAGCCGCCGCCttt                                                     >  2:459733/1‑149 (MQ=255)
                                                                                                      aTCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaa                                                <  1:174791/151‑1 (MQ=255)
                                                                                                        cacaAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacg                                              <  2:463130/151‑1 (MQ=255)
                                                                                                        cacaAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTcgcaacg                                              <  1:556634/151‑1 (MQ=255)
                                                                                                           aaaCGCGATATGGTACGGCGGG‑CAGGCCGCAGGACCCAGGGGACGCATTTTTTTGACGAGGAAGGTTTTTCGGTTGCCGGGAGGCAGACGGGCTTTGGTTTTCTGGTTGAGGTTCGGTTTTTTTGCAGAGCCGCCCCCTTTcgccacg                                              >  1:371225/1‑148 (MQ=255)
                                                                                                            aaCGCGATATGGTACGGCGGG‑CCGGCTGCAGTACCCGGGGGACGCATTTTTTCGACGAGGGAGTTTTTTAGTTTGCCGGGGGGCAGCAGGGCTTTGGTTTTCTGGGAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTtgcaacgcaa                                           >  2:264667/1‑150 (MQ=255)
                                                                                                               gcgATATGGTACGGCGGG‑CAGGCTGCCGTACCGCGGGGACGCCATTTTTTGCCGAGGGAGTTTTTTAGTTTGCCGGGGGGCAGCAGGGCCTTGGTTTCCTGGTAGAGATACGTTTTTTTGGCAGCGCCGCCCCCCTTcgccacgc                                             >  1:310282/1‑145 (MQ=255)
                                                                                                                 gATATGGTACGGCGGG‑CCGGGTGCAGGACCGAGGGGGAGCAATTTTTTGACGGGGAAGGTTTTTTGTTTGCCGGGGGTGAGCCGGGGTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCCTTCGCAACGCCAAgg                                        >  1:358265/1‑148 (MQ=255)
                                                                                                                 gATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAAt                                     <  2:616963/151‑1 (MQ=255)
                                                                                                                  atatGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGTCGCCGCCTTTCGCAACGCAAAGGAAtt                                    <  1:637016/151‑1 (MQ=255)
                                                                                                                      ggTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGt                                 <  1:896341/150‑1 (MQ=255)
                                                                                                                         aCGGCGGG‑GAGGCCGGAGTACCCGGGGGACCCATTTTTTCGCCGAGGAAGGTTTTTAGTTTTCCGGGAGTCAGCCGGGCTTTGGGTTCCTGGGAGAGGTACGTTTTGTTGGCAGAGCCGCCGCCCTTTGCAACGCAACGGAATTTGTACTc                             >  1:713242/1‑151 (MQ=255)
                                                                                                                                   aGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATc                    >  2:497166/1‑151 (MQ=255)
                                                                                                                                     gCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTa                  <  1:787027/151‑1 (MQ=255)
                                                                                                                                        gCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTAcct               <  1:247717/151‑2 (MQ=255)
                                                                                                                                         cAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTa                                >  2:264405/1‑133 (MQ=255)
                                                                                                                                         cAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTa                                <  1:264375/133‑1 (MQ=255)
                                                                                                                                          aGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcgcg             >  2:67327/1‑151 (MQ=255)
                                                                                                                                           gTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACcgcg             <  1:805481/150‑1 (MQ=255)
                                                                                                                                             aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTAc          <  2:734550/151‑1 (MQ=255)
                                                                                                                                             aCCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTAc          >  1:79840/1‑151 (MQ=255)
                                                                                                                                              ccGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACa         <  1:574887/151‑1 (MQ=255)
                                                                                                                                                  gggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTc     <  2:478995/151‑1 (MQ=255)
                                                                                                                                                   ggTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCg    >  1:242252/1‑151 (MQ=255)
                                                                                                                                                     tACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAt  >  1:354889/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AACGCCTGACCATGTTCGTTCCCTTCCGTCGGCAGTTCATCGTAATAGTGAGCGCTTGCTAACTTGACGGTTTTCAGGTTGGTTTCCGCAGACGTACCGCCAATCACAAACGCGATATGGTACGGCGGG‑CAGGCTGCAGTACCGAGGGTACGCATTTTCTCGACGAGGAAGTTTTTCAGTTTGCCGGGAGTCAGCAGGGCTTTGGTTTCCTGGTAGAGATACGTTTTGTTGGCAGAGCCGCCGCCTTTCGCAACGCAAAGGAATTTGTACTCATCGCCATCTACCGCGTACAGGTCGAT  >  NC_000913/4346527‑4346825

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑