Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 394,153 C→A 10.0% pseudogene (1184/1449 nt) yaiT → pseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913394,1530CA10.0% 312.2 / 14.3 124pseudogene (1184/1449 nt)yaiTpseudogene, autotransporter family;putative structure; Not classified; interrupted by IS3; putative flagellin structural protein
Reads supporting (aligned to +/- strand):  ref base C (53/57);  new base A (11/3);  total (64/60)
Fisher's exact test for biased strand distribution p-value = 4.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

CGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGG  >  NC_000913/394004‑394302
                                                                                                                                                     |                                                                                                                                                     
cGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGc                                                                                                                                                       >  2:1231842/1‑150 (MQ=255)
cGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGc                                                                                                                                                       <  2:652311/150‑1 (MQ=255)
 gACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCa                                                                                                                                                      >  2:1623840/1‑150 (MQ=255)
    ttCAAACTGGGGGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGt                                                                                                                                                   <  2:835311/150‑1 (MQ=255)
      aaaaCTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGAGTTTCTGTTCTGTTCCAGTAAGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTCTGAACTTGGTGTAGATGCAGGtgt                                                                                                                                                 >  2:1142928/2‑148 (MQ=255)
        aaCTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTtata                                                                                                                                               <  2:711877/150‑1 (MQ=255)
        aaCTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTtata                                                                                                                                               <  2:718624/150‑1 (MQ=255)
                gATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAATGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAcc                                                                                                                                       <  1:305302/150‑1 (MQ=255)
                 aTGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCt                                                                                                                                      <  2:432840/150‑1 (MQ=255)
                  tGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTa                                                                                                                                     >  1:584041/1‑150 (MQ=255)
                  tGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTa                                                                                                                                     >  2:1156277/1‑150 (MQ=255)
                   gCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGATCAGTCTTACGTCAGCATGCGTTATGAACTGGGTGTAGATGAAGGTTATACCTTCACCGAc                                                                                                                                    >  2:964902/1‑150 (MQ=255)
                      ggTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGc                                                                                                                                 >  2:383992/1‑150 (MQ=255)
                      ggTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGc                                                                                                                                 >  2:1775639/1‑150 (MQ=255)
                       gtgaCGGGACTCCTTACCGCAGCGTTTCTGGGCTCGTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCg                                                                                                                                <  2:1232018/147‑1 (MQ=255)
                       gTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCg                                                                                                                                >  1:351043/1‑150 (MQ=255)
                          aCGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCATTCTTGTGATGACTACCAGCTGAGAAACGAAATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAg                                                                                                                             <  2:906163/150‑1 (MQ=255)
                           cGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATCAAGGTTATACCATCAACTAAAgataaaa                                                                                                                            >  2:1362849/1‑144 (MQ=255)
                            gTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGaa                                                                                                                           >  2:999573/1‑149 (MQ=255)
                             tGACTCCTTGCGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATc                                                                                                                          >  1:1799601/1‑150 (MQ=255)
                             tGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGaac                                                                                                                          >  2:1686155/1‑148 (MQ=255)
                              gACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCa                                                                                                                         <  1:145150/150‑1 (MQ=255)
                               aCTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAg                                                                                                                        <  2:533948/150‑1 (MQ=255)
                               aCTCCTTACGGAAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAg                                                                                                                        <  2:1379141/150‑1 (MQ=255)
                                  ccTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAACTACAGCGAAGATCAGGct                                                                                                                     >  1:1405444/1‑150 (MQ=255)
                                  ccTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGAAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGct                                                                                                                     >  1:498020/1‑150 (MQ=255)
                                      aCGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGa                                                                                                                 >  1:447395/1‑150 (MQ=255)
                                         gCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGGCTc                                                                                                              >  1:686530/1‑150 (MQ=255)
                                         gCAGCGTTTCTGGTCTGTTCCAGTCTGGAGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTAACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACACCGAAGATCAGGCTCTGACTc                                                                                                              >  1:216717/1‑150 (MQ=255)
                                               tttCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACt                                                                                                        <  2:593046/150‑1 (MQ=255)
                                               tttCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACt                                                                                                        <  1:1835645/150‑1 (MQ=255)
                                                 tCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGTGTGTAGATGCAGGTTATACCATAACCTACAACAAAGATCAGGCTCTGACTCCGTACtta                                                                                                      >  2:376957/1‑149 (MQ=255)
                                                  cTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCa                                                                                                     >  1:748845/1‑150 (MQ=255)
                                                   ttgtatGTTCCAGTCTGGTGATGACTACCAGCTTATAAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGTACTGGGTGTAGATGAAGGTAATACATTCACCTACAACCAAGAAATATCTCTGAATCCTTACTTCaa                                                                                                    >  2:1000603/6‑150 (MQ=255)
                                                   tggtctggtcCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCaa                                                                                                    <  2:1598882/150‑1 (MQ=255)
                                                   tGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATCCCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTAATTCaa                                                                                                    >  1:272998/1‑150 (MQ=255)
                                                   tGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCaa                                                                                                    <  2:228041/150‑1 (MQ=255)
                                                   tGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGAAAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCaa                                                                                                    >  2:474168/1‑150 (MQ=255)
                                                         gTTCCAGGCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGc                                                                                              <  1:1560003/150‑1 (MQ=255)
                                                          ttCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATGCCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCt                                                                                             <  1:664418/150‑1 (MQ=255)
                                                           tCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCtt                                                                                            >  1:650614/1‑150 (MQ=255)
                                                           tCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCtt                                                                                            >  1:1841985/1‑150 (MQ=255)
                                                            ccAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTa                                                                                           >  1:781258/1‑150 (MQ=255)
                                                            ccAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGGCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTa                                                                                           <  1:192566/150‑1 (MQ=255)
                                                             cAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTAc                                                                                          >  1:1664030/1‑150 (MQ=255)
                                                              aGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACg                                                                                         <  2:1766719/150‑1 (MQ=255)
                                                               gTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACACCAGGCGATATGAACTGGGTGTAGATGCAGGTTATACCTTCACCGACAGCGAATATCAGGCTCTGACTCCTTACTTCAAACTGGCTTACGt                                                                                        >  1:1761484/1‑150 (MQ=255)
                                                                  tGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTa                                                                                     <  2:618544/150‑1 (MQ=255)
                                                                   ggTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTAAACTAAAGATAATATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTac                                                                                    >  2:218789/1‑150 (MQ=255)
                                                                    gTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGCCTCCGTACTTCAAACTGGCTTACGTCTacg                                                                                   >  1:913887/1‑150 (MQ=255)
                                                                    gTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTacg                                                                                   <  2:1552508/150‑1 (MQ=255)
                                                                     tgatgaCTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTacga                                                                                  <  2:748845/150‑1 (MQ=255)
                                                                       atgaCTACCAGCTGAGCAACGACATGAAAGATGAAGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTacgacg                                                                                <  2:1718828/150‑1 (MQ=255)
                                                                        tgaCTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTAcgacga                                                                               >  1:1330883/1‑150 (MQ=255)
                                                                          aCTACCATCTGATCAACGACATGAAAGTTGACGGTCATTCTTTAGACAGCATGCGGTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGAct                                                                             <  2:1685586/150‑1 (MQ=255)
                                                                           cTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATCAAGGTTATACCATAAACAAAAAAAAAAATAAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGActc                                                                            >  2:875678/1‑150 (MQ=255)
                                                                              ccAGCTGAGCAACGACATGTAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATCAAGGTTATACCATAAACTAAAAAAAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTaa                                                                         >  2:641030/1‑150 (MQ=255)
                                                                              ccAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTaa                                                                         >  1:464950/1‑150 (MQ=255)
                                                                              ccAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCATCACCTACAGCGAAGATAAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTaa                                                                         >  2:153500/1‑150 (MQ=255)
                                                                                    gAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaccga                                                                   >  1:389322/1‑150 (MQ=255)
                                                                                    gAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacga                                                                   <  1:802654/150‑1 (MQ=255)
                                                                                      tcAACGACATGAAAGTTTACGGGAAGTCTTACTACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgata                                                                 <  1:1205562/149‑1 (MQ=255)
                                                                                      ccaccGACATGAAAGTTGACGGTCAGTCTAACGACAGCATGCGTTATGAACTGGGTGTTGATGCAGGTTATACCTTCACCTACAGCGAAGATCCGGCTCTGACTCCGTACTTCAAACTGGCTTACGGCTACGACGACTCTAACaacgata                                                                 <  1:1576155/146‑1 (MQ=255)
                                                                                          cGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgataacga                                                             <  1:1156277/150‑1 (MQ=255)
                                                                                          cGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGAAGGTTATACCATCCCCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTACCaacgataacgc                                                             >  2:1691283/1‑149 (MQ=255)
                                                                                              aTGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGAtgtg                                                         <  2:963884/150‑1 (MQ=255)
                                                                                               tGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATCAAGGTTATCCCATAAACAAAAAAAAAAAAAAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGa                                                        >  2:1307646/1‑150 (MQ=255)
                                                                                                gAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGaa                                                       >  2:1804838/1‑150 (MQ=255)
                                                                                                 aaaGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAAc                                                      >  1:184967/1‑150 (MQ=255)
                                                                                                       gACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGAt                                                >  2:1114749/1‑150 (MQ=255)
                                                                                                        ccGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGAtt                                               <  2:1534198/149‑1 (MQ=255)
                                                                                                        ccGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGAtt                                               <  1:721325/149‑1 (MQ=255)
                                                                                                             cAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATc                                          <  1:1623840/150‑1 (MQ=255)
                                                                                                              aGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGTTGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCTTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCg                                         <  2:1663996/150‑1 (MQ=255)
                                                                                                              aGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCTTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCg                                         <  2:1109959/150‑1 (MQ=255)
                                                                                                                 cTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGATATACATTCACATACAGCGAAGAACAGGCTCTTACTCCGAACTTCAAACTAGCTTACGTCTACGAAGACTCTAACAACGATAACCAAGTGAACGCCGATACCATCCATa                                      >  1:227098/1‑150 (MQ=255)
                                                                                                                     cGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACgg                                  >  1:1601744/1‑150 (MQ=255)
                                                                                                                     cGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACgg                                  >  2:750034/1‑150 (MQ=255)
                                                                                                                       aCAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTa                                <  2:681624/150‑1 (MQ=255)
                                                                                                                          gCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTg                             <  1:1752987/150‑1 (MQ=255)
                                                                                                                           cATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGa                            >  1:1801673/1‑150 (MQ=255)
                                                                                                                             tGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAg                          <  2:747534/150‑1 (MQ=255)
                                                                                                                             tGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGAGCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAg                          <  2:1369701/150‑1 (MQ=255)
                                                                                                                              gCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAgg                         >  1:1295273/1‑150 (MQ=255)
                                                                                                                               cGTTATTAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGGCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAggg                        <  1:1641496/150‑1 (MQ=255)
                                                                                                                               cGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCTAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAggg                        >  2:1215705/1‑150 (MQ=255)
                                                                                                                               cGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTCGCTTACGTCTACGCCGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATACCGCTACTGAAggg                        >  1:1348460/1‑150 (MQ=255)
                                                                                                                                    tGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCTAAGATCAGGCTCTGACTCCGTACTTACAACTGGCTTACGTCTACGACGACTCTACCACCGCTAACGATGTGAACGGCGCTTCCATCGATAACGGTACTGAAGGGTCTGc                   >  1:101753/1‑150 (MQ=255)
                                                                                                                                     gAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCg                  >  2:1601830/1‑150 (MQ=255)
                                                                                                                                      aaCTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCgg                 <  2:210788/150‑1 (MQ=255)
                                                                                                                                        cTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTa               >  1:731757/1‑150 (MQ=255)
                                                                                                                                           gctgtAGGTGAAGGTATAACCTGCATCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACgt            >  1:1441090/3‑150 (MQ=255)
                                                                                                                                            ctgtAGGTGAAGGTATAACCTGCATCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACgtg           >  2:8822/2‑150 (MQ=255)
                                                                                                                                             tttttaTGCTGGTAATACTTTCACATACAGCTAAGATCTGGATCTTACTCGTTATTTAAAACTGGCTTACTTCTACGACGACTATAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACgtgt          <  2:1841985/145‑1 (MQ=255)
                                                                                                                                             tgtAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACgtgt          >  1:3370/1‑150 (MQ=255)
                                                                                                                                               tAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTg        >  1:1119990/1‑150 (MQ=255)
                                                                                                                                                 gATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGCCGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTggg      >  1:31194/1‑149 (MQ=255)
                                                                                                                                                 gATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTggg      >  1:193454/1‑149 (MQ=255)
                                                                                                                                                  aTGCAGGTTATACCTTCACATACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTc     <  1:1454563/150‑1 (MQ=255)
                                                                                                                                                   tGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCt    >  1:743761/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGG  >  NC_000913/394004‑394302

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