Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 394,153 C→A 11.6% 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,1530CA11.6% 230.4 / 11.2 92pseudogene (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 (32/48);  new base A (9/3);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 3.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTG  >  NC_000913/394006‑394301
                                                                                                                                                   |                                                                                                                                                    
aCTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGAGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCAGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAg                                                                                                                                                    >  2:1240260/1‑150 (MQ=255)
 cTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAgg                                                                                                                                                   >  1:420415/1‑150 (MQ=255)
 cTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAgg                                                                                                                                                   >  1:1197556/1‑150 (MQ=255)
      aaCTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTtata                                                                                                                                              >  1:399111/1‑150 (MQ=255)
       aCTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATAc                                                                                                                                             >  2:1122422/1‑150 (MQ=255)
        cTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGCCATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATAcc                                                                                                                                            >  2:248931/1‑150 (MQ=255)
        cTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATAcc                                                                                                                                            <  2:1235104/150‑1 (MQ=255)
        cTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATAcc                                                                                                                                            >  1:1107641/1‑150 (MQ=255)
            gTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCa                                                                                                                                        <  2:1116871/150‑1 (MQ=255)
             tGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAc                                                                                                                                       <  2:276625/150‑1 (MQ=255)
              gATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAcc                                                                                                                                      >  2:1462337/1‑150 (MQ=255)
              gATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAcc                                                                                                                                      <  2:1143305/150‑1 (MQ=255)
               aTGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCt                                                                                                                                     <  2:667765/150‑1 (MQ=255)
                tGCTGGATACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGATGACGCTCAGTCTTACGACAGCATGCGTTATGAACTGGCTGTAGAAGCAGGATATACCTTCACCTa                                                                                                                                    >  2:874332/1‑150 (MQ=255)
                 gCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACTGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTAc                                                                                                                                   <  2:887813/150‑1 (MQ=255)
                    ggTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGc                                                                                                                                >  2:920814/1‑150 (MQ=255)
                     gTTACGTGACTCCTTATGGCTGGGGTTCTGGTCTGTTCCAGTCTTGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGTGTGTAGGTGAAGGTATAACCTGCATCTACaccc                                                                                                                               <  2:868632/150‑4 (MQ=255)
                       tACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGaa                                                                                                                             <  2:1087023/150‑1 (MQ=255)
                        aCGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAACTACAGCGAAg                                                                                                                            >  2:1209022/1‑150 (MQ=255)
                             aCTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAg                                                                                                                       <  2:1521131/150‑1 (MQ=255)
                             aCTCCTTACGGAAGCGTTGCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACTTTAAAGTTTACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAg                                                                                                                       <  2:1493433/150‑1 (MQ=255)
                              cTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAgg                                                                                                                      <  1:435258/150‑1 (MQ=255)
                               tCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGc                                                                                                                     <  1:487912/150‑1 (MQ=255)
                               tCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGc                                                                                                                     <  1:998407/150‑1 (MQ=255)
                                ccTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGCTCAGGct                                                                                                                    >  1:1149345/1‑150 (MQ=255)
                                    aCGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCAACTACAGCGAAGATCAGGCTCTGa                                                                                                                >  2:545767/1‑150 (MQ=255)
                                     cGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTACGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGAc                                                                                                               <  1:1519573/150‑1 (MQ=255)
                                        cAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTcc                                                                                                            <  1:859124/150‑1 (MQ=255)
                                            gTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTAc                                                                                                        >  1:1411919/1‑150 (MQ=255)
                                            gTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTAc                                                                                                        <  2:964889/150‑1 (MQ=255)
                                             tttCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATCCCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTCCt                                                                                                       >  2:327560/1‑150 (MQ=255)
                                             ggtCTGGTCTGTTCCAGTCTGGGGATGACTACCCGCAGAGCCACGACATCACAGTTGACGCTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCCCCTACAGCGAAGATCAGGCTCTGACTCCGTACt                                                                                                       <  2:329139/148‑1 (MQ=255)
                                              ttCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACtt                                                                                                      <  1:700884/150‑1 (MQ=255)
                                                 tGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCaa                                                                                                   >  1:1554650/1‑150 (MQ=255)
                                                  ggTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGATATGAACTGGGTGTAGATCAAGGTTATACCTTAAACTACAACAAAAAAAAGGCTCTGACTCCGTCCTTCaca                                                                                                  >  2:1159333/1‑148 (MQ=255)
                                                    tCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACt                                                                                                <  2:422312/150‑1 (MQ=255)
                                                        ttCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCt                                                                                            <  1:545767/150‑1 (MQ=255)
                                                        ttCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCt                                                                                            <  1:676447/150‑1 (MQ=255)
                                                          ccAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTa                                                                                          <  2:506847/150‑1 (MQ=255)
                                                             gTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGt                                                                                       <  1:1180763/150‑1 (MQ=255)
                                                                 ggTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGCCTCCGTACTTCAAACTGGCTTACGTCTac                                                                                   >  1:1022034/1‑150 (MQ=255)
                                                                 ggTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTac                                                                                   >  1:641764/1‑150 (MQ=255)
                                                                  gTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTacg                                                                                  <  1:1637829/150‑1 (MQ=255)
                                                                  gTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTACAAACTGGCTTACGTCTacg                                                                                  >  2:1066750/1‑150 (MQ=255)
                                                                     atgaCTACCAGCTGAGCAACGACATGAATGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTacgacg                                                                               <  2:433383/150‑1 (MQ=255)
                                                                        aCTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATCAAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGAct                                                                            >  2:572779/1‑150 (MQ=255)
                                                                         cTACCAGCTGAGCAACGACATGAAAGTTGCCGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGActc                                                                           >  1:82882/1‑150 (MQ=255)
                                                                         cTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGActc                                                                           <  1:853966/150‑1 (MQ=255)
                                                                         cTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGActc                                                                           <  2:1480861/150‑1 (MQ=255)
                                                                          aaCCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGActct                                                                          <  2:620651/149‑1 (MQ=255)
                                                                             cAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGCCTTACGTCTACGACGACTCTaac                                                                       >  1:659368/1‑150 (MQ=255)
                                                                              aGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATCAAGGTTATACCATAAACAAAAGAAAAAATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTaaca                                                                      >  2:818701/1‑150 (MQ=255)
                                                                                   cgCAACGACATGAAAGTTGACGGTCAGTCTTAGGAAAGCATGCGTTATGAACTGGGTGTAGATGCAGGTAATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgat                                                                 <  1:321489/149‑1 (MQ=255)
                                                                                   aGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTAAAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgat                                                                 >  2:1416458/1‑150 (MQ=255)
                                                                                    gCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGAAGGTTATACCATAAACTAAAAAAAAAAAAAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgatt                                                                >  2:1238841/1‑149 (MQ=255)
                                                                                     cAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgataa                                                               >  1:1076666/1‑150 (MQ=255)
                                                                                      aaCGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgataac                                                              <  1:166054/150‑1 (MQ=255)
                                                                                       aCGACATTAATGTTGACGGTCAGTCTTACGACAGCTTGGGTTATGAACTGGGTGTAGATGCAGTTTTTACCTTCACCTACTGCGAAGATCTGGCTCTGACTCCTTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgataacg                                                             <  2:1410682/150‑1 (MQ=255)
                                                                                       aCGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgataacg                                                             <  1:920814/150‑1 (MQ=255)
                                                                                        cGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACaacgataacga                                                            <  1:1122422/150‑1 (MQ=255)
                                                                                          aCATGTAATTTGACGGTCAGTCTTACGACTGCATGCGTTATGATCTGGTTGTAGATGCAGGTTATACATTCAGCTACAGCGAATATCAGGCTCTGACTCCGTACTTCAAACTGGATTACGTCTAGGACGACTCTAACAACGATAACGAtg                                                          <  2:1692585/150‑1 (MQ=255)
                                                                                           cATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCAACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGAtgt                                                         >  2:1399531/1‑150 (MQ=255)
                                                                                            aTGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGAtgtg                                                        <  2:562070/150‑1 (MQ=255)
                                                                                            aTGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGAAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACCAtgtg                                                        >  1:121670/1‑150 (MQ=255)
                                                                                                  gTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCTTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGc                                                  <  2:1650492/150‑1 (MQ=255)
                                                                                                  gTTGACGGTCAGTCTTACGACAGCATGCGTGATGAACTGGGTGTAGATGCAGGTTATAACTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGc                                                  >  1:689210/1‑150 (MQ=255)
                                                                                                     gACGGTCATTCTTACGACAGCATCCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGAt                                               <  2:1183247/150‑1 (MQ=255)
                                                                                                        ggTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTcc                                            >  1:394553/1‑150 (MQ=255)
                                                                                                        ggTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTcc                                            >  1:1158402/1‑150 (MQ=255)
                                                                                                           cAGTCTTAGTCCAGCATGCGTTATGAACTGGGTGTAGATGCCGGTTATACCTTCCCCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATc                                         <  2:62948/150‑1 (MQ=255)
                                                                                                           cAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATc                                         <  2:926022/150‑1 (MQ=255)
                                                                                                               tttACTACAGCATGCGTTTTTAACTGGGTGTAGATGCAGTTTATACCTTCACCTACAGCTAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATa                                     <  2:1325697/149‑1 (MQ=255)
                                                                                                               cTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATa                                     >  1:358420/1‑150 (MQ=255)
                                                                                                                ttACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATaa                                    <  1:473673/150‑1 (MQ=255)
                                                                                                                 tACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAAc                                   <  1:961041/150‑1 (MQ=255)
                                                                                                                 tACGACAGCATGCGTTATGAACTGGGTGTAGATCAAGGTTATACCATAAACTAAATCAAAAATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTCCGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGCTacc                                   >  2:1098228/1‑148 (MQ=255)
                                                                                                                   cGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACgg                                 <  1:1240260/150‑1 (MQ=255)
                                                                                                                              gTTATGAACTGGGTGTAGATCAAGGTTATACCATAAACAAAAAAAAAAAAAAGGCTCTGACTCCGAACTTCAAACTGGCTTACTTCTACGACGACTCTAACAACGAAAACGATGTGAACTGCTATTCACTCGATAACGGTAATGAAGGGt                      >  2:1204287/1‑150 (MQ=255)
                                                                                                                               ttATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTc                     <  2:1197556/150‑1 (MQ=255)
                                                                                                                                  tGAACTGGGTGTAGATGCAGGTTATACCTTCACATACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGc                  <  2:1654091/150‑1 (MQ=255)
                                                                                                                                   gtaCTGGGTGTAGATGCAGGGTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCg                 <  2:1107641/148‑1 (MQ=255)
                                                                                                                                   gAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCg                 <  2:420415/150‑1 (MQ=255)
                                                                                                                                     cttcgctgtAGGTGAAGGTATAACCTGCATCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGt               >  2:904305/7‑150 (MQ=255)
                                                                                                                                         ggTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACgt           >  2:1647611/1‑150 (MQ=255)
                                                                                                                                         gctgtAGGTGAAGGTATAACCTGCATCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGAGGTGACCGGCGATTCCATCGATAACGGTACTGCAGGGTCTGCTGTACgt           >  2:1562873/3‑150 (MQ=255)
                                                                                                                                          ggttAGATGATGGTTATACCTTCACCTACAGCGAAGAGCTGGCTCTGACTCCGTACTTCAATCTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACgtg          <  1:695636/147‑1 (MQ=255)
                                                                                                                                            gtAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGtt        <  2:1411919/150‑1 (MQ=255)
                                                                                                                                              aGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGCCGACTCTAACAACGATAACGATGTGAACGGCGATTCCCTCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTgg      >  1:1516079/1‑150 (MQ=255)
                                                                                                                                              aGATGAAGGTTATACCTTCACCTAAAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTATACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTgg      <  2:1490760/150‑1 (MQ=255)
                                                                                                                                               gATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGt     >  1:310555/1‑150 (MQ=255)
                                                                                                                                                  gCAGGTTATACCTTCACCTACAGCTAAGATCAGGCGCTGACTCCGTACTTCAACCTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTg  <  1:592805/150‑1 (MQ=255)
                                                                                                                                                  gCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTg  >  2:163173/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ACTTCAAACTGGGTGATGCTGGTTACGTGACTCCTTACGGCAGCGTTTCTGGTCTGTTCCAGTCTGGTGATGACTACCAGCTGAGCAACGACATGAAAGTTGACGGTCAGTCTTACGACAGCATGCGTTATGAACTGGGTGTAGATGCAGGTTATACCTTCACCTACAGCGAAGATCAGGCTCTGACTCCGTACTTCAAACTGGCTTACGTCTACGACGACTCTAACAACGATAACGATGTGAACGGCGATTCCATCGATAACGGTACTGAAGGGTCTGCGGTACGTGTTGGTCTG  >  NC_000913/394006‑394301

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: