Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091329,4860CA43.1% 117.3 / 41.1 72intergenic (+291/‑165)dapB/carAdihydrodipicolinate reductase/carbamoyl phosphate synthetase small subunit, glutamine amidotransferase
Reads supporting (aligned to +/- strand):  ref base C (33/8);  new base A (0/31);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 4.14e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTG  >  NC_000913/29339‑29630
                                                                                                                                                      |                                                                                                                                                
ttgttgTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACaa                                                                                                                                                >  2:529708/1‑150 (MQ=255)
 tgttgtTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACaaa                                                                                                                                               <  1:457093/150‑1 (MQ=255)
           aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                     >  2:568536/1‑150 (MQ=255)
           aTGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaa                                                                                                                                     >  1:607076/1‑150 (MQ=255)
            tGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaa                                                                                                                                    >  2:23220/1‑150 (MQ=255)
              tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATaaaaaa                                                                                                                                  <  2:805173/150‑1 (MQ=255)
              tAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                 >  1:56576/1‑151 (MQ=255)
               aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAAt                                                                                                                                 <  1:312474/150‑1 (MQ=255)
               aaaTTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATc                                                                                                                                >  1:136879/1‑151 (MQ=255)
                    ttGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                           >  1:704319/1‑151 (MQ=255)
                     tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                           >  1:469746/1‑150 (MQ=255)
                     tGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGcc                                                                                                                           >  1:373771/1‑150 (MQ=255)
                         cATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                       >  2:110/1‑150 (MQ=255)
                          aTTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTa                                                                                                                       >  1:302015/1‑149 (MQ=255)
                           tttGTTCCTTTTTTTTCTGCTCGTTTTTATTCAATCCATTCTTCTTGCTGCGCAGGCTTTTCCCAAAACAGGTTAGATAATTTTTTTTTCTTTTAATCCCTAAAAAACATGCATGCGCCAAAAAAAAATATAAAAAACCCCGCCATTAAg                                                                                                                     <  2:739451/150‑1 (MQ=255)
                            ttGGTCCTCTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCTTTTTCAAAAACAGTTTAGATAATTTTTTTGTCGCTAAATGCCTGAAAAACATGCATGAGCCACAAAATAATATAAAAAACCCCGCCATTAAGt                                                                                                                    <  2:224001/150‑1 (MQ=255)
                             tGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTg                                                                                                                  >  1:745994/1‑151 (MQ=255)
                                tCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACt                                                                                                               >  1:526123/1‑151 (MQ=255)
                                  cACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACttt                                                                                                             >  1:672491/1‑151 (MQ=255)
                                       ttttCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAgcg                                                                                                        >  2:828440/1‑151 (MQ=255)
                                                cGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                               >  2:394327/1‑151 (MQ=255)
                                                 gTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtctc                                                                                              >  2:503641/1‑151 (MQ=255)
                                                  ttttttttTCATGCAATTTTCTTCTTGCGAAGGTTTTTTCAAAACCAGTTAAAATGTTTTTTTTGTCCTTTAATGCTTTTAAAACATGCATGACCCAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCatat                                                                                                 <  1:640971/147‑1 (MQ=255)
                                                  tttttATTTCATGAAATCTTCTGGCTGCGCAAGCTTTTCCCAGAACGGTTTAGATGATTTTTTTGTCGCTAAATGCCTGAAAAACTTGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATAtct                                                                                               <  1:720134/149‑1 (MQ=255)
                                                    tttATTTCATGAATATTTTTTGCTCGGCAATGTTTTTCAAGACCAGTTTAAATGATTTTTTGTCCGTTAAATCCTTTAAAAAAATGTATGCCCCAAAAAATAATATAAAACACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCa                                                                                            <  2:685976/150‑1 (MQ=255)
                                                      tATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGa                                                                                          >  1:22882/1‑150 (MQ=255)
                                                       aTTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGAATTTTAGCGCCCATATCTCCAGAAt                                                                                        >  2:566989/1‑151 (MQ=255)
                                                        tttCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                       >  2:85110/1‑151 (MQ=255)
                                                        tttCAGGCATTCTTCTTGCGGCGCAGGCTTTTCCCAGAAAAGGTAAGATGATTTTTTTGTCTTTAAAGCCCGGAAAAACAGGCAGGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAAt                                                                                        <  1:764687/150‑1 (MQ=255)
                                                            aTGCAATCTTTTTGCGGCGCAAGCTTTTTCAAGACAAGGTTAAATAATTTTTTTTTCGTTTAATCCCTAAAAAACATGCTTGAGCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccg                                                                                    <  2:601395/150‑1 (MQ=255)
                                                              gCATTCTTCTTGCTGCACAAGCTTTTTCAAGAACAGGTTAAATGATTTTTTTGTCTCTTAATGCCTGTAAAACATTCATGAGCAACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATg                                                                                       <  2:89663/145‑1 (MQ=255)
                                                                attttccttGCTGCGAAAGCTTTTCCAAAAACAGTTAAGATGATTTTTTTGTCGCTAAATCCCGGAAAAAAATGCATGACCCACAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgcc                                                                                     <  2:239023/139‑1 (MQ=255)
                                                                attcttctTGCTGCGCAAGCGTTTACAAAAACAGTTTAGATAATTTTTTGGTCGTTTAACGCCTGAAAAAAATTAAGGACCCACAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATgccgcc                                                                                  <  2:49172/146‑1 (MQ=255)
                                                                       ttGCTGCGCAAGCGTTTTCAAAAAAAGGTTAAATAATCTTTTTGTCTCTTAATGCCTTAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAg                                                                         <  2:83029/150‑1 (MQ=255)
                                                                       ttGCTGCACATGTTTTTTCCAGAGCAGTTAAGATTATTTTTTTTTCGCTAAATCCCTATAAAACAGGCATGCGCCAAAAAATAATAAAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCa                                                                          <  2:68209/149‑1 (MQ=255)
                                                                            gcaaaggctttttCAAGAACGGTTAAGATGATTTTTTTGTCGTTTAATGCCTAAAAAACATGCATGACCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGaa                                                                       <  2:577322/138‑1 (MQ=255)
                                                                             cgaaacctttttCAAAAACGGTTTAGATGATTTTTTTGCCTTTAAATCCTTAAAAAACATTCATGAGCAAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAAt                                                                     <  2:314747/140‑1 (MQ=255)
                                                                                aggctttttCAAAAACGGGTTAGATGTTTTTTTTGCCGTTAATTCCTTGAAAAACAGGTGTGACAAACAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                <  1:126702/145‑1 (MQ=255)
                                                                                aagttttttCCAAAACGGGTTAAATGATCTTTTTTTCGTTAATTGCATTAAAAACATGCATGTCCCAAAAAAAAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATtcgtc                                                                <  2:221437/145‑1 (MQ=255)
                                                                                     ttttCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTa                                                            >  1:652513/1‑149 (MQ=255)
                                                                                      tttCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAg                                                          >  2:48550/1‑150 (MQ=255)
                                                                                         caaGAACAGGTTAAATGATTTTTTTGTCTTTTAATGCTTGAAAAACATGCATGAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                      <  1:474660/149‑1 (MQ=255)
                                                                                           aGAACAGTTTAGATGATTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAtt                                                    <  1:220789/151‑1 (MQ=255)
                                                                                           aGAACAGTTTAGATGATTTTTTTGTCGCTTAATGCTTGAAAAACATGCAGGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGAt                                                     <  1:820699/150‑1 (MQ=255)
                                                                                            aaaaaGGTTAAGATGATTTTTTTTTCTTTTATTCCCTGAAAAAATTGCTTGACCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGa                                                      <  2:566108/145‑1 (MQ=255)
                                                                                              aCAGTTTAGATCATCTTTTTTTCGATTAATGCCTGAAAAACATGCATGGCCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCa                                                        <  2:136842/144‑1 (MQ=255)
                                                                                                      aGATTATTTTTTTTTCTCTTACTGCCTAAAAAACATGCATGAACCAAAAAATAATATAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCAtt                                              <  2:114054/147‑1 (MQ=255)
                                                                                                      aGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGCTTTGCATTGAtt                                          >  2:798088/1‑151 (MQ=255)
                                                                                                         tgttttttttGTCTTTTAATGCCTAAAAAACATGCAGGAGCCAAAAAATAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                        <  1:742625/145‑1 (MQ=255)
                                                                                                         tgatTTTTTTGTCGTTTAATGCCTGAAAAACATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTa                                        <  1:188717/150‑1 (MQ=255)
                                                                                                         tgatCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTAc                                       >  1:120667/1‑151 (MQ=255)
                                                                                                           atCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGt                                     >  1:786668/1‑151 (MQ=255)
                                                                                                            tCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGtc                                    >  2:457018/1‑151 (MQ=255)
                                                                                                              ttttttTCTCTTAATCCCTAAAAAAAATTAATGACCAAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGAttt                                         <  2:591716/144‑1 (MQ=255)
                                                                                                                 ttGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                               >  2:736392/1‑151 (MQ=255)
                                                                                                                 ttGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTcatcat                               >  1:438929/1‑151 (MQ=255)
                                                                                                                  tGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCAtt                              >  1:87830/1‑151 (MQ=255)
                                                                                                                    tCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATtgt                            >  1:257190/1‑151 (MQ=255)
                                                                                                                         tAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAAtt                       >  2:434008/1‑151 (MQ=255)
                                                                                                                                 ctaaaaaaCATTGATGACCCAAAAAAAAATAAAAAAAACCCCCGCATTTAATTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAat                    <  2:704096/144‑1 (MQ=255)
                                                                                                                                     aaacCATGCATGACCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAata                   <  1:538588/145‑1 (MQ=255)
                                                                                                                                     aaaaCTTGCATGACCAAAAAAAAAAAAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATg                 <  1:85139/147‑1 (MQ=255)
                                                                                                                                     aaaaCATGAATGGCCCAAAAAATAAAAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCaaa             <  1:200710/151‑1 (MQ=255)
                                                                                                                                      aaaCATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAAt            >  1:650173/1‑151 (MQ=255)
                                                                                                                                       aaCATGCATGACCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATa           <  1:799452/151‑1 (MQ=255)
                                                                                                                                          atgcatgACCCAAAAAAAAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa         <  2:541862/150‑1 (MQ=255)
                                                                                                                                            gaatgAGCCAAAAAAAAATAAAAAAAACCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt       <  1:182175/148‑1 (MQ=255)
                                                                                                                                             cttgAGCAAAAAAAAAATAAAAAAAACGCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagt       <  2:501314/147‑1 (MQ=255)
                                                                                                                                             catgAGCCAAAAAATAATAAAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa         <  1:303876/147‑1 (MQ=255)
                                                                                                                                              attggCCAAAAAAAAAAAAAAAAAACCCCGCCCTTAAGTTGGCTTTTAGCGGCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATaaa         <  1:791009/142‑1 (MQ=255)
                                                                                                                                              atgAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtg      >  1:641531/1‑149 (MQ=255)
                                                                                                                                                   ccAAAAAATAAAATAAAAATTCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAagtgagtg  <  2:25364/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
TTGTTGTTTTAATGTAAATTTTGACCATTTGGTCCACTTTTTTCTGCTCGTTTTTATTTCATGCAATCTTCTTGCTGCGCAAGCGTTTTCCAGAACAGGTTAGATGATCTTTTTGTCGCTTAATGCCTGTAAAACATGCATGAGCCACAAAATAATATAAAAAATCCCGCCATTAAGTTGACTTTTAGCGCCCATATCTCCAGAATGCCGCCGTTTGCCAGAAATTCGTCGGTAAGCAGATTTGCATTGATTTACGTCATCATTGTGAATTAATATGCAAATAAAGTGAGTG  >  NC_000913/29339‑29630

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