Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,418 A→T 5.0% intergenic (‑436/+167) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,4180AT 5.0% 796.5 / 23.7 230intergenic (‑436/+167)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (109/106);  new base T (6/6);  total (118/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 236 total aligned reads displayed.

TATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGG  >  NC_000913/987268‑987567
                                                                                                                                                      |                                                                                                                                                     
 atGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCaa                                                                                                                                                       <  2:2013514/150‑1 (MQ=255)
     tttCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAg                                                                                                                                                   >  1:5340202/1‑150 (MQ=255)
     tttCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCAAAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAg                                                                                                                                                   <  1:4635347/150‑1 (MQ=255)
      ttCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGa                                                                                                                                                  >  1:3453717/1‑150 (MQ=255)
        ccTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAAt                                                                                                                                                <  2:5784087/150‑1 (MQ=255)
         cTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATg                                                                                                                                               <  2:313538/150‑1 (MQ=255)
         cTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATg                                                                                                                                               >  2:3558933/1‑150 (MQ=255)
          ttATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAATATCTTGTTTCAGTCATCCAGAATGCCATAATTTAGTATGAAAATTGTGCTATTTGCTTAGTAAATTCTATTCCGGTATTAGAATgg                                                                                                                                              >  2:177430/1‑149 (MQ=255)
          ttATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGc                                                                                                                                              <  1:843797/150‑1 (MQ=255)
            atatTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCgg                                                                                                                                            <  1:3765456/150‑1 (MQ=255)
               tttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAtt                                                                                                                                         <  2:5686876/150‑1 (MQ=255)
               tttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAtt                                                                                                                                         >  2:178280/1‑150 (MQ=255)
               tttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAtt                                                                                                                                         >  1:3580420/1‑150 (MQ=255)
                ttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTa                                                                                                                                        <  2:3050628/150‑1 (MQ=255)
                ttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTa                                                                                                                                        <  1:965142/150‑1 (MQ=255)
                ttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTa                                                                                                                                        >  2:9236/1‑150 (MQ=255)
                ttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTa                                                                                                                                        <  2:5112437/150‑1 (MQ=255)
                ttAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTa                                                                                                                                        <  2:1956064/150‑1 (MQ=255)
                  aaaTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGt                                                                                                                                      <  1:2802829/150‑1 (MQ=255)
                  aaaTTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTACAGACATCCAGAATGCAAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGt                                                                                                                                      <  1:2659329/150‑1 (MQ=255)
                   aaTTAACTACCCGGAATTTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTCTGAAAATTGTGCTATTAGCTAAGAAATAACTAATCCGCAATAAGACTGCGGATTAGtt                                                                                                                                     <  2:3774096/150‑1 (MQ=255)
                     ttAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGtttt                                                                                                                                   >  1:73689/1‑150 (MQ=255)
                       aaCTAAACGGAAAGTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGCGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTATTGCGGATTAGtttttt                                                                                                                                 <  2:4737055/150‑2 (MQ=38)
                           aaaCGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGcgcg                                                                                                                             >  1:2170494/1‑150 (MQ=255)
                            aaCGGAATTTTTGTGTGTGCTGGAAACTTTAACGACCTTGCACAAAATCTTCTTTCTGACATCAATAATGCCAAAAATTAGTATTAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCATTCTTTTTGCGGTTTAGttttttcttag                                                                                                                            <  1:3904838/150‑7 (MQ=19)
                            aaCGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTGTTTGCGCGt                                                                                                                            >  2:4165982/1‑150 (MQ=255)
                               ggAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGcgcggtct                                                                                                                         >  1:5009841/1‑146 (MQ=255)
                               ggAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTaaa                                                                                                                         <  2:470090/150‑1 (MQ=255)
                               ggAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAACAGTTCATGcgcgactc                                                                                                                         >  2:1867015/1‑146 (MQ=255)
                                gAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAAt                                                                                                                        >  2:5101506/1‑150 (MQ=255)
                                    ttttGTTTCTGATTTAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTg                                                                                                                    <  1:1295655/150‑1 (MQ=255)
                                    ttttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCTGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTg                                                                                                                    <  2:3382788/150‑1 (MQ=255)
                                    ttttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTg                                                                                                                    <  1:1010205/150‑1 (MQ=255)
                                      ttGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAGTAGTTTTTGCGCGTAAATGTTGCa                                                                                                                  >  2:4228820/1‑150 (MQ=255)
                                         tttCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATg                                                                                                               <  2:3241395/150‑1 (MQ=255)
                                           tCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAACAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTa                                                                                                             <  2:2423022/150‑1 (MQ=255)
                                           tCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTa                                                                                                             >  1:3684199/1‑150 (MQ=255)
                                           tCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTa                                                                                                             >  1:761646/1‑150 (MQ=255)
                                           tCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTa                                                                                                             <  2:946514/150‑1 (MQ=255)
                                            cTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAATTAGTATGAAAATTGTGCTATTAGCTAAGAAATAACTAATCAGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTaa                                                                                                            >  2:4864977/1‑150 (MQ=255)
                                                 ggAAACTTTATCGACCTGGCACAAAATCTTCTTTCATACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTc                                                                                                       <  1:5040382/150‑1 (MQ=255)
                                                  gAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCACTAAGAATGCGGATTAGTTTTTGCGCGTCAATGTTGCAATGTAATATTcg                                                                                                      >  2:5009677/1‑149 (MQ=255)
                                                  gAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTcc                                                                                                      <  2:2972445/150‑1 (MQ=255)
                                                   aaaCTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAATATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCg                                                                                                     <  2:410033/150‑1 (MQ=255)
                                                      cTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTaa                                                                                                  <  1:5272210/150‑1 (MQ=255)
                                                       tttATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAAc                                                                                                 <  1:5738806/150‑1 (MQ=255)
                                                       tttATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAAc                                                                                                 >  1:4427934/1‑150 (MQ=255)
                                                       tttATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAAc                                                                                                 >  2:2204511/1‑150 (MQ=255)
                                                        ttATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACa                                                                                                <  1:1772968/150‑1 (MQ=255)
                                                          aTCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAgg                                                                                              >  2:2762895/1‑150 (MQ=255)
                                                             gACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGgatg                                                                                           >  1:2862749/1‑150 (MQ=255)
                                                              aCCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCCATGTAATATTCCGTAACAGgatga                                                                                          >  1:5536219/1‑150 (MQ=255)
                                                                cTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTTCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATc                                                                                        <  2:643400/150‑1 (MQ=255)
                                                                 tGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCg                                                                                       >  2:3958708/1‑150 (MQ=255)
                                                                     acaAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTAt                                                                                   >  1:1363872/1‑150 (MQ=255)
                                                                      caAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATg                                                                                  >  1:5702207/1‑150 (MQ=255)
                                                                           tcttctTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTCATATTCCGTAACAGGATGATCGTGAGGGCggg                                                                             >  1:4497359/1‑150 (MQ=255)
                                                                                tttCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGAGGATTAGTTTTTTAGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGc                                                                        >  2:4398771/1‑150 (MQ=255)
                                                                                  tCAGACATACAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaa                                                                      <  1:4931093/150‑1 (MQ=255)
                                                                                    aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAAt                                                                    <  2:460456/150‑1 (MQ=255)
                                                                                    aGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGAGTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCaact                                                                    >  2:1336264/1‑148 (MQ=255)
                                                                                           cAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTACCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGc                                                             >  1:4150779/1‑150 (MQ=255)
                                                                                             gactGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCtt                                                           <  2:104668/147‑1 (MQ=255)
                                                                                             gAATGCCAAAAAATAGTAAGAAAATAGAGCTATTACCTAACAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGATGCAAATAGATTGCtt                                                           <  1:2957982/150‑1 (MQ=255)
                                                                                               aTGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGa                                                         <  1:208259/150‑1 (MQ=255)
                                                                                                 gCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACt                                                       <  2:4925308/150‑1 (MQ=255)
                                                                                                  ccAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACtt                                                      >  2:468828/1‑150 (MQ=255)
                                                                                                   cAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTa                                                     >  2:1422518/1‑150 (MQ=255)
                                                                                                   cAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTa                                                     >  1:2992908/1‑150 (MQ=255)
                                                                                                   cAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTa                                                     <  2:2332088/150‑1 (MQ=255)
                                                                                                    aaaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaa                                                    <  1:3231158/150‑1 (MQ=255)
                                                                                                     aaaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaac                                                   <  1:3729633/150‑1 (MQ=255)
                                                                                                      aaaaTAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaaca                                                  >  1:1221144/1‑150 (MQ=255)
                                                                                                          tAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTaacaaaca                                              <  1:4271048/150‑1 (MQ=255)
                                                                                                             tATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGAAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATATATTGCTTGACTAAACAAACAGAt                                           >  1:1569447/1‑150 (MQ=255)
                                                                                                                 aaaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGGTGCAAATAGATTGCTTGACTTAACAAACAGATAACt                                       >  2:4800594/1‑150 (MQ=255)
                                                                                                                  aaaTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACtt                                      >  2:388131/1‑150 (MQ=255)
                                                                                                                    aTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGa                                    <  1:615403/150‑1 (MQ=255)
                                                                                                                    aTTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGa                                    <  1:5101506/150‑1 (MQ=255)
                                                                                                                           tATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGacaa                             >  2:4019547/1‑150 (MQ=255)
                                                                                                                           tATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGaaaa                             <  2:3567237/150‑1 (MQ=255)
                                                                                                                                gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCTATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaa                        <  2:2537575/150‑1 (MQ=255)
                                                                                                                                gCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaa                        <  2:73689/150‑1 (MQ=255)
                                                                                                                                  tAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaa                      >  1:4584773/1‑150 (MQ=255)
                                                                                                                                  tAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaa                      >  2:1279060/1‑150 (MQ=255)
                                                                                                                                   aaGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaa                     <  1:4300278/150‑1 (MQ=255)
                                                                                                                                    aGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaaa                    <  2:2480574/150‑1 (MQ=255)
                                                                                                                                     gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTTCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaaaa                   <  1:178200/150‑1 (MQ=255)
                                                                                                                                     gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaaaa                   <  2:1461197/150‑1 (MQ=255)
                                                                                                                                     gAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATaaaaaaa                   <  1:178280/150‑1 (MQ=255)
                                                                                                                                        aaaaCTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAACGATAAAAAAagag                >  2:5490131/1‑150 (MQ=255)
                                                                                                                                             tAATCCGCATTCTTATTGCGGATTAGTTTTTTCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGc           >  1:1451514/1‑150 (MQ=255)
                                                                                                                                             tAATCCGCATTCTTATTGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCACATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGc           >  1:1918987/1‑150 (MQ=255)
                                                                                                                                             tAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTg                                                          >  1:4289124/1‑103 (MQ=255)
                                                                                                                                             tAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTg                                                          <  2:4289124/103‑1 (MQ=255)
                                                                                                                                                tCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCggg        >  2:2938096/1‑150 (MQ=255)
                                                                                                                                                tCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCggg        >  2:1661040/1‑150 (MQ=255)
                                                                                                                                                   gCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcgc     >  2:547560/1‑150 (MQ=255)
                                                                                                                                                     aaTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTg   >  2:5365306/1‑150 (MQ=255)
                                                                                                                                                      aTAAGATTGCGGATTAGTTTTTGCTTTTAATTGTTGATTTTTATTTTTCCGTAACTGGTTTCTGGTTTTGGCGGGTTTGAAAATTTTTTGCTTGAGTTAACAAACTCATAACTTTACTTAAAATATAAAAAAAGAGCCAGCGGGCGCTgg  <  1:4062167/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TATGATTTCCTTATATTTAAATTAACTAAACGGAAATTTTGTTTCTGATGGAAACTTTATCGACCTGGCACAAAATCTTCTTTCAGACATCCAGAATGCCAAAAAATAGTATGAAAATTGTGCTATTAGCTAAGAAAAAACTAATCCGCAATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGG  >  NC_000913/987268‑987567

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: