Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 524,394 A→C 5.3% A378A (GCA→GCC 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,3940AC 5.3% 464.9 / 15.3 131A378A (GCA→GCCrhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (62/62);  new base C (3/4);  total (65/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

Only 100 of 131 total aligned reads displayed.

TATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCA  >  NC_000913/524248‑524539
                                                                                                                                                  |                                                                                                                                                 
tatGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAgg                                                                                                                                                 <  2:1489307/149‑1 (MQ=255)
  tGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGt                                                                                                                                               <  2:340509/149‑1 (MQ=255)
   gACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGtt                                                                                                                                              <  2:110534/149‑1 (MQ=255)
    aCCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTa                                                                                                                                             >  4:150495/1‑149 (MQ=255)
     ccGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCTGCCGGGttca                                                                                                                                            <  3:1249012/149‑3 (MQ=255)
     ccGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                            <  1:577388/149‑1 (MQ=255)
     ccGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                            <  3:1060632/149‑1 (MQ=255)
     ccGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                            <  4:1220224/149‑1 (MQ=255)
       gcagcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCTGCCGGGttcagt                                                                                                                                          >  2:1439495/1‑145 (MQ=255)
       gcagcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCTGCCGGGTTCAgt                                                                                                                                          <  2:1477244/149‑2 (MQ=255)
         agcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTa                                                                                                                                        >  1:1482085/1‑149 (MQ=255)
          gcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCTGCCGGGTTCAgttgc                                                                                                                                       <  1:1439495/149‑5 (MQ=255)
          gcaATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCTGCCGGGTTCAgttgc                                                                                                                                       <  1:267522/149‑5 (MQ=255)
            aaTACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCg                                                                                                                                     >  3:1602880/1‑149 (MQ=255)
            aaTACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCg                                                                                                                                     >  3:1535812/1‑149 (MQ=255)
                 gCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                            >  4:76644/1‑137 (MQ=255)
                 gCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATc                                                                                                                                <  2:744785/149‑1 (MQ=255)
                 gCAGGTGCGCGCTGTAACGGATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTaa                                                                                                                                            <  3:76644/137‑1 (MQ=255)
                  cAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCt                                                                                                                               <  1:852656/149‑1 (MQ=255)
                     gTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTa                                                                                                                            <  4:1362220/149‑1 (MQ=255)
                      tGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTAt                                                                                                                           <  1:67938/149‑1 (MQ=255)
                        cgcgcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTAc                                                                                                                                       >  2:1142829/1‑135 (MQ=255)
                        cgcgcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGa                                                                                                                         <  3:1207317/149‑1 (MQ=255)
                           gcTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCa                                                                                                                      <  3:235029/149‑1 (MQ=255)
                            cTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAg                                                                                                                     <  2:1614244/149‑1 (MQ=255)
                             tttCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAgg                                                                                                                    >  1:47227/1‑149 (MQ=255)
                             tttCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAgg                                                                                                                    >  3:538161/1‑149 (MQ=255)
                             tttCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAgg                                                                                                                    >  3:1146442/1‑149 (MQ=255)
                              ttCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGa                                                                                                                   <  3:40848/149‑1 (MQ=255)
                               tCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGAc                                                                                                                  >  3:377044/1‑149 (MQ=255)
                                cACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGAcc                                                                                                                 >  3:1126471/1‑149 (MQ=255)
                                cACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGAcc                                                                                                                 <  1:1392289/149‑1 (MQ=255)
                                 aCGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCg                                                                                                                >  4:225679/1‑149 (MQ=255)
                                    tATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCAt                                                                                                             >  2:344994/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                                                                                              >  2:1133638/1‑149 (MQ=255)
                                                   ccGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAg                                                                                              >  4:668044/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)
                                                                  gcgcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGa                                                                               <  2:280127/149‑1 (MQ=255)
                                                                    gcACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGAGTTAAGCTAc                                                                                                                                       <  1:1427490/91‑1 (MQ=255)
                                                                            aCGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCAt                                                                     <  4:1535812/149‑1 (MQ=255)
                                                                                gggAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGc                                                                                   <  3:120572/131‑1 (MQ=255)
                                                                                gggAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGc                                                                                   >  4:120572/1‑131 (MQ=255)
                                                                                     ggCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAg                                                                              <  3:994334/131‑1 (MQ=255)
                                                                                     ggCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGc                                                            <  1:1378541/149‑1 (MQ=255)
                                                                                      gCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCg                                                           >  2:599116/1‑149 (MQ=255)
                                                                                       ccGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgg                                                          <  2:183939/149‑1 (MQ=255)
                                                                                       ccGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCgg                                                          >  4:70335/1‑149 (MQ=255)
                                                                                        cGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCggg                                                         >  1:900226/1‑149 (MQ=255)
                                                                                        cGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCggg                                                         <  1:1134125/149‑1 (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                                                  <  4:1669944/149‑1 (MQ=255)
                                                                                               gcgcTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGc                                                  <  2:1482085/149‑1 (MQ=255)
                                                                                                 gctaccgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTg                                                >  4:821250/1‑149 (MQ=255)
                                                                                                 gctaccgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTg                                                >  4:454783/1‑149 (MQ=255)
                                                                                                    accgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGaaa                                             <  1:599116/149‑1 (MQ=255)
                                                                                                      cgctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACg                                           >  1:1719062/1‑149 (MQ=255)
                                                                                                       gctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAg                                                                 <  2:1383432/126‑1 (MQ=255)
                                                                                                        ctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACggg                                         >  3:1303676/1‑149 (MQ=255)
                                                                                                        ctacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACggg                                         >  2:1594905/1‑149 (MQ=255)
                                                                                                         tacGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggt                                        <  3:821250/149‑1 (MQ=255)
                                                                                                          acgacgATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGggtg                                       <  1:932744/149‑1 (MQ=255)
                                                                                                                   cGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGaaaaaa                              <  1:1414/149‑1 (MQ=255)
                                                                                                                    ggggCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGc                                                  <  4:249489/128‑1 (MQ=255)
                                                                                                                    ggggCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGc                                                  >  3:249489/1‑128 (MQ=255)
                                                                                                                    ggggCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAg                             >  1:655142/1‑149 (MQ=255)
                                                                                                                    ggggCGGGTGGTGACGCAACTGAACCCGGCAGGGTTAAGCTAGCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAg                             <  2:655142/149‑1 (MQ=255)
                                                                                                                     gggCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGa                            >  1:108801/1‑149 (MQ=255)
                                                                                                                         gggTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTg                        >  4:1066432/1‑149 (MQ=255)
                                                                                                                          ggtggtGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTgg                       <  1:59199/149‑1 (MQ=255)
                                                                                                                          ggtggtGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTgg                       <  1:1714505/149‑1 (MQ=255)
                                                                                                                              gtggAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGa                   >  3:620452/1‑149 (MQ=255)
                                                                                                                                  aGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGTCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGc               >  3:1016127/1‑149 (MQ=255)
                                                                                                                                   gCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCTGACGGCa              >  3:1452860/1‑149 (MQ=255)
                                                                                                                                   gCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCa              >  3:495971/1‑149 (MQ=255)
                                                                                                                                   gCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCa              >  4:934196/1‑149 (MQ=255)
                                                                                                                                   gCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCa              >  1:395521/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)
                                                                                                                                       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)
                                                                                                                                         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)
                                                                                                                                                  |                                                                                                                                                 
TATGACCGCAGCAATACGCAGGTGCGCGCTTTCACGTATGACGCGCAGCACCCGGGCCGGATGGTGGCGCACCGTTACGCGGGAAGGCCGGAGATGCGCTACCGCTACGACGATACGGGGCGGGTGGTGGAGCAACTGAACCCGGCAGGGTTAAGCTACCGCTATCTTTATGAGCAGGACCGCATCACCGTCACCGACAGCCTGAACCGGCGTGAGGTGCTGCATACAGAAGGCGGGGCCGGGCTGAAACGGGTGGTGAAAAAAGAACTGGCGGACGGCAGCGTCACGCGCA  >  NC_000913/524248‑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: