Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,622,422 C→G 5.9% A1077A (GCC→GCG rhsB → 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_0009133,622,4220CG 5.9% 3608.0 / 178.6 1048A1077A (GCC→GCGrhsBRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (556/430);  new base G (28/33);  total (585/463)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 1097 total aligned reads displayed.

CGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGA  >  NC_000913/3622272‑3622571
                                                                                                                                                      |                                                                                                                                                     
 ggCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGcc                                                                                                                                                       >  1:9557673/1‑150 (MQ=255)
  gCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGccc                                                                                                                                                      <  1:8163929/150‑1 (MQ=255)
    tGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCtt                                                                                                                                                    >  1:5711644/1‑150 (MQ=255)
    tGGGACGGCGACAGGCTGACCACGATACAGAACGACAGCAGCCGCATCCAGACCAATTATCAGCCGGGTAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGCATTACCtt                                                                                                                                                    >  1:11448897/1‑150 (MQ=14)
     gggACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTc                                                                                                                                                   <  2:1690390/150‑1 (MQ=255)
      ggACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTca                                                                                                                                                  <  1:119821/150‑1 (MQ=255)
       gACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcag                                                                                                                                                 <  2:4990947/150‑1 (MQ=255)
          ggCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTcagcag                                                                                                                                              >  1:6080787/1‑150 (MQ=255)
            cGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTc                                                                                                                                            >  2:2123520/1‑150 (MQ=255)
             gACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTcc                                                                                                                                           >  2:2345094/1‑150 (MQ=255)
             gACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTcc                                                                                                                                           >  2:5826634/1‑150 (MQ=255)
              aCCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCcg                                                                                                                                          >  2:1206809/1‑150 (MQ=255)
                  gCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCcggcgg                                                                                                                                      >  2:3861037/1‑150 (MQ=255)
                      aCCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGaa                                                                                                                                  >  2:4316411/1‑150 (MQ=255)
                       ccACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAg                                                                                                                                 <  2:11830018/150‑1 (MQ=255)
                        cACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGa                                                                                                                                >  1:11589107/1‑150 (MQ=255)
                        cACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGa                                                                                                                                >  1:10751857/1‑150 (MQ=255)
                        cACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGa                                                                                                                                >  1:5170938/1‑150 (MQ=255)
                         aCGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGAc                                                                                                                               >  1:2505376/1‑150 (MQ=255)
                         aCGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGAc                                                                                                                               >  1:4192431/1‑150 (MQ=255)
                           gATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACgg                                                                                                                             >  1:2398038/1‑150 (MQ=255)
                            aTACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGt                                                                                                                            <  2:1506378/150‑1 (MQ=255)
                              aCAGAACTACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCATGGTCGAAAATGCCACCGGATATCTGTCGAAAAAGCAGCTCAGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAATACGCTgg                                                                                                                          >  2:10927237/1‑150 (MQ=14)
                                  aaCGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAgt                                                                                                                      >  1:2000301/1‑150 (MQ=255)
                                     gACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGtgtg                                                                                                                   <  2:6525278/150‑1 (MQ=255)
                                     gACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGtgtg                                                                                                                   <  1:6748132/150‑1 (MQ=255)
                                      aCAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTgtgg                                                                                                                  >  1:1738307/1‑150 (MQ=255)
                                        aGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTgtggtg                                                                                                                <  2:7570179/150‑1 (MQ=40)
                                        aGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTgtggtg                                                                                                                >  1:8375982/1‑150 (MQ=40)
                                        aGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTgtggtg                                                                                                                >  2:2987156/1‑150 (MQ=40)
                                          gAGCCGCATACAGACGATTTATCAGCCGGGGATCATCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGtt                                                                                                              >  1:3250602/1‑150 (MQ=255)
                                               gCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCccgc                                                                                                         >  1:6560114/1‑150 (MQ=37)
                                                 aTCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCccgccg                                                                                                       >  2:9043056/1‑150 (MQ=37)
                                                      gACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCt                                                                                                  >  1:1229912/1‑150 (MQ=37)
                                                      gACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCt                                                                                                  >  1:9151113/1‑150 (MQ=37)
                                                       aCGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTg                                                                                                 <  1:10091537/150‑1 (MQ=37)
                                                         gATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGt                                                                                               <  1:4728443/150‑1 (MQ=21)
                                                              aTCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGa                                                                                          >  2:6761051/1‑150 (MQ=37)
                                                               tCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGAt                                                                                         <  1:5433153/150‑1 (MQ=37)
                                                                  gCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCt                                                                                      >  1:5849010/1‑150 (MQ=37)
                                                                   ccGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTc                                                                                     <  2:10263294/150‑1 (MQ=37)
                                                                       ggAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGAcc                                                                                 >  1:11127303/1‑150 (MQ=37)
                                                                            ttCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTg                                                                            >  1:10975668/1‑150 (MQ=37)
                                                                             tCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTgg                                                                           >  1:6911764/1‑150 (MQ=37)
                                                                               aCGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGaa                                                                         >  2:6182522/1‑150 (MQ=21)
                                                                                cGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGaaa                                                                        <  1:3485118/150‑1 (MQ=21)
                                                                                cGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGaaa                                                                        <  2:7470268/150‑1 (MQ=37)
                                                                                cGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGaaa                                                                        <  1:9917710/150‑1 (MQ=37)
                                                                                   cACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTg                                                                     >  1:9283876/1‑150 (MQ=37)
                                                                                    aCTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGa                                                                    >  2:11448807/1‑150 (MQ=37)
                                                                                    aCTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGa                                                                    >  1:10244246/1‑150 (MQ=37)
                                                                                     cTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGaa                                                                   >  1:11293601/1‑150 (MQ=37)
                                                                                      tcatcaGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCGCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGaaa                                                                  <  2:1044641/150‑1 (MQ=21)
                                                                                        atcaGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATc                                                                >  2:190324/1‑150 (MQ=37)
                                                                                         tcaGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATcc                                                               >  1:2155200/1‑150 (MQ=37)
                                                                                         tcaGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATcc                                                               >  1:1455717/1‑150 (MQ=37)
                                                                                          caGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCt                                                              <  1:2856665/150‑1 (MQ=37)
                                                                                          caGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCt                                                              <  2:8496302/150‑1 (MQ=37)
                                                                                           aGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTg                                                             <  1:8966605/150‑1 (MQ=37)
                                                                                            gggTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                            <  2:1219928/150‑1 (MQ=37)
                                                                                            gggTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTgg                                                            >  2:404918/1‑150 (MQ=37)
                                                                                              gTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCt                                                          >  2:11121117/1‑150 (MQ=37)
                                                                                               tCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTg                                                         >  1:7433722/1‑150 (MQ=35)
                                                                                                  aaaCTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGAcc                                                      >  1:7153508/1‑150 (MQ=34)
                                                                                                    aCTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCgg                                                    >  1:7655802/1‑150 (MQ=34)
                                                                                                          aCCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtgagt                                              >  2:10696721/1‑150 (MQ=32)
                                                                                                            cGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGgtgagtga                                            <  2:5629988/150‑1 (MQ=32)
                                                                                                             ggTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGtgagtgag                                           <  2:6313655/150‑1 (MQ=32)
                                                                                                              gTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAgg                                          <  2:5235426/150‑1 (MQ=32)
                                                                                                              gTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAgg                                          <  2:2712917/150‑1 (MQ=32)
                                                                                                                    tGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgcc                                    >  1:10556234/1‑150 (MQ=32)
                                                                                                                    tGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgcc                                    >  2:473064/1‑150 (MQ=32)
                                                                                                                    tGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgcc                                    <  1:8725582/150‑1 (MQ=32)
                                                                                                                      gCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAgccgc                                  >  2:727019/1‑150 (MQ=32)
                                                                                                                        gAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                                >  2:3419744/1‑150 (MQ=32)
                                                                                                                        gAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGccgccg                                >  2:2539948/1‑150 (MQ=32)
                                                                                                                            aCGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctgg                            >  1:4047231/1‑150 (MQ=32)
                                                                                                                              gCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggct                          >  2:9548180/1‑150 (MQ=32)
                                                                                                                              gCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggct                          >  2:4517364/1‑150 (MQ=32)
                                                                                                                               cAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCTGTGTGGTGTTCCCGCGGGTGCTGGTGCAGATGCGCGACCGGCTGTAAAGTTAAATCCTGGCTGACCGGGTGAGTGAGTAAAGCCGCCgctggctg                         <  1:3730786/150‑1 (MQ=11)
                                                                                                                               cAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCgctggctg                         <  2:6048924/150‑1 (MQ=32)
                                                                                                                                aGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGctggctgg                        <  2:3263739/150‑1 (MQ=32)
                                                                                                                                aGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGctggctgg                        <  2:872341/150‑1 (MQ=32)
                                                                                                                                    ccgcAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATc                    <  2:7076195/150‑1 (MQ=32)
                                                                                                                                     cgcAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCg                   >  1:6258494/1‑150 (MQ=32)
                                                                                                                                     cgcAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCg                   >  2:8890829/1‑150 (MQ=32)
                                                                                                                                      gcAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGt                  <  2:9309589/150‑1 (MQ=32)
                                                                                                                                       cAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTg                 >  1:4954295/1‑150 (MQ=32)
                                                                                                                                          ccTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAACTCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg              <  2:8309590/150‑1 (MQ=17)
                                                                                                                                          ccTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCgg              >  2:11615043/1‑150 (MQ=32)
                                                                                                                                           cTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGc             <  1:7852730/150‑1 (MQ=32)
                                                                                                                                            tGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGcc            <  2:8108088/150‑1 (MQ=32)
                                                                                                                                              gCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTg          >  1:11883252/1‑150 (MQ=32)
                                                                                                                                               cGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGa         >  1:4466184/1‑150 (MQ=32)
                                                                                                                                               cGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGa         <  2:1338942/150‑1 (MQ=32)
                                                                                                                                                 gcggCCCTGCCGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGGTCCCGCCGGTGCGGGTGCAGCTGCGCGACCGGCTGGAACGTGCACTCCGGGCGGACCGGGTGAGGGCGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGAcc       <  2:9725975/147‑1 (MQ=11)
                                                                                                                                                 gATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCAGCGTGCGGCCTGAcc       >  2:7878136/1‑150 (MQ=17)
                                                                                                                                                      ccTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGa  <  2:7734412/150‑1 (MQ=11)
                                                                                                                                                      ccTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGa  <  2:5395930/150‑1 (MQ=11)
                                                                                                                                                      ccTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGa  >  2:3728744/1‑150 (MQ=11)
                                                                                                                                                      |                                                                                                                                                     
CGGCTGGGACGGCGACCGGCTGACCACGATACAGAACGACAGGAGCCGCATCCAGACGATTTATCAGCCGGGGAGCTTCACGCCACTCATCAGGGTCGAAACTGCCACCGGTGAGCTGGCGAAAACGCAGCGCCGCAGCCTGGCGGATGCCCTTCAGCAGTCCGGCGGCGAAGACGGTGGCAGTGTGGTGTTCCCGCCGGTGCTGGTGCAGATGCTCGACCGGCTGGAAAGTGAAATCCTGGCTGACCGGGTGAGTGAGGAAAGCCGCCGCTGGCTGGCATCGTGCGGCCTGACCGTGGA  >  NC_000913/3622272‑3622571

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