Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,578,1430AC42.6% 91.4 / 13.8 55L79V (TTA→GTA) maeBmalic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (23/8);  new base C (0/23);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 7.37e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTG  >  NC_000913/2577997‑2578289
                                                                                                                                                         |                                                                                                                                                   
ttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑cccc                                                                                                                                                >  2:378765/1‑151 (MQ=255)
  cGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCa                                                                                                                                               <  1:691542/150‑1 (MQ=255)
   gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGc                                                                                                                                             <  2:721852/151‑1 (MQ=255)
   gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGc                                                                                                                                             <  2:297466/151‑1 (MQ=255)
        cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCAccgc                                                                                                                                        >  1:903849/1‑151 (MQ=255)
        cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCAccgc                                                                                                                                        <  2:859124/151‑1 (MQ=255)
          gcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCATGTTCAAATACATCAAGCCCGGCGAATTTCTT‑AAAC‑AGAC‑CGCCCTTCCCTTCCATCACGGC‑TTGG‑CCTG‑CCAGGGCGCCAATGTAGCCTA‑A‑CCCCTGCCccgccg                                                                                                                                      <  1:268987/151‑1 (MQ=255)
           cgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGt                                                                                                                                     >  2:58620/1‑151 (MQ=255)
            gACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGt                                                                                                                                     >  2:709781/1‑150 (MQ=255)
                     aaTAAATTTGCCGGGGTCGAGTCCGTCAACTTAAATGTAAAATACACCAACCCCGGCAATTTTCTA‑AAAC‑AGAC‑CCCCCTTCCCTCCCACCCCCGT‑TTTG‑CCTC‑CCAGCGCGCAAATGTGGCCTA‑C‑CCCCAGCACCGCCGTACCGTTa                                                                                                                              <  2:510795/148‑1 (MQ=255)
                      aTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTag                                                                                                                             >  1:878895/1‑148 (MQ=255)
                        aaaTTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTag                                                                                                                             >  2:21830/1‑146 (MQ=255)
                              gTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGcc                                                                                                                  >  1:92609/1‑151 (MQ=255)
                                  gggTGGTGTTCTCCAACTTGAATGCAAAATACGCAAATCCGGGAGAATTTTTA‑AAAC‑AAAA‑CCCCTTTGCTTCCCACCCCTGG‑TTTG‑CCTC‑CCGCCGC‑CCAACGTTGCCCA‑ACCCCCGGCTCCGCCGTACCGTTAGAGATCACCGccacc                                                                                                               <  1:468924/150‑1 (MQ=255)
                                     tCGAGTTCGTCACCTTAAATTTAAAATCCAACACCCCGGGCGAATTTCTA‑AAAC‑AGAC‑CGCCCTTCCCTCCCACCACGGT‑TTTG‑CCTG‑CCAGCGCCCAAATTTTCCCTA‑C‑CCCAAGCACCGCCGTACCGTTAGAGATCACCGccacca                                                                                                              <  1:1000105/148‑1 (MQ=255)
                                       gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                          >  1:206691/1‑150 (MQ=255)
                                         tttCGTCAATTTCAATGTAAAATACACCCCTCCCGGCGAATTTCTA‑AAAC‑AACC‑CCCCCTTCCTTTCCACCCCCGT‑TTTG‑CCGG‑CCCGCGCGCATATGTTGCCTA‑C‑CCCCGGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTa                                                                                                         <  1:134065/148‑1 (MQ=255)
                                                accTTCAATGCAAAAAACATCAACCCCGGCGAATTTCTA‑AAAC‑AGAC‑CGCCCTTCCCTTCCACCCCCGT‑TTTG‑CCTG‑CCGCCGCGCCAATTTTGCCAA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCggg                                                                                                  <  2:196247/147‑1 (MQ=255)
                                                  cTTAAATGTCAAATACATCACCCCGGGCAAATTT‑TT‑AAACAAAAA‑CCCCCTTCCCTCCCA‑CCCCGG‑TTTG‑GCCG‑CCCGGCCCCCCATTTTTCCAA‑ACCCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCg                                                                                                <  2:938859/149‑1 (MQ=255)
                                                   ttCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                              >  2:512162/1‑150 (MQ=255)
                                                      aaTGTCAAATACACCAACCCCGGCGAATTTCTA‑AAAC‑AGAC‑CGCCCTTCCCTTCCACCCCCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑C‑CCCAAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                            <  2:973122/149‑1 (MQ=255)
                                                        tGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                        >  2:917692/1‑151 (MQ=255)
                                                          tCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                                      >  2:555838/1‑151 (MQ=255)
                                                            aaaTACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                                      >  2:149942/1‑149 (MQ=255)
                                                             aaatacaccaccccGGGCAATTTCCT‑AAAA‑AAAACCCCCCTTCCCTCCCACCCCCGT‑TTTG‑CCTC‑CCGGCGCCCATATTTTCCCTA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                        <  1:530265/136‑1 (MQ=255)
                                                                       ccccGGCAATTTTCTA‑AAAC‑AG‑A‑CCCCCTTCCCCTCCCACCCCCG‑TTTGCCCCG‑CCCGCGC‑CCCATGTTCCCTACC‑CCCCCGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                          <  1:605280/149‑1 (MQ=255)
                                                                        cccGGAGATTTTTTA‑AAAC‑AAAC‑CCCCCTTCCCTCCCACCCCGGG‑TTGG‑CCGC‑CCAGACCCACAATGTTGCATA‑C‑CCCAGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc                                                                           <  1:379411/148‑1 (MQ=255)
                                                                         ccGGGGAATTTCTTAAAAA‑AGAC‑CCCCTTTCCCCCCCATCCCCGT‑TTTC‑CCGC‑CCGCCGCCCAAATTTGCCCTC‑C‑CCCCGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                               <  1:645972/144‑1 (MQ=255)
                                                                         ccGGCTAATTTATA‑AAAC‑AAAA‑CG‑CCTTTCCCTCCATCCCCCG‑TTTGCCCTC‑CCAGCGCGCAAATTTTGCCTC‑C‑CCCAAGACCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                              <  2:630061/144‑1 (MQ=255)
                                                                              gAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                                  >  2:1012915/1‑151 (MQ=255)
                                                                                 tttCTA‑AAAC‑AGAC‑CCCCTTTCCCTCCCCTCCCCGG‑TTTC‑CCTC‑CCAGCGCCCAAATTTGCCCTA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                <  1:123899/150‑1 (MQ=255)
                                                                                 tttCTA‑AAAA‑AGAC‑CCCCTTTCCCTTCC‑TCCCCGGTTTTC‑CCCC‑CCAGCCCCCCAATGTTGCCTA‑C‑CCCCGCCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                                <  1:108728/150‑1 (MQ=255)
                                                                                  ttCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                               >  1:541440/1‑150 (MQ=255)
                                                                                  ttCTA‑AAAC‑AGAC‑CCCCCTTGCCTTCAACCCCGGT‑TTTG‑CCTG‑CAAGCGCGCCATTGTTCAAAA‑C‑CCCCAGCTCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                               <  2:21537/150‑1 (MQ=255)
                                                                                   tCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                             >  1:435385/1‑151 (MQ=255)
                                                                                       aaaaC‑AGAA‑CCCCCTTGCTTCCCACCCCGGG‑TTTG‑CCTG‑CCGCCCCCCTATTGTGGCCTA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                          <  1:784297/150‑1 (MQ=255)
                                                                                          aC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                         >  2:480820/1‑149 (MQ=255)
                                                                                          aC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                         >  1:25728/1‑149 (MQ=255)
                                                                                             aGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGAc                                                     >  2:142614/1‑151 (MQ=255)
                                                                                                  cGCCCTTCCCTTCCACCCCGGG‑TTGG‑CCCG‑CAAGCCCGCAAATTTTCCCAA‑A‑CCCAAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                                 <  2:891095/151‑1 (MQ=255)
                                                                                                   gCCCTTGCCTTCCGCAACAGC‑ATTA‑TCTG‑CTATCACGCCAATCTGGCGGA‑T‑CCCCAGCACCGCCGTACCGTGAGAGACCTCCGCCACCAGGTTACCTCGGGCCGTATATCTGAAGGCTTTTAACGGGTCTTTTTCGACTTCAAGACAAgg                                                 >  1:492966/1‑150 (MQ=255)
                                                                                                   gCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                                >  2:144350/1‑151 (MQ=255)
                                                                                                    cctttcccTCCCACCCCCGT‑TTTG‑CCTG‑CCGGCGCGCAATTTTTCCCTA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGa                                                      <  1:336034/138‑1 (MQ=255)
                                                                                                       ttcccTCCCACCCCCGT‑TTTG‑CCTG‑CCGGCGCCCCATTGTTCCCTA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGc                                              <  1:225129/146‑1 (MQ=255)
                                                                                                         tccTTCCCACCCCGG‑TTTT‑CCCGCCCGGCGCCCAAATGTCCCCAA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                            <  2:981811/149‑1 (MQ=255)
                                                                                                              ccccccccGG‑TTGC‑CCTC‑CCGCCGCGCAATTTTTGCCAA‑C‑CCCCCGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                        <  1:822777/142‑1 (MQ=255)
                                                                                                                  ctccGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACACCTgg                                  >  1:164755/3‑150 (MQ=255)
                                                                                                                   ccccg‑TTTGCCCTC‑CCGGCCCGCCAATGTCGCATA‑C‑CCCCATCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACg                                       <  1:189997/141‑1 (MQ=255)
                                                                                                                      gg‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAg                              >  1:753034/1‑150 (MQ=255)
                                                                                                                      gg‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                             >  2:541738/1‑151 (MQ=255)
                                                                                                                                            cAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg          >  2:992638/1‑151 (MQ=255)
                                                                                                                                             attGTTGCCAA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAg                <  2:564744/142‑1 (MQ=255)
                                                                                                                                             attGTTCCCAA‑C‑CCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcg            <  2:656555/146‑1 (MQ=255)
                                                                                                                                                 ttGCCTA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCt        <  2:773341/148‑1 (MQ=255)
                                                                                                                                                    ccTA‑C‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTg  <  1:199827/151‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
TTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTT‑AAAC‑AGAA‑CGCCCTTGCCTTCCATCACCGG‑TTTG‑CCTG‑CCAGCGCGCCAATGTTGCCTA‑A‑CCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGTGTTG  >  NC_000913/2577997‑2578289

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