Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 394,153 C→A 8.4% 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,1530CA 8.4% 257.2 / 4.8 97pseudogene (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 (37/50);  new base A (8/2);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 4.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTGG  >  NC_000913/394004‑394302
                                                                                                                                                     |                                                                                                                                                     
cGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGc                                                                                                                                                       >  1:593365/1‑150 (MQ=255)
cGACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGc                                                                                                                                                       <  1:562301/150‑1 (MQ=255)
   cTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAgg                                                                                                                                                    <  1:1160392/150‑1 (MQ=255)
   cTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAgg                                                                                                                                                    >  1:81822/1‑150 (MQ=255)
    ttCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGt                                                                                                                                                   >  2:1157358/1‑150 (MQ=255)
      cAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGGCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTta                                                                                                                                                 <  2:459545/150‑1 (MQ=255)
        aaCTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAgg                                                                                                                                                    >  2:648149/1‑145 (MQ=255)
        aaCTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAgg                                                                                                                                                    <  1:648149/145‑1 (MQ=255)
        aaCTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTtata                                                                                                                                               >  1:1279990/1‑150 (MQ=255)
        aaCTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTtata                                                                                                                                               <  2:1111956/150‑1 (MQ=255)
              gTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCa                                                                                                                                         >  2:293975/1‑150 (MQ=255)
              gTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCa                                                                                                                                         <  1:293975/150‑1 (MQ=255)
              gTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCa                                                                                                                                         >  2:879534/1‑150 (MQ=255)
               tGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAc                                                                                                                                        <  2:32289/150‑1 (MQ=255)
               tGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAc                                                                                                                                        <  2:133471/150‑1 (MQ=255)
                taTGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTTGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACTACATGAAAGTTGACTGTCAGTCTTACGACAGCATGCTTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAcc                                                                                                                                       <  1:464207/149‑1 (MQ=255)
                gATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAcc                                                                                                                                       >  2:27092/1‑150 (MQ=255)
                     tGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAg                                                                                                                                  >  2:31423/1‑150 (MQ=255)
                      ggTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGc                                                                                                                                 <  1:1575599/150‑1 (MQ=255)
                       gTTACGTGACTCCTTACGGCAGCGTTTCTGGGCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGGTGAAGGTATAACCTGCATCTACaccc                                                                                                                                <  1:327629/150‑4 (MQ=255)
                         tACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGaa                                                                                                                              >  1:1686049/1‑150 (MQ=255)
                         tACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGaa                                                                                                                              <  2:955761/150‑1 (MQ=255)
                          aCGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAg                                                                                                                             <  1:1460199/150‑1 (MQ=255)
                            gTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGAt                                                                                                                           <  2:749026/150‑1 (MQ=255)
                             tGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATc                                                                                                                          >  2:1643864/1‑150 (MQ=255)
                              gACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCa                                                                                                                         <  1:206623/150‑1 (MQ=255)
                               aCTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAg                                                                                                                        <  2:1000300/150‑1 (MQ=255)
                               aCTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAg                                                                                                                        <  1:700881/150‑1 (MQ=255)
                                cTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAgg                                                                                                                       >  1:1298422/1‑150 (MQ=255)
                                 tCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGc                                                                                                                      <  2:563760/150‑1 (MQ=255)
                                   cTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGctc                                                                                                                    <  2:570978/150‑1 (MQ=255)
                                    ttACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGct                                                                                                                     <  2:1687135/148‑1 (MQ=255)
                                    ttACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGct                                                                                                                     >  1:1687135/1‑148 (MQ=255)
                                       cGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGAc                                                                                                                >  1:691441/1‑150 (MQ=255)
                                        ggCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACt                                                                                                               >  1:681970/1‑150 (MQ=255)
                                              gTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTAc                                                                                                         <  2:1251388/150‑1 (MQ=255)
                                                ttatGGTGTGTTCCAGTCTGGTGATGACTACCATATGAGAAAAGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACtt                                                                                                       <  2:84586/147‑1 (MQ=255)
                                                ttCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACtt                                                                                                       <  2:452575/150‑1 (MQ=255)
                                                   tGGTCTGTTCCAGTCTGGTGATGACTACCAGCGGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCaa                                                                                                    <  2:1171107/150‑1 (MQ=255)
                                                    ggTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCCCCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCaaa                                                                                                   >  1:299436/1‑150 (MQ=255)
                                                     gTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAAc                                                                                                  <  2:700515/150‑1 (MQ=255)
                                                     gTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAAc                                                                                                  <  2:906419/150‑1 (MQ=255)
                                                      tCTGTTCCAGTCTGGTGATGACTTCCATGTGAGCAACGTCATGTAAGTTTACGGGCAGTCTTACGAAAGCATTCGGTATTAACTGGGTGTAGATTCAGGGTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACt                                                                                                 <  2:1265204/150‑1 (MQ=255)
                                                         gTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGc                                                                                              <  2:1260672/150‑1 (MQ=255)
                                                           tCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCtt                                                                                            <  1:937867/150‑1 (MQ=255)
                                                           tCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCtt                                                                                            >  1:1282713/1‑150 (MQ=255)
                                                           tCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTAAAGCGAAGATCAGGCACTGACTCCGTACGGCAAACGGGCtt                                                                                            >  2:291644/1‑150 (MQ=255)
                                                                  tGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTATATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTa                                                                                     >  2:761262/1‑150 (MQ=255)
                                                                   ggTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTac                                                                                    <  1:361593/150‑1 (MQ=255)
                                                                    gTGATGACTACCAGCTGAGCAACGACATGAAAGTTGAAGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATCACGGTTATCCCATAAACAAAAACAAAAATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTacg                                                                                   >  2:1554598/1‑150 (MQ=255)
                                                                       atgaCTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTacgacg                                                                                <  2:1474589/150‑1 (MQ=255)
                                                                        tgaCTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTAcgacga                                                                               >  1:1586324/1‑150 (MQ=255)
                                                                        tgaCTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTAcgacga                                                                               >  1:593681/1‑150 (MQ=255)
                                                                        tgaCTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCCAAGATCAGGCTCTGAGTCAGTACTTCAATCTGGCTTACGTCTAcgacga                                                                               >  2:457264/1‑150 (MQ=255)
                                                                            tACCAGCTGTGAAAGGTCATGAATGTTGACGGTCAGTCTTACGAAAGCATGCGTTATGAAGTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTTACTCCGTACTTCAAAGTGGCTTACGTATACGACGActct                                                                           <  2:1540585/150‑1 (MQ=255)
                                                                            caCCAGCTGTGGAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCGGTACTTCAAACTGGCTTACGTCTACGACGActct                                                                           <  1:266340/149‑1 (MQ=255)
                                                                             aCCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTa                                                                          >  1:713433/1‑150 (MQ=255)
                                                                               ctgCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAAGTGGGTGTATATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTAGGTCTACGACGACTCTaac                                                                        <  2:527154/148‑1 (MQ=255)
                                                                                 gCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTaacaa                                                                      >  1:50480/1‑150 (MQ=255)
                                                                                     aGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgat                                                                  <  2:65452/150‑1 (MQ=255)
                                                                                      gCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCTAACTGGCTTACGTCTACGACGACTCTAACaacgata                                                                 >  1:1095373/1‑150 (MQ=255)
                                                                                        aaCGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGAAGGTTATACCTTAAACTAAATCGAATATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTTCGACGACTCTAACaacgataac                                                               >  2:1091802/1‑150 (MQ=255)
                                                                                          cGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATCAAGGTTATACCATAAACAAAAAAAAAAAACAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgataacga                                                             >  2:978940/1‑150 (MQ=255)
                                                                                           gACATGAATGTTGAGGGGCAGTCTTACGACAGCTTGCGTTATGAAGTGGGTGTAGATGCCGGTTATACCTTCACCTACAGCGAAGGTCTGGCTCTGACTCCGTACTTCAAAGTGGCTTACGTCTACGACGACTCTAACaacgataacgat                                                            <  2:1438007/150‑1 (MQ=255)
                                                                                           gACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgataacgat                                                            >  2:1384560/1‑150 (MQ=255)
                                                                                             cATTAAAGTTGAGGGTCAGTATTACGAAAGCATGCGTTTTGAACTGGGTTTAGATGAAGTTTATACCTTCACCTACAGCGAATTTCAGTCTCTGACTACGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGAtgt                                                          <  1:1440350/150‑1 (MQ=255)
                                                                                             cATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGAAGGTTATACCATAACCTAAACAGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGAtgt                                                          >  2:508325/1‑150 (MQ=255)
                                                                                               tGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGa                                                        >  1:996941/1‑150 (MQ=255)
                                                                                               tGAAAGTTGACGGTCAGTATTACGACAGCATGCGTTTTGAACTGGGTGTAGATGCTGGTTTTACCTTCACCTACATCGAAGTTCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGa                                                        <  2:1171772/150‑1 (MQ=255)
                                                                                               tGAAAGTTGACGGGCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGTTTACGTCTACGACGACTCTAACAACGATAACGATGTGa                                                        <  1:879534/150‑1 (MQ=255)
                                                                                                 aaaGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAAc                                                      >  1:1388042/1‑150 (MQ=255)
                                                                                                    gTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGc                                                   >  1:774373/1‑150 (MQ=255)
                                                                                                         cGGTCATGCTTACGCAAGCATGCGTTCTTCACTGGGTGTAGATGCAGGTTATACCTTCCCCTCCAGCTAAGATCAGGCGCGGACTCCGTACTTAAAACTGGCGTACGGCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTc                                              <  2:882537/150‑1 (MQ=255)
                                                                                                         cGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGAtta                                              >  1:175040/1‑149 (MQ=255)
                                                                                                             ctgTCTTACTACAGCTTGCGTTATGACCTGGGTGTAGATGCAGGTTATACCTTCACCTACCGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATc                                          <  1:1511447/148‑1 (MQ=255)
                                                                                                              attCTTACGACAGCATTCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCg                                         <  2:1153822/148‑1 (MQ=255)
                                                                                                                 attACTACTGCATGCTTTATGAACGGGGTGTATATGCTGGTTATACCTTAACCTACTGCGAATATCAGGCTCTGTCTCCGTACTTCAAACTGGCTTACGTTTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATa                                      <  2:431595/149‑1 (MQ=255)
                                                                                                                      gACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGt                                 >  1:268228/1‑150 (MQ=255)
                                                                                                                       aCAGCATGCGTTATGAACTGGGTGTAGATCAAGGTTATACCTTAACCTAAAAAGAATATCAGGCTCTGACGCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGAGAACGGTa                                >  2:169035/1‑150 (MQ=255)
                                                                                                                          gCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTg                             >  2:385505/1‑150 (MQ=255)
                                                                                                                          gCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTg                             <  1:31423/150‑1 (MQ=255)
                                                                                                                           cATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGCCTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGa                            >  1:1184231/1‑150 (MQ=255)
                                                                                                                               ggTTATGAATTGGGTGTATATGCAGTTTATACCTTCACCTACAGCGAATATCAGGCTCTGACTCCGTACTTCAAATTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAggg                        <  2:175040/149‑1 (MQ=255)
                                                                                                                               cGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAggg                        >  1:1440364/1‑150 (MQ=255)
                                                                                                                               cGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAggg                        <  2:774373/150‑1 (MQ=255)
                                                                                                                                   aTGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTg                    >  1:1611611/1‑150 (MQ=255)
                                                                                                                                    tGAACTGGGTGTAGATCAAGGTTATACCTTAAACTAAAAATAAAAACAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGc                   >  2:55512/1‑150 (MQ=255)
                                                                                                                                     taaCTGGTTGTAGATGCAGGTAATACCTTCACCTACAGCGAAGATCAGGCTAGTACTCCGTACTTCAAACTGGCTTCCGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCg                  <  1:1643864/149‑1 (MQ=255)
                                                                                                                                     gAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCg                  >  1:928840/1‑150 (MQ=255)
                                                                                                                                      aaCTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCgg                 <  1:893106/150‑1 (MQ=255)
                                                                                                                                       aCTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGt                <  2:81822/150‑1 (MQ=255)
                                                                                                                                        ttcgctgtAGGTGAAGGTATAACCTGCATCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTa               >  2:434056/6‑150 (MQ=255)
                                                                                                                                        cTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTa               <  2:1586324/150‑1 (MQ=255)
                                                                                                                                          gggTGTAGATCAAGGTTATACCATAAACAAAAAAAAAAATAAGGCTCTGACTCAGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATGCCATCGATACCGGTACTGAAGGGTCTGCGGTACg             >  2:250296/1‑150 (MQ=255)
                                                                                                                                              gtAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGtt         <  2:1686049/150‑1 (MQ=255)
                                                                                                                                               tAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTg        <  1:559017/150‑1 (MQ=255)
                                                                                                                                                 gATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGt      >  1:948282/1‑150 (MQ=255)
                                                                                                                                                  aTGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTc     <  2:1703614/150‑1 (MQ=255)
                                                                                                                                                   tGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCt    >  1:1447158/1‑150 (MQ=255)
                                                                                                                                                     cAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTgg  <  1:1157358/150‑1 (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: