Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,381,504 T→C E272G (GAA→GGA)  frdA ← anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,381,5040TC99.1% 398.0 / ‑5.0 107E272G (GAA→GGA) frdAanaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (54/51);  total (55/52)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 109 total aligned reads displayed.

GCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGG  >  NC_000913/4381356‑4381652
                                                                                                                                                    |                                                                                                                                                    
gCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTccc                                                                                                                                                    >  3:400265/1‑151 (MQ=255)
 ccGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCg                                                                                                                                                   >  3:456709/1‑151 (MQ=255)
   gAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCggg                                                                                                                                                 >  2:132705/1‑151 (MQ=255)
   gAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCAGGTTCCATATATTTGTTTTTCGGCCCCCCCAGCGGAGTTCCCggg                                                                                                                                                 >  2:800034/1‑151 (MQ=255)
        gACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCAt                                                                                                                                            <  4:354187/151‑1 (MQ=255)
        gACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCAt                                                                                                                                            <  3:129027/151‑1 (MQ=255)
         aCGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATg                                                                                                                                           >  4:40122/1‑151 (MQ=255)
           gCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGcc                                                                                                                                         >  4:42025/1‑151 (MQ=255)
            cAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCg                                                                                                                                        <  2:1169141/151‑1 (MQ=255)
             aaGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGt                                                                                                                                       <  2:1209072/151‑1 (MQ=255)
              aGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGt                                                                                                                                       >  2:1284223/1‑150 (MQ=255)
                 cGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATc                                                                                                                                   >  2:313509/1‑151 (MQ=255)
                    gaTAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTg                                                                                                                                <  1:1132394/151‑1 (MQ=255)
                        aaCCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGa                                                                                                                            <  2:137743/151‑1 (MQ=255)
                          ccACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGAt                                                                                                                           >  1:255086/1‑150 (MQ=255)
                           cacaTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATaa                                                                                                                         <  2:1201383/151‑1 (MQ=255)
                           cacaTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATa                                                                                                                          >  2:823414/1‑150 (MQ=255)
                             caTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAAc                                                                                                                        <  2:130817/150‑1 (MQ=255)
                                cGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTa                                                                                                                    <  3:118443/151‑1 (MQ=255)
                                 gCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGt                                                                                                                     >  1:231056/1‑149 (MQ=255)
                                     cgcgGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAt                                                                                                               <  1:295712/151‑1 (MQ=255)
                                      gcggcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAtt                                                                                                              <  3:440738/151‑1 (MQ=255)
                                        ggcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAtttt                                                                                                            >  4:377257/1‑151 (MQ=255)
                                        ggcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAtttt                                                                                                            >  1:464261/1‑151 (MQ=255)
                                        ggcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAtttt                                                                                                            >  4:348431/1‑151 (MQ=255)
                                         gcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAttttt                                                                                                           >  1:1091475/1‑151 (MQ=255)
                                         gcgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAttttt                                                                                                           <  1:702301/151‑1 (MQ=255)
                                          cgTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCAttttt                                                                                                           >  2:5340/1‑150 (MQ=255)
                                           gTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgt                                                                                                         <  2:1053213/151‑1 (MQ=255)
                                           gTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgt                                                                                                         <  2:1228729/151‑1 (MQ=255)
                                           gTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTttgt                                                                                                         >  2:999359/1‑151 (MQ=255)
                                                      gtTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCGGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAAt                                                                                               <  2:272326/150‑1 (MQ=255)
                                                           cTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGcc                                                                                         >  3:46664/1‑151 (MQ=255)
                                                             ttACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCtt                                                                                       >  2:46792/1‑151 (MQ=255)
                                                              tACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCTGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCtt                                                                                       <  4:233775/150‑1 (MQ=255)
                                                               aCGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGAAGCCATTTTTGTTGACCAGAATACCGCCtt                                                                                       >  4:182023/1‑149 (MQ=255)
                                                               aCGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTACAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTc                                                                                      >  2:468377/1‑150 (MQ=255)
                                                                 gCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCAcc                                                                                   >  1:887420/1‑151 (MQ=255)
                                                                  ccATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCAcc                                                                                   >  2:421509/1‑150 (MQ=255)
                                                                   cATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACcgc                                                                                 <  1:165815/151‑1 (MQ=255)
                                                                     ttCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGTAGTTCCCGGGCCCATGCTGTAATCTTGCAGATAACGGTAGCCCTTTTTGTTGACCAGAATACCGCCTTCACcgcg                                                                                <  3:171974/150‑1 (MQ=255)
                                                                      tCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCAATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGc                                                                              >  2:112224/1‑151 (MQ=255)
                                                                       cGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGc                                                                              <  1:731431/150‑1 (MQ=255)
                                                                        gTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCa                                                                             >  1:1072045/1‑150 (MQ=255)
                                                                            cAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCtt                                                                        <  2:181599/151‑1 (MQ=255)
                                                                              gAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCg                                                                      <  1:312838/151‑1 (MQ=255)
                                                                              gAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCg                                                                      <  1:951197/151‑1 (MQ=255)
                                                                               aaGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCAATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCgg                                                                     >  3:262855/1‑151 (MQ=255)
                                                                                 ggCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGtc                                                                   >  1:406922/1‑151 (MQ=255)
                                                                                  gCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGtc                                                                   >  1:1222454/1‑150 (MQ=255)
                                                                                    cTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGtcatc                                                                >  1:1081973/1‑151 (MQ=255)
                                                                                      gagagaCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGtcatca                                                               <  1:296710/150‑1 (MQ=255)
                                                                                      gagagaCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAg                                                              <  1:1066722/151‑1 (MQ=255)
                                                                                       agagaCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCGGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAgg                                                             >  4:255157/1‑151 (MQ=255)
                                                                                        gagaCTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGa                                                            <  1:1098923/151‑1 (MQ=255)
                                                                                                   gcgTGGACCCAGTTCCATATCTTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCCTGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCt                                                 <  4:125961/151‑1 (MQ=255)
                                                                                                   gcgTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCt                                                 >  3:256096/1‑151 (MQ=255)
                                                                                                     gTGGACCCAGTTACATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGACCATGCCGTAATATTGCAGATAACGGGAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTg                                                >  2:366878/1‑150 (MQ=255)
                                                                                                      tGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTgg                                               <  1:1072947/150‑1 (MQ=255)
                                                                                                      aggACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAACCATTTTTGTTGACCAGAATACCGCCTTCAACGCGGCAACCCTCGGTCATCAGGATACCGGAACCTg                                                >  2:839765/2‑149 (MQ=255)
                                                                                                       ggACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCa                                             <  2:694401/151‑1 (MQ=255)
                                                                                                       ggACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCa                                             <  1:284785/151‑1 (MQ=255)
                                                                                                       ggACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCa                                             <  2:1254606/151‑1 (MQ=255)
                                                                                                        gACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGc                                              >  2:511695/1‑149 (MQ=255)
                                                                                                        gACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAAACTGGc                                              >  2:701056/1‑149 (MQ=255)
                                                                                                         aCCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGa                                           <  1:447221/151‑1 (MQ=255)
                                                                                                          cccAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGAc                                          >  1:1020771/1‑151 (MQ=255)
                                                                                                           ccAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGAcc                                         <  4:377238/151‑1 (MQ=255)
                                                                                                             aGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCgg                                       <  2:210823/151‑1 (MQ=255)
                                                                                                              gTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCg                                        >  3:445005/1‑149 (MQ=255)
                                                                                                              gTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGt                                      >  2:714749/1‑151 (MQ=255)
                                                                                                               ttCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGt                                      <  1:31182/150‑1 (MQ=255)
                                                                                                               ttCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGt                                      <  2:664907/150‑1 (MQ=255)
                                                                                                                  cATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTggg                                  <  4:26772/151‑1 (MQ=255)
                                                                                                                        tttGTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATa                             <  4:173715/150‑1 (MQ=255)
                                                                                                                        tttGTTTTTCGCCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGAt                              <  1:1023630/149‑1 (MQ=255)
                                                                                                                           gTTTTTCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTg                          >  2:404649/1‑150 (MQ=255)
                                                                                                                            tttttCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGa                         >  1:365322/1‑150 (MQ=255)
                                                                                                                              tttCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACg                      >  1:730056/1‑151 (MQ=255)
                                                                                                                               ttCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACg                      >  3:103546/1‑150 (MQ=255)
                                                                                                                               ttCGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGGGATACTGAACg                      >  4:126455/1‑150 (MQ=255)
                                                                                                                                 tggCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAAt                   >  2:623735/2‑151 (MQ=255)
                                                                                                                                 cGGCTCGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAAt                   <  4:5915/151‑1 (MQ=255)
                                                                                                                                      cGCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCAt              <  1:18126/151‑1 (MQ=255)
                                                                                                                                       gCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATg             >  4:165853/1‑151 (MQ=255)
                                                                                                                                       gCCCAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATg             >  3:393018/1‑151 (MQ=255)
                                                                                                                                        cccAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATg             <  3:177378/150‑1 (MQ=255)
                                                                                                                                        cccAGCGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGt            <  1:467120/151‑1 (MQ=255)
                                                                                                                                           aGCGGAGTTACCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCCTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCAc         >  3:198584/1‑151 (MQ=255)
                                                                                                                                             cGGAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCAc         <  2:71316/149‑1 (MQ=255)
                                                                                                                                              ggAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGc       <  2:661333/150‑1 (MQ=255)
                                                                                                                                              ggAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGc       >  4:50119/1‑150 (MQ=255)
                                                                                                                                              ggAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCa      <  1:701200/151‑1 (MQ=255)
                                                                                                                                              ggAGTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCa      <  1:535976/151‑1 (MQ=255)
                                                                                                                                                aGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGc    <  1:361518/151‑1 (MQ=255)
                                                                                                                                                 gTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCg   >  3:268293/1‑151 (MQ=255)
                                                                                                                                                 gTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCg   >  3:165040/1‑151 (MQ=255)
                                                                                                                                                 gTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCg   >  3:202889/1‑151 (MQ=255)
                                                                                                                                                 gTTCCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCg   <  2:1115100/151‑1 (MQ=255)
                                                                                                                                                   ttCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCgg  <  1:884051/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGG  >  NC_000913/4381356‑4381652

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