Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,578,1430AC50.6% 86.5 / 164.4 87L79V (TTA→GTA) maeBmalic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (38/5);  new base C (0/44);  total (38/49)
Fisher's exact test for biased strand distribution p-value = 1.45e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-15

GTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGC  >  NC_000913/2577996‑2578282
                                                                                                                                                         |                                                                                                                                            
gttCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAAccc                                                                                                                                          >  2:690216/1‑151 (MQ=255)
  tCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAAcccc                                                                                                                                         >  1:835340/1‑150 (MQ=255)
    gAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGc                                                                                                                                      >  2:497739/1‑151 (MQ=255)
     aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGc                                                                                                                                      >  2:708534/1‑150 (MQ=255)
     aGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCa                                                                                                                                     <  1:1222304/151‑1 (MQ=255)
      gcgcGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAc                                                                                                                                    >  2:325680/1‑151 (MQ=255)
       cgcgGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAcc                                                                                                                                   >  1:141312/1‑151 (MQ=255)
          ggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAccgcc                                                                                                                                <  2:163862/151‑1 (MQ=255)
                   ttaaaTAAATTTGTCGGGGTCGATTTCGTAAACTTAAATGTAAAATACATAAACCCCGGCGAATTTCTAAAACAGACCGCCTTTGCCTCCCATCCCCGTTTTGCCTGCCAGCGCGCCAATGTTGCCTACCCCCAGCACCGCCGTACCGtt                                                                                                                        <  2:1178609/147‑1 (MQ=255)
                     cAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGtt                                                                                                                        >  1:735020/1‑148 (MQ=255)
                       aTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTag                                                                                                                      >  1:468574/1‑148 (MQ=255)
                         aaaTTTGTCCGGGTCGAGTTCGTCAACTTCAATGTAAAATACATCAACCCCGGCAAATTTCTAAAACAGACCGCCCTTGCCTCCCATCCCCGTTTTGCCTGCCAGCGCGCCAATTTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGAt                                                                                                                  <  2:1017629/150‑1 (MQ=255)
                          attttGCCCGGGTCGGTTTCCACATCTACAATTAAAAAACCATAAACCCCGGCAATTTTCTAAAACAGACCCCCTTTGCCTCCCATCCCCGTTTTCCCTGCCGGCGCCCAAATTTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCa                                                                                                                <  1:403516/149‑1 (MQ=255)
                              tGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGc                                                                                                            >  2:770443/1‑151 (MQ=255)
                               gTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGcc                                                                                                           >  1:862386/1‑151 (MQ=255)
                                tCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGcc                                                                                                           >  1:1466683/1‑150 (MQ=255)
                                 ccGGGTCGATTTCGTCAACTTAAATGAAAAATACATAAACCCCGGATATTTTCTAAAACAGACCCCCCTTGCCTTCCACCCCCGTTTTCCCTCCCAGCGCGCAAATGTTGCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCg                                                                                                             <  1:1636001/147‑1 (MQ=255)
                                      tCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAgg                                                                                                     >  1:361464/1‑150 (MQ=255)
                                      tCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGt                                                                                                    >  2:355098/1‑151 (MQ=255)
                                        gttttCGTAAACTTAAATTAAAATTCCATCAACCCGGGCAATTTTTTAAAACAGACCCCCCTTCCCTCCCACCCCCGTTTTCCCGCCCGGCCCGCAATTTTTCCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAg                                                                                                      <  1:976344/144‑1 (MQ=255)
                                        gAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                   >  2:218383/1‑150 (MQ=255)
                                         aGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTAc                                                                                                 >  2:1173121/1‑151 (MQ=255)
                                           ttCGTCAACTTCAATGTAAAATCCACCACCCCGGGGAAATTTCTTAAAGAGACCGCCTTTGCCTTCCACCCCCGGTTTGCCTGCCAGCGCGCCAATGTTGCTTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                   <  2:1198172/147‑1 (MQ=255)
                                           ttCGTCAACTTAATTTAAAAATACAACACCCCGGGCAATTTTCTAAAAAAAACCCCCTTCCCCTCCCACCCCGGTTTTGCCTCCCCGCGCGCCATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                   <  1:1573943/147‑1 (MQ=255)
                                            tCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTc                                                                                              >  1:747893/1‑151 (MQ=255)
                                             cGTCACCTTCAATGTAAAAAACATCACCCCCGGCAATTTTTTAAAACAGACCCCCCTTCCCTCCCATCACCGGTTTCCCCCCCCGCGCCCAAATGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCg                                                                                             <  2:706544/151‑1 (MQ=255)
                                               tCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCggg                                                                                           >  2:953137/1‑151 (MQ=255)
                                                 aaCTTCAATGCAAAAACCACCACCCCCGGCAATTTTCTAAAACAGACCCCCCTTCCCTCCCACCCCCTTTTTCCCTCCCAGCCCCCCATTTTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                          <  1:1229692/150‑1 (MQ=255)
                                                  aCTTCAAATGTCAAAAACATCAACCCCGGCGATTTTCTTAAACAGACCCCCCTTCCTTCCCACCCCCGTTTTCCCTCCCAGCGCCCCATTGTTCCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                          <  2:1109339/150‑1 (MQ=255)
                                                   cTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                       >  2:730737/1‑151 (MQ=255)
                                                    ttCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                       >  1:165886/1‑150 (MQ=255)
                                                    ttCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                       >  2:345862/1‑150 (MQ=255)
                                                      cAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                     >  1:1568090/1‑150 (MQ=255)
                                                       aaaTGTAAAACACACAACCCCCGCGGATTTTCTAAAAAAAACCCCCCTTCCCTCCCCCCCCTGTTTCCCCGGCCAGCCCCCCATTTTGCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGt                                                                                       <  1:1241636/147‑1 (MQ=255)
                                                          tGTCAAATACATCAACCCCGGCAATTTTCTAAAACAGACCCCCCTTCCCTCCCATCCCCGTTTGCCCTGCCAGCGCCCAATTTTTCCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                   <  2:1405063/149‑1 (MQ=255)
                                                           gtaaaaaccATAAACCCCGGCGATTTTCTAAAACAGACCCCCTTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCGCCAATTTTGCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAttt                                                                                 <  2:793608/142‑1 (MQ=255)
                                                              aaaTCCACCAACCCCGGCAATTTTCTAAAACAGACCCCCCTTCCCTCCCACCCCGGGTTTGCCGCCCAGCGCCCAATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGt                                                                               <  2:523664/149‑1 (MQ=255)
                                                              aaaTCCACCAACCCCGGCAAATTTTTAAAACAAACCCCCCTTCCCTTCCACCCCGGTTTTGCCTCCCAGCGCGCAAATGTTCCCTACCCCCGGAACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg                                                                                <  2:607315/148‑1 (MQ=255)
                                                              aaaTACACCACCCCCGGCAATTTTTTAAAAAAGACCCCCTTTCCCTCCCACCCCGGTTTTCCCTCCCAGCGCGCAAATGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAg                                                                             <  2:1493620/151‑1 (MQ=255)
                                                              aaaCACATCACCCCCGGCAATTTTCTAAAAAAGACCCCCTTTCCCTCCCACCCCGGTTTTCCTGGCCACCGCCCCATTGTTCCCACCCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAg                                                                             <  1:72408/151‑1 (MQ=255)
                                                                   cATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtt                                                                         >  1:41823/1‑150 (MQ=255)
                                                                    aTCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                       >  1:1338819/1‑151 (MQ=255)
                                                                    aTCAACCCCGGCAATTTTCTAAAACAAACCCCCCTTCCCTTCCACCCCGGGTTTGCCTGCCAGCGCCCAAATGTCGCCAACCCCAAGAACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGc                                                                           <  2:1143383/147‑1 (MQ=255)
                                                                     tCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                       >  1:600769/1‑150 (MQ=255)
                                                                      cAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTa                                                                      >  1:1392312/1‑150 (MQ=255)
                                                                      cAACCCCGGCGAATTTCTAAAACAGAACGCCCTTGCCTCCCACCCCGGGTTTGCCTGCCAGCGCCCAAATGTGGCCTACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                       <  2:835402/149‑1 (MQ=255)
                                                                         tCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACgg                                                                  >  2:1705670/1‑151 (MQ=255)
                                                                              ggcaaTTTTTTAAAACAGACCCCCCTTCCCTCCCACCCCGGTTTTCCCTCCCCGCGCCCAATTGTTCCCTACCCCCACCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCt                                                              <  2:1335168/148‑1 (MQ=255)
                                                                                 aaTTTCTTAAACAAACCGCCTTCCCCTTCCATCCCGGTTTTCCCTGCCGGCGCCCCATTTTGGCTTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtttt                                                           <  2:621003/150‑1 (MQ=255)
                                                                                    tttttaaaacaaacCCCCCTTCCCTCCCACCCCCGTTTTGCCGGCCAGCGCCCAATTGTTCCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                         <  2:1145841/144‑1 (MQ=255)
                                                                                    ttCTAAACCAGACCCCCCTTCCCTCCCCCCCCGGGTTGCCCTCCCGGCGCCCAAATTTGCCCAACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                          <  1:1676774/148‑1 (MQ=255)
                                                                                     ttctaaacaagaccccccTTGCCTTCCACCCCCGGTTGCCCTCCCGGCGCCCAAGTTTCCCCTACCCCAAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                        <  2:534667/135‑1 (MQ=255)
                                                                                        tAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                   >  2:1569982/1‑151 (MQ=255)
                                                                                         aaaacaaacCCCCTTTCCCTCCCATCCCCGTTTTGCCGCCCAGCGCCCAATTGTTGCTAACCCCCAGCGCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                   <  1:800685/150‑1 (MQ=255)
                                                                                         aaaaCAGACCCCCCTTCCCTCCCCTCCCGGGTTTGCCTCCCGGCGCCCCAATTTTCCACACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTc                                                   <  1:1052549/150‑1 (MQ=255)
                                                                                          aaacaaacCCCCCTTCCCTCCCACCCCCGGTTTGCCTGCCAGCCCCCCAATGTTCCCAACCCCCACCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt                                                    <  2:696944/149‑1 (MQ=255)
                                                                                          aaaaaaaccccccTTGCCTCCCACCCCCGGTTTCCCGCCCAGCGCCCAAATGTTCCCTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                      <  2:165895/137‑1 (MQ=255)
                                                                                          aaaCAGACCCCCTTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCGCCAATGTTCCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                      <  2:945541/147‑1 (MQ=255)
                                                                                          aaaCAGACCCCCCTTGCCTCCCACCCCCGTTTGGCCGGCCAGCGCCCCAATTTTGCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt                                                    <  2:1514537/149‑1 (MQ=255)
                                                                                          aaaCAGACCCCCCTTCCCTCCCACCCCCGTTTTGCCGCCCAGCGCCCAAATGTTCCCAACCCCAAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                      <  1:199277/147‑1 (MQ=255)
                                                                                          aaaCAGACCCCCCTTCCCTCCCACCCCCGTTTTCCCTGCCAGCGCCCAATTTTTGCCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAt                                                      <  2:667878/147‑1 (MQ=255)
                                                                                          aaaCAGACCCCCTTTCCCTCCCTCCTCGGTTTTCCCGCCAAGCGCGCAATTTTCCCCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAtt                                                     <  1:1404783/148‑1 (MQ=255)
                                                                                                  cccctttcccTCCCACCCCCGGTTTCCCTCCCAGCCCCCAAATGTTCCCAACCCCAACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                          <  1:729565/143‑1 (MQ=255)
                                                                                                  cccccTTGCCTCCCACCCCCGTTTTCCCTCCCAGCCCGCCAGTTTTGCCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCa                                                  <  1:323693/141‑1 (MQ=255)
                                                                                                  cGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                          >  1:313492/1‑151 (MQ=255)
                                                                                                   gCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                         >  2:1286831/1‑151 (MQ=255)
                                                                                                     ccttcccttcCATCCCCGTTTGGCCTGCCAGCGCCCAATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGc                                       <  2:1679460/151‑1 (MQ=255)
                                                                                                         cccTCCCACCCCCGTTTTCCCTGCCAGCGCCCAATTGTTCCCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAgg                                          <  1:770380/143‑1 (MQ=255)
                                                                                                          cctcccaccccgggTTTCCCCGCCGGCGCGCAATTGTTCCCTACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCa                                   <  1:323921/138‑1 (MQ=255)
                                                                                                          ccTCCCACCCCCGTTTTGCCGGCCAGCGCGCAAATGTTCCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCg                                      <  2:1418268/147‑1 (MQ=255)
                                                                                                           cTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                 >  2:1320786/1‑151 (MQ=255)
                                                                                                              ccacccccGTTTTCCCTCCCAGCGCCCAAATGTTCCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGc                                       <  2:718639/136‑1 (MQ=255)
                                                                                                              ccacccccGGTTTCCCTGCCAGCGCGCAATTTTTGCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                 <  1:923000/142‑1 (MQ=255)
                                                                                                                      tgttttGCGAGCCCGCGCACCATTGTTCCCCCACCCGGGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg                         <  1:1108578/146‑1 (MQ=255)
                                                                                                                      ggTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAg                       >  1:1428548/1‑150 (MQ=255)
                                                                                                                       gTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGt                      >  2:479073/1‑150 (MQ=255)
                                                                                                                         tttGCCTCCCAGCGCGCAATTTTTCCCAACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGa                        <  2:1415812/147‑1 (MQ=255)
                                                                                                                           tgcctgccAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGc                  >  1:931218/1‑151 (MQ=255)
                                                                                                                              cctcccAGCGCCCCATTTTTCCCTACCCCCAGCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGa                        <  1:215862/143‑1 (MQ=255)
                                                                                                                               ctgccAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGcc                 >  1:356873/1‑149 (MQ=255)
                                                                                                                                    cAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc           >  2:807249/1‑151 (MQ=255)
                                                                                                                                      aGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc           >  2:530937/1‑150 (MQ=255)
                                                                                                                                      aGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc           >  1:583127/1‑150 (MQ=255)
                                                                                                                                       gcgcgcAATTTTTCCCAACCCCCACCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCa          <  1:702519/150‑1 (MQ=255)
                                                                                                                                          cccaaaTTTTCCCTACCCCCACCCCCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGAt       <  2:1308845/146‑1 (MQ=255)
                                                                                                                                              aattctcccTACCCCCACCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATGTTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGa        <  1:1668002/138‑1 (MQ=255)
                                                                                                                                              aaTGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCgcgc  >  1:1579098/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                            
GTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGC  >  NC_000913/2577996‑2578282

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: