Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 524,391 G→T 5.3% P377P (CCG→CCT)  rhsD → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913524,3910GT 5.3% 463.7 / 12.2 132P377P (CCG→CCT) rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (62/63);  new base T (3/4);  total (65/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 133 total aligned reads displayed.

GTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCA  >  NC_000913/524245‑524539
                                                                                                                                                  |                                                                                                                                                    
gTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGc                                                                                                                                                    <  2:213025/149‑1 (MQ=255)
gTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGc                                                                                                                                                    >  1:983016/1‑149 (MQ=255)
     tGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGt                                                                                                                                               <  2:340509/149‑1 (MQ=255)
       aCCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTa                                                                                                                                             >  4:150495/1‑149 (MQ=255)
        ccGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                            <  4:1220224/149‑1 (MQ=255)
        ccGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                            <  3:1060632/149‑1 (MQ=255)
        ccGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                            <  1:577388/149‑1 (MQ=255)
          gcagcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCTGCCGGGTTCAgt                                                                                                                                          <  2:1477244/149‑2 (MQ=255)
          gcagcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGc                                                                                                                                          >  3:307026/1‑149 (MQ=255)
            agcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTa                                                                                                                                        >  1:1482085/1‑149 (MQ=255)
             gcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCTGCCGGGTTCAgttgc                                                                                                                                       <  1:1439495/149‑5 (MQ=255)
             gcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTAc                                                                                                                                       <  1:89876/149‑1 (MQ=255)
               aaTACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCg                                                                                                                                     >  3:1602880/1‑149 (MQ=255)
                    gCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                            >  4:76644/1‑137 (MQ=255)
                    gCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATc                                                                                                                                <  2:744785/149‑1 (MQ=255)
                         tGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCTCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTAt                                                                                                                           <  3:1060037/149‑1 (MQ=255)
                         tGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTAt                                                                                                                           <  1:67938/149‑1 (MQ=255)
                         tGCGCGCTTTCACGTATGACGCCCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTAt                                                                                                                           <  2:1084489/149‑1 (MQ=255)
                           cgcgcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGa                                                                                                                         <  3:1207317/149‑1 (MQ=255)
                            gcgcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAg                                                                                                                        >  2:1714505/1‑149 (MQ=255)
                             cgcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGc                                                                                                                       <  4:1361362/149‑1 (MQ=255)
                               cTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAg                                                                                                                     <  2:1614244/149‑1 (MQ=255)
                                tttCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAgg                                                                                                                    >  3:538161/1‑149 (MQ=255)
                                tttCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAgg                                                                                                                    >  3:1615131/1‑149 (MQ=255)
                                tttCACGTATGACGCGCACCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAgg                                                                                                                    >  2:613312/1‑149 (MQ=255)
                                  tCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGAc                                                                                                                  >  3:377044/1‑149 (MQ=255)
                                   cACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGAcc                                                                                                                 <  1:1392289/149‑1 (MQ=255)
                                   cACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGAcc                                                                                                                 >  3:1126471/1‑149 (MQ=255)
                                    aCGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCg                                                                                                                >  2:1134125/1‑149 (MQ=255)
                                    aCGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCg                                                                                                                >  4:225679/1‑149 (MQ=255)
                                          gACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtcac                                                                                                          >  2:148585/1‑149 (MQ=255)
                                           aCGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtcacc                                                                                                         <  4:644101/149‑1 (MQ=255)
                                                 agcaCCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtcaccgtcacc                                                                                                   >  1:1156292/1‑149 (MQ=255)
                                                  gcaCCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAtcaccgtcaccg                                                                                                  >  1:115858/1‑149 (MQ=255)
                                                   caCCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGa                                                                                                 <  3:1453482/149‑1 (MQ=255)
                                                      ccGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAg                                                                                              >  4:668044/1‑149 (MQ=255)
                                                      ccGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAg                                                                                              >  2:1133638/1‑149 (MQ=255)
                                                       cGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGc                                                                                             <  3:668044/149‑1 (MQ=255)
                                                       cGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGc                                                                                             <  1:1088826/149‑1 (MQ=255)
                                                        gggCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGcc                                                                                            <  2:47227/149‑1 (MQ=255)
                                                            cGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGaa                                                                                        <  2:1105379/149‑1 (MQ=255)
                                                               aTGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCg                                                                                     >  3:1669944/1‑149 (MQ=255)
                                                                    ggCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTg                                                                                <  1:625599/149‑1 (MQ=255)
                                                                     gcgcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGa                                                                               <  3:834211/149‑1 (MQ=255)
                                                                     gcgcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGa                                                                               <  2:280127/149‑1 (MQ=255)
                                                                     gcgcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCGTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGa                                                                               >  3:1389193/1‑149 (MQ=255)
                                                                       gcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAgg                                                                             >  1:556858/1‑149 (MQ=255)
                                                                       gcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGAGTTAAGCTAc                                                                                                                                       <  1:1427490/91‑1 (MQ=255)
                                                                       gcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGAGTTAAGCTAc                                                                                                                                       >  2:1427490/1‑91 (MQ=255)
                                                                             ttACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGtgctgc                                                                       <  4:1146442/149‑1 (MQ=255)
                                                                               aCGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCAt                                                                     <  4:1535812/149‑1 (MQ=255)
                                                                                   gggAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGc                                                                                   <  3:120572/131‑1 (MQ=255)
                                                                                        ggCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGc                                                            <  1:1378541/149‑1 (MQ=255)
                                                                                          ccGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgg                                                          >  4:70335/1‑149 (MQ=255)
                                                                                           cGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCggg                                                         <  1:1134125/149‑1 (MQ=255)
                                                                                           cGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCggg                                                         >  1:900226/1‑149 (MQ=255)
                                                                                            ggAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgggg                                                        <  3:1402946/149‑1 (MQ=255)
                                                                                              agaTGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGcc                                                      <  4:1295389/149‑1 (MQ=255)
                                                                                               gaTGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAggg                                                                                                                                                <  1:1040568/58‑1 (MQ=255)
                                                                                               gaTGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAggg                                                                                                                                                >  2:1040568/1‑58 (MQ=255)
                                                                                                aTGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCgg                                                    <  4:1389193/149‑1 (MQ=255)
                                                                                                 tGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCggg                                                   >  3:683812/1‑149 (MQ=255)
                                                                                                  gcgcTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGc                                                  <  2:1482085/149‑1 (MQ=255)
                                                                                                  gcgcTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGc                                                  <  4:1669944/149‑1 (MQ=255)
                                                                                                    gctaccgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTg                                                >  4:454783/1‑149 (MQ=255)
                                                                                                    gctaccgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTg                                                >  4:821250/1‑149 (MQ=255)
                                                                                                       accgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGaaa                                             <  1:344994/149‑1 (MQ=255)
                                                                                                         cgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACg                                           >  1:1719062/1‑149 (MQ=255)
                                                                                                          gctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAg                                                                 <  2:1383432/126‑1 (MQ=255)
                                                                                                           ctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACggg                                         >  2:1594905/1‑149 (MQ=255)
                                                                                                           ctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACggg                                         >  3:1303676/1‑149 (MQ=255)
                                                                                                            tacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggt                                        >  3:1192060/1‑149 (MQ=255)
                                                                                                            tacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggt                                        <  3:821250/149‑1 (MQ=255)
                                                                                                                      cGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGaaaaaa                              <  1:1414/149‑1 (MQ=255)
                                                                                                                       ggggCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGc                                                  >  3:249489/1‑128 (MQ=255)
                                                                                                                       ggggCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGc                                                  <  4:249489/128‑1 (MQ=255)
                                                                                                                       ggggCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAg                             >  1:655142/1‑149 (MQ=255)
                                                                                                                       ggggCGGGTGGTGACGCAACTGAACCCGGCAGGGTTAAGCTAGCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAg                             <  2:655142/149‑1 (MQ=255)
                                                                                                                        gggCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGa                            >  1:108801/1‑149 (MQ=255)
                                                                                                                           cGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACt                         >  3:706637/1‑149 (MQ=255)
                                                                                                                            gggTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTg                        >  4:1066432/1‑149 (MQ=255)
                                                                                                                             ggtggtGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTgg                       <  1:1714505/149‑1 (MQ=255)
                                                                                                                             ggtggtGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTgg                       <  1:59199/149‑1 (MQ=255)
                                                                                                                                      gCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCa              >  4:934196/1‑149 (MQ=255)
                                                                                                                                      gCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCa              >  1:395521/1‑149 (MQ=255)
                                                                                                                                      gCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCa              >  3:495971/1‑149 (MQ=255)
                                                                                                                                        aaCTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGc            <  3:780068/149‑1 (MQ=255)
                                                                                                                                         aCTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggtg                                       >  2:1384951/1‑121 (MQ=255)
                                                                                                                                         aCTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggtg                                       <  1:1384951/121‑1 (MQ=255)
                                                                                                                                         aCTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCg           >  3:1482523/1‑149 (MQ=255)
                                                                                                                                          cttaaCCCTGCCGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAGAAAGAACTGGCGGACGGCAGCGt          >  1:1597247/4‑149 (MQ=255)
                                                                                                                                          cTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGATGGCAGCGt          >  3:1534320/1‑149 (MQ=255)
                                                                                                                                           ttaaCCCTGCCGGGTTAAGCTACCGCTATATTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTc         >  4:1047684/3‑149 (MQ=255)
                                                                                                                                           tGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGaa                           <  4:532996/131‑1 (MQ=255)
                                                                                                                                           tGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGaa                           >  3:532996/1‑131 (MQ=255)
                                                                                                                                            gAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggt                                        <  2:365133/117‑1 (MQ=255)
                                                                                                                                            gAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggt                                        >  1:365133/1‑117 (MQ=255)
                                                                                                                                             aaCCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCAc       <  2:556858/149‑1 (MQ=255)
                                                                                                                                               cccGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCAcgc     >  1:93092/1‑149 (MQ=255)
                                                                                                                                                  ggCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCa  >  4:333913/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GTATATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCA  >  NC_000913/524245‑524539

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