Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4960TA47.0% 63.5 / 40.9 66intergenic (+301/‑155)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base T (28/7);  new base A (0/31);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 2.29e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTG  >  NC_000913/29347‑29641
                                                                                                                                                     |                                                                                                                                                 
ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCACTCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                                   >  1:490727/1‑150 (MQ=255)
ttAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAatat                                                                                                                                                   <  2:44666/150‑1 (MQ=255)
 tAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAAtata                                                                                                                                                  <  2:952824/150‑1 (MQ=255)
       aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                           <  1:792110/151‑1 (MQ=255)
              gACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCa                                                                                                                                     >  1:123030/1‑150 (MQ=255)
                   tttGGTCCATTTTTTTCTGCTCTTTTTTATTTAATGCAATCTTTTTGCTGCGAAAGCTTTTTCAAAAACAGGTTAGATAATTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGGGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAg                                                                                                                                <  1:864406/150‑1 (MQ=255)
                         ccccTTTTTTCTGCTCTTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCTTTTCCCAGAACAGGTTAGATGATTTTTTTGTCGCTTAATGCCTGTAAAACATGCATGACCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACt                                                                                                                          <  2:293011/147‑1 (MQ=255)
                             ttttttCTGCTCTTTTTTTTTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAAAACAGGTTAGATGATTTTTTTGTCGCTAAATGCCTGAAAAACATGCATGAGCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTa                                                                                                                      <  2:1020817/150‑1 (MQ=255)
                             ttttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACtttt                                                                                                                       >  2:924510/1‑149 (MQ=255)
                              tttttCTGCTCGTTTTTTTTTCATGAATTCTTCTTGCTGCGCAAGCTTTTCCCAGAACAGGTTAGATGATCTTTTTGTCGTTTAATGCTTGTAAAACATGCATGACCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAg                                                                                                                     <  2:123716/150‑1 (MQ=255)
                                     gCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCat                                                                                                              >  1:808961/1‑150 (MQ=255)
                                      ctctTTTTTTTTTCATCCAATCTTCTTCCTGCGCAAGCTTTTTCCAAAACAGTTTAGATGATTTTTTTGTCGCTTAAGGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCa                                                                                                               <  1:821472/148‑1 (MQ=255)
                                      ctctTTTTTATTTCATGAATTCTTCTTGCTGCGCAAGCTTTTCCAAGAACAGTTTAGATAATTTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGccc                                                                                                                <  1:1007231/147‑1 (MQ=255)
                                      cTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                            >  1:489652/1‑151 (MQ=255)
                                        cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGTGCCCATAtc                                                                                                           >  1:681605/1‑150 (MQ=255)
                                          tttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                        >  1:438726/1‑151 (MQ=255)
                                           ttttttttCATGCATTCTTCTTGCTGCGAAAGCTTTTTCCAGAACAGTTTAAATGTTTTTTTTGCCGTTAAATGCCTGTAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                        <  1:542434/150‑1 (MQ=255)
                                           ttttATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTcc                                                                                                        >  2:144572/1‑150 (MQ=255)
                                             ttATTTCATGCAATTTTCTTGCTGCGCAAGCTTTTTCCAGAAAAGGTTAGATGATTTTTTTGTCGTTAAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAg                                                                                                      <  2:830303/150‑1 (MQ=255)
                                                        aaTCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                          >  1:921169/1‑151 (MQ=255)
                                                        aaTCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGtt                                                                                          >  1:1002845/1‑151 (MQ=255)
                                                         aTCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGttt                                                                                         >  2:673813/1‑151 (MQ=255)
                                                            ttcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                      >  1:51070/1‑151 (MQ=255)
                                                            ttcttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                      >  1:210309/1‑151 (MQ=255)
                                                             tcttGCTGCGCAAGCTTTTTCCAGAACAGGTTAGATGATTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                      <  2:898389/150‑1 (MQ=255)
                                                               ttGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGcc                                                                                      >  2:322814/1‑148 (MQ=255)
                                                                        aaGCTTTTTCAAGAACAGTTTAGATGTTTTTTTTGTCGTTTAATCCCTAAAAAACATGCAGGACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                          <  1:981133/151‑1 (MQ=255)
                                                                        aaGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtcg                                                                          >  2:483199/1‑151 (MQ=255)
                                                                           cGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGt                                                                        >  2:388464/1‑150 (MQ=255)
                                                                             ttttCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTaa                                                                      <  2:171506/150‑1 (MQ=255)
                                                                             ttttCCAAAACAGTTTAGATGTTCTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCgg                                                                         <  1:151253/147‑1 (MQ=255)
                                                                              tttCCAGAACAGGTTGGATGTTTTTTTTGTCGCTTAATGCCTGAAAAACATGCATGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                                     <  2:584482/150‑1 (MQ=255)
                                                                               ttcaaaaaaagtttAGATTATTTTTTTGTCTTTTAATGCCTGAAAAACATGCAGGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                    <  1:444071/138‑1 (MQ=255)
                                                                               ttCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGc                                                                    >  1:467579/1‑150 (MQ=255)
                                                                                   aaaaCAGGTTAGTTGTTTTTTTTGTCGCTTAATGCCTGAAAAACATGAATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                                 <  2:380854/147‑1 (MQ=255)
                                                                                   aGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                               >  2:627121/1‑151 (MQ=255)
                                                                                     aaaaGGTTAGATGATTTTTTTGTCGCTTAATCCCTGAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTg                                                             <  1:223850/148‑1 (MQ=255)
                                                                                      aCAGTTTAGATAATTTTTTTGTCGTTTAATCCCTGAAAAACATGCAGGACCAAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAttt                                                              <  1:588405/149‑1 (MQ=255)
                                                                                         ggTTAGATGATTTTTTTGTCGCTTAATGCCTGAAAAACATCCAGGAGCCACAAAAAAAAATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                         <  2:99546/151‑1 (MQ=255)
                                                                                         ggTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                                         >  1:943513/1‑151 (MQ=255)
                                                                                                tgatTTTTTTGTCTTTTAATGCCTGAAAAACATGCATGACCCACAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGa                                                       <  1:359307/146‑1 (MQ=255)
                                                                                                tgatTTTTTTGTCGCTAAATGCCTGAAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                                   <  2:175265/150‑1 (MQ=255)
                                                                                                 gatTTTTTTGTCGTTTAATCCCTAAAAAACATGCATGACCCAAAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                  <  1:663035/150‑1 (MQ=255)
                                                                                                   tttttttGTCGCTAAATGCCTGAAAAACATGCAGGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAtt                                                     <  2:566820/143‑1 (MQ=255)
                                                                                                   tCTTTTTGTCGCTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                  <  2:429706/148‑1 (MQ=255)
                                                                                                    cTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                               >  1:119033/1‑150 (MQ=255)
                                                                                                    cTTTTTGTAGCTTCCTGCATGTAAAACATGCATGAACCACAAACTAGTATAAAAAATCGCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATACCGCCGTTTGACAGAAATTCGGCGGTAAGCAGCTTTGCATTGATTTActc                                                >  1:656862/1‑147 (MQ=255)
                                                                                                     tttttGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtca                                              <  1:245602/150‑1 (MQ=255)
                                                                                                     tttttGTCGCTTAATGCCTTTAAAACATGCATGAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                                  <  2:500849/146‑1 (MQ=255)
                                                                                                     tttttGTCGATTAATGCCTGAAAAACATACTTGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                               <  1:689754/149‑1 (MQ=255)
                                                                                                        ttGTCGCTTAATCCCTAAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtcatca                                           <  1:804665/150‑1 (MQ=255)
                                                                                                          gTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                        >  1:956409/1‑151 (MQ=255)
                                                                                                           tCGCTTAATGCCTGAAAAAAATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtg                                        <  1:29290/150‑1 (MQ=255)
                                                                                                            cGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgtg                                      >  2:62050/1‑151 (MQ=255)
                                                                                                               ttAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                   >  1:636224/1‑151 (MQ=255)
                                                                                                                 aaTGCCTGAAAAACAGGCATGACCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAt                                   <  1:63567/149‑1 (MQ=255)
                                                                                                                  aTGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTa                                 >  2:965683/1‑150 (MQ=255)
                                                                                                                  aTGCCTGAAAAAAAGGCATGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTaa                                <  1:302758/151‑1 (MQ=255)
                                                                                                                   tGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                               >  1:856184/1‑151 (MQ=255)
                                                                                                                           aaaCAGCCATCAGCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt                       <  1:310147/151‑1 (MQ=255)
                                                                                                                              ctggaaagaccaaaaaaaTAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaa                     <  1:168624/136‑1 (MQ=255)
                                                                                                                                     gacccACAAAATAAGAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg             <  1:962359/148‑1 (MQ=255)
                                                                                                                                       cccAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAgtgagtga            <  1:227809/149‑1 (MQ=255)
                                                                                                                                          acaAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAatat        >  1:345553/1‑151 (MQ=255)
                                                                                                                                          aaaaaaaaATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGaa           <  1:905573/146‑1 (MQ=255)
                                                                                                                                          aaaaaaTAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAata         <  1:134632/148‑1 (MQ=255)
                                                                                                                                            aaaaTAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAata         <  2:427302/148‑1 (MQ=255)
                                                                                                                                             aaaTAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATtct     >  2:704783/1‑151 (MQ=255)
                                                                                                                                                taataTAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTg  >  1:203971/1‑151 (MQ=255)
                                                                                                                                                  atatAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTctct   <  2:721956/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTGAATATTCTCTG  >  NC_000913/29347‑29641

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