Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 563,450 A→G T41A (ACC→GCC)  sfmF → putative fimbrial protein SfmF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913563,4500AG97.6% 289.7 / ‑4.4 84T41A (ACC→GCC) sfmFputative fimbrial protein SfmF
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (43/39);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCGGGTTCAGTTG  >  NC_000913/563319‑563585
                                                                                                                                   |                                                                                                                                       
aGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAac                                                                                                                                 >  2:1069747/1‑140 (MQ=255)
aGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAac                                                                                                                                 >  1:282456/1‑140 (MQ=255)
 ggTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAaca                                                                                                                                >  2:93938/1‑140 (MQ=255)
 ggTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAaca                                                                                                                                >  2:698/1‑140 (MQ=255)
   tACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGGTGATTTCTCATGTGCCGTAAAcaca                                                                                                                              >  1:488317/1‑140 (MQ=255)
    aCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAg                                                                                                                             <  2:906747/140‑1 (MQ=255)
    aCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAg                                                                                                                             <  2:826861/140‑1 (MQ=255)
    aCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAg                                                                                                                             >  2:1061023/1‑140 (MQ=255)
    aCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAg                                                                                                                             >  2:105005/1‑140 (MQ=255)
     cGGAGAATGAGAAGAGTTCTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGt                                                                                                                                     <  1:993848/131‑1 (MQ=255)
     cGGAGAATGAGAAGAGCACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGGACGGTAACGTTGTTGATTTCTCCTGTGCCGt                                                                                                                                     >  2:993848/1‑131 (MQ=255)
       gagaATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCgg                                                                                                                          <  1:958076/140‑1 (MQ=255)
         gaATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGat                                                                                                                        <  1:93938/140‑1 (MQ=255)
         gaATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGat                                                                                                                        <  2:801965/140‑1 (MQ=255)
         gaATGAGAAGAGTACCCTTTAGTTGATTCTGCGGGCTCCTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGat                                                                                                                        <  1:614508/140‑1 (MQ=255)
          aaTGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGata                                                                                                                       >  1:42352/1‑140 (MQ=255)
            tGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAg                                                                                                                             >  1:17618/1‑132 (MQ=255)
            tGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAg                                                                                                                             <  2:17618/132‑1 (MQ=255)
                    gTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGAc                                                                                                             >  2:1148699/1‑140 (MQ=255)
                     tACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTCCTCCTGTACCGTAAGCACAGCAGATATTGATAAGACg                                                                                                            >  1:653941/1‑140 (MQ=255)
                      aCTCTTTAGCTGTTTCTTCGGGTTACTGTGGGGTTTCAGCGGATGGGCCGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACgg                                                                                                           <  2:410956/140‑1 (MQ=255)
                      aCTCTTTAGCTGTTTCTTCGGGCTACTGTGGAGTTTCAGGGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACgg                                                                                                           <  2:488317/140‑1 (MQ=255)
                      aCTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACgg                                                                                                           <  2:213545/140‑1 (MQ=255)
                            tAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGAtt                                                                                                     >  2:95325/1‑140 (MQ=255)
                             aGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGAttt                                                                                                    >  2:1157660/1‑140 (MQ=255)
                             aGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAAAACAGCGAATATTGGTAAGACGGTAGAttt                                                                                                    >  1:410956/1‑140 (MQ=255)
                               cTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAg                                                                                                        <  1:1228863/134‑1 (MQ=255)
                               cTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAg                                                                                                        >  2:1228863/1‑134 (MQ=255)
                               cTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAg                                                                                                  >  1:791392/1‑140 (MQ=255)
                                tGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAgg                                                                                                 <  1:698/140‑1 (MQ=255)
                                  tttCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCa                                                                                               >  2:1257744/1‑140 (MQ=255)
                                  tttCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCa                                                                                               >  1:367920/1‑140 (MQ=255)
                                   ttCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAg                                                                                              <  1:947958/140‑1 (MQ=255)
                                    tCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGa                                                                                             <  2:42352/140‑1 (MQ=255)
                                    tCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGa                                                                                             <  1:1069747/140‑1 (MQ=255)
                                    tCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGa                                                                                             >  1:354009/1‑140 (MQ=255)
                                    tCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGa                                                                                             <  1:105005/140‑1 (MQ=255)
                                         gggCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGcc                                                                                        <  2:354009/140‑1 (MQ=255)
                                              aCTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGa                                                                                   <  1:1366116/140‑1 (MQ=255)
                                                 gtgGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGacac                                                                                <  1:1275181/140‑1 (MQ=255)
                                                       ttCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAactac                                                                          >  2:288886/1‑140 (MQ=255)
                                                       ttCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAactac                                                                          >  2:1175072/1‑140 (MQ=255)
                                                       ttCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTCGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAactac                                                                          >  1:805778/1‑140 (MQ=255)
                                                         ccAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTg                                                                        >  1:1451518/1‑140 (MQ=255)
                                                         ccAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTg                                                                        >  1:296344/1‑140 (MQ=255)
                                                          cAGTGGATGGGCAGTGGACCCTATAAGAATGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGa                                                                       <  1:515281/140‑1 (MQ=255)
                                                             tGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACg                                                                    <  2:367920/140‑1 (MQ=255)
                                                             tGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACg                                                                    <  1:1058881/140‑1 (MQ=255)
                                                               gATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCt                                                                  >  1:912426/1‑140 (MQ=255)
                                                                 tGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGCTTTAGGCAAATGGCCTACGACCCAACTACTAAACGCTgg                                                                >  2:740648/1‑140 (MQ=255)
                                                                     cAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGAt                                                            <  2:296344/140‑1 (MQ=255)
                                                                     cAGTTGACCCTTTAGGAACGATTAATATCAATTCGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGAt                                                            <  2:805778/140‑1 (MQ=255)
                                                                      aGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATa                                                           >  2:1348748/1‑140 (MQ=255)
                                                                         tGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTATATTTAGGCGGATGGCCCGCGACACAACTACGGAACGCTGGCGATACCa                                                        >  2:754061/1‑140 (MQ=255)
                                                                              cTTTAGGACCGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCa                                                   <  1:288886/140‑1 (MQ=255)
                                                                              cTTTAGGAACGATTAGTATCAATTTGCACAGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCa                                                   <  1:754061/140‑1 (MQ=255)
                                                                                 tAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTc                                                <  2:1451518/140‑1 (MQ=255)
                                                                                  aGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCg                                               >  2:9220/1‑140 (MQ=255)
                                                                                   ggAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGt                                              <  1:1157660/140‑1 (MQ=255)
                                                                                         aTTAATATCAATTTGCACGGTAACGTTGTTGTTTTCTCCGGTGCCGTATACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCttt                                        <  2:522341/140‑1 (MQ=255)
                                                                                          ttAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACAACTGAACACTGGCGATACCACCGCACTCGTCCCtttt                                       >  1:522341/1‑140 (MQ=255)
                                                                                             atatCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGc                                    >  1:800331/1‑140 (MQ=255)
                                                                                              tatCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCGCTCGTCCCTTTTAGcc                                   >  1:965849/1‑140 (MQ=255)
                                                                                               atCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCt                                  <  1:1348748/140‑1 (MQ=255)
                                                                                                   aTTTGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAGCTACTGAACGCTGGCGGAACCACGGTACTCGGCCCTTTTCGCCTGCgg                              >  2:710257/1‑140 (MQ=255)
                                                                                                    tttGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCGCTCGTCCCTTTTAGCCTGCGGc                             >  1:705909/1‑140 (MQ=255)
                                                                                                    tttGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGc                             >  1:1461806/1‑140 (MQ=255)
                                                                                                    tttGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGc                             >  2:799547/1‑140 (MQ=255)
                                                                                                      tGCACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCGCTCGTCCCTTTTAGCCTGCGGCTg                           <  2:965849/140‑1 (MQ=255)
                                                                                                        cACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTg                                 >  2:43883/1‑132 (MQ=255)
                                                                                                        cACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTg                                 <  1:43883/132‑1 (MQ=255)
                                                                                                        cACGGTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACAGAACGCTGACGATACCACGGCACTCGTCCCATTTAGCCTGCGGCTGGa                         >  2:470009/1‑140 (MQ=255)
                                                                                                            gTAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGa                     <  1:236038/140‑1 (MQ=255)
                                                                                                             tAACGTTGTTGATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGAt                    <  1:1175072/140‑1 (MQ=255)
                                                                                                                  ttgttgATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGtcct               <  2:791392/140‑1 (MQ=255)
                                                                                                                     ttgATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCg            >  1:655612/1‑140 (MQ=255)
                                                                                                                       gATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCggg          >  1:39942/1‑140 (MQ=255)
                                                                                                                       gATTTCTCCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCggg          >  1:385972/1‑140 (MQ=255)
                                                                                                                         tttCTCCTGTGCCGTAAACACAGTGGATATTTATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCGGGtt        >  1:337268/1‑140 (MQ=255)
                                                                                                                         tttCTCCTGTGCCGTAAACACACCGTATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCGGGtt        <  1:740648/140‑1 (MQ=255)
                                                                                                                           tctcCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCGGGTTCa      >  1:803718/1‑140 (MQ=255)
                                                                                                                           tctcCTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCGGGTTCa      >  1:870712/1‑140 (MQ=255)
                                                                                                                               cTGTGCCGTAAACACAGCGGATCTTTATAAGACGGTATATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCGGGTTCAGTTg  <  1:587982/140‑1 (MQ=255)
                                                                                                                               cTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCGGGTTCAGTTg  <  2:655612/140‑1 (MQ=255)
                                                                                                                               cTGTGCCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCGGGTTCAGTTg  <  1:915618/140‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
AGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATATCAATTTGCACGGTAACGTTGTTGATTTCTCCTGTACCGTAAACACAGCGGATATTGATAAGACGGTAGATTTAGGCAGATGGCCTACGACACAACTACTGAACGCTGGCGATACCACGGCACTCGTCCCTTTTAGCCTGCGGCTGGAGGGATGTCCTCCGGGTTCAGTTG  >  NC_000913/563319‑563585

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