Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,578,1430AC46.6% 90.5 / 17.9 58L79V (TTA→GTA) maeBmalic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (26/5);  new base C (0/27);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 8.16e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCG  >  NC_000913/2577994‑2578281
                                                                                                                                                      |                                                                                                                                          
tggttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTaa                                                                                                                                            <  2:299566/150‑1 (MQ=255)
  gttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAAccc                                                                                                                                         >  1:190456/1‑151 (MQ=255)
     cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAAcccc                                                                                                                                        >  1:430165/1‑149 (MQ=255)
          gcggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCAccg                                                                                                                                 <  1:842813/151‑1 (MQ=255)
            ggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCAccgcc                                                                                                                               <  1:501182/151‑1 (MQ=255)
            ggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCAccgcc                                                                                                                               <  1:462623/151‑1 (MQ=255)
              cgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGt                                                                                                                             >  1:390964/1‑151 (MQ=255)
               gACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGt                                                                                                                             >  2:393902/1‑150 (MQ=255)
                aCAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTAATTGTCAAATCCACCAACCCCGGCGATTTTCTTAAACAAAGCGCCCTGGCCTTCCATCCCCGATTTGCCTCCCAGCGCCCCAATGTTG‑CCAACCCCCACCACCGCCGt                                                                                                                             <  2:492083/149‑1 (MQ=255)
                        aaTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTag                                                                                                                     >  1:4632/1‑149 (MQ=255)
                           aaaTTTGTCCGGGTCGAGTTCGTCAACTTAAATGTCAAATACATCAACCCCGGCGATTTTCTTAAACAGACCGCCCTTGCCTCCCACCCCCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATc                                                                                                                <  1:647896/151‑1 (MQ=255)
                                     gggTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGccacc                                                                                                       >  1:457409/1‑150 (MQ=255)
                                          gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAgg                                                                                                    >  2:344002/1‑148 (MQ=255)
                                          gAGTTCCTCAACTTAAATGTAAAATACATCAACCCCGGCGATTTTCTTAAACAGACCCCCTTTGCCTTCCACCCCCGTTTTCCCTCCCGGCGCGCAAATGTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTa                                                                                                 <  1:413138/151‑1 (MQ=255)
                                           atttCGTCAACTTAAATGTCAAATACACCACCCCCGGCAAATTTCTTAAACAGACCGCCCTTCCCTCCCATCACCGGTTTGCCTCCCAGCGCCCCAATGTCG‑CCAACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGt                                                                                                   <  2:296197/146‑1 (MQ=255)
                                             ttCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCt                                                                                              >  2:360265/1‑151 (MQ=255)
                                                   aaCTTCAATGTAAAATACACCACCCCCGGCGAATTTCTTAAACAGAACCCCCTTCCCTTCCATCGCCGTTTTCCCTGCCAGCGCGCAATTGTTC‑CCAACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCgg                                                                                           <  1:612656/148‑1 (MQ=255)
                                                   aaCTTAAATGTCAAATACATCACCCCCGGCAATTTTCTAAAACAGACCCCCCTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCGCAATTGTCC‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTc                                                                                             <  2:216270/146‑1 (MQ=255)
                                                   aaCTCCAATGTAAAATACATCAACCCCGGCGATTTTCTAAAACAGACCCCCCTTCCCTCCCACCCCCGGTTGGCCTGCCAGCGCGCAAATTTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                         <  1:443110/150‑1 (MQ=255)
                                                    aCTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCgg                                                                                       >  2:46520/1‑151 (MQ=255)
                                                    aCTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCgg                                                                                       >  2:783868/1‑151 (MQ=255)
                                                      ttCAATGTAAAATACATCAACCCCGGCGATTTTCTTAAACAGACCCCCCTTGCCTCCCATCCCCGGTTTGCCTGCCGGCGCGCCAATGTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCgg                                                                                       <  1:247102/149‑1 (MQ=255)
                                                       tCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGTCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                    >  2:435142/1‑151 (MQ=255)
                                                        cAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                    >  2:37314/1‑150 (MQ=255)
                                                         aaTGTAAAATACATCAACCCCGGCAAATTTCTTAAACAGACCGCCCTTGCTTTCCATCCCCGGTTTGCCTGCCAGCGCGCAAATTTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                  <  1:574817/151‑1 (MQ=255)
                                                         aaTGTAAAATACACCAACCCCGGCGATTTTCTAAAACAGACCCCCCTTCCCTCCCACCCCCGGTTTCCCTGCCAGCGCCCAAATGTGC‑CCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                  <  2:405477/151‑1 (MQ=255)
                                                            gtaaaaTCCATCAACCCCGGCGATTTTCTAAAACAGACCGCCCTTGCCTTCCACCCCCGTTTTGCCTGCCAGCGCCCCATTGTTC‑CCTACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAtt                                                                                 <  2:136339/146‑1 (MQ=255)
                                                               aaaaccacaaaccccGGCAATTTTCTAAAACAGACCCCCCTTGCCTCCCACCCCCGGTTGGCCGGCCGGCCCCCAAATGTTC‑CCAACCCCAAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAg                                                                            <  2:35073/139‑1 (MQ=255)
                                                               aaaTCCATCATCCCGGGCAAATTCCTTAAACAGACCCCCCTTCCCTTCCACCCCCGGTTGGCCGGCCAG‑GCCCCCATGGTGCCCTACCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg                                                                               <  2:457474/148‑1 (MQ=255)
                                                                  tACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCt                                                                         >  1:379739/1‑151 (MQ=255)
                                                                     accAACCCCGGCAAATTTCTAAAAAAAAACGCCTTTGCCTCCAACCACGGTTTGGCCGGCAAGCGCGCAAATGTGG‑CCAAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGc                                                                          <  2:485573/145‑1 (MQ=255)
                                                                        aaTCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc                                                                   >  1:953759/1‑151 (MQ=255)
                                                                        aaTCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc                                                                   >  2:545778/1‑151 (MQ=255)
                                                                        aaTCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc                                                                   >  2:798794/1‑151 (MQ=255)
                                                                        aaTCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc                                                                   >  2:945261/1‑151 (MQ=255)
                                                                          ccccGGCGATTTTTTTAAACAGACCGCCCTTGCCTCCCACCCCGGGTTTGCCTGCCAGCGCCCAAATTTTC‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACgg                                                                 <  1:736454/150‑1 (MQ=255)
                                                                           cccGGCGATTTTCTTAAACAGACCGCCCTTGCCTTCCATCCCCGTTTTGCCTGCCAGCGCGCAAATTTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                  <  1:853691/149‑1 (MQ=255)
                                                                             cGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTc                                                              >  2:156324/1‑151 (MQ=255)
                                                                                  aaTTTCTTAAACAGACCGCCCTTGCCTTCCATCCCCGGTTTGCCTCCCGGCGCGCAAATTTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt                                                            <  1:535918/148‑1 (MQ=255)
                                                                                          aaaaAGACCGCCTTTCCTTCCCACCCCCGGTTTGCCTCCCAGCGCCCAATTTTTC‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                  <  2:918872/150‑1 (MQ=255)
                                                                                          aaaCAGACCCCCCTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCGCCAATGTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                     <  1:838842/147‑1 (MQ=255)
                                                                                          aaaCAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                  >  1:802035/1‑150 (MQ=255)
                                                                                           aaCAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                 >  2:560428/1‑150 (MQ=255)
                                                                                            aCAGAACGCCCTTGCCTCCCACCCCCGGTTTGCCTGCCAGCGCGCAATTGTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAg                                               <  2:404275/151‑1 (MQ=255)
                                                                                                accccccttcccttCCACCCCCGTTTTGCCTGCCAGCGCCCCAATTTGC‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                            <  1:584060/146‑1 (MQ=255)
                                                                                                aaCGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                            >  2:923546/1‑150 (MQ=255)
                                                                                                    cccTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTg                                       >  1:563099/1‑151 (MQ=255)
                                                                                                          ccTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCa                                  >  1:794780/1‑150 (MQ=255)
                                                                                                          ccTCCCACCCCCGGTTGCCCTCCCGGCGCCCCAATTTTC‑CCAACCCCAAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCg                                     <  1:203090/147‑1 (MQ=255)
                                                                                                                accccggttttCCCTGCCAGCGCCCCAATGTTC‑CCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                <  2:913083/138‑1 (MQ=255)
                                                                                                                      ggTTTGCCTCCCAGCGCGCAAATTTTC‑CCTACCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGa                       <  1:977829/149‑1 (MQ=255)
                                                                                                                      ggTTTGCCGGCCGGCGCGCAAATGTTG‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGa                       <  1:928794/149‑1 (MQ=255)
                                                                                                                           gcctgccAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc                >  1:723974/1‑151 (MQ=255)
                                                                                                                                ccAGCGCCCCAATGTTG‑CCTACCCCCGGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAgcgc           <  1:300138/151‑1 (MQ=255)
                                                                                                                                        ccAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgc   >  1:376211/1‑151 (MQ=255)
                                                                                                                                        ccAATGTTC‑CCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGa       <  1:154469/147‑1 (MQ=255)
                                                                                                                                         cAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg  >  1:114630/1‑151 (MQ=255)
                                                                                                                                         cAATGTCC‑CCTAACCCAAACGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgc   <  2:953643/150‑1 (MQ=255)
                                                                                                                                               ttcccTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg  <  2:821346/143‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                          
TGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTG‑CCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCG  >  NC_000913/2577994‑2578281

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