Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 525,739 G→C 5.9% E827Q (GAG→CAG) ‡ 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_000913525,7390GC 5.9% 446.4 / 15.3 135E827Q (GAG→CAG) ‡rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (65/62);  new base C (4/4);  total (69/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 = 2.79e-01

Only 100 of 140 total aligned reads displayed.

ACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCT  >  NC_000913/525590‑525887
                                                                                                                                                     |                                                                                                                                                    
acGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTg                                                                                                                                                      <  2:739261/150‑1 (MQ=255)
acGCATACAACGAGCAGGCGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTg                                                                                                                                                      <  1:1779061/150‑1 (MQ=255)
 cGCATCCAACGAGCAGGGGCTGGCAAACCGCGCCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTga                                                                                                                                                     <  2:1566368/150‑1 (MQ=255)
     tACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTAGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACg                                                                                                                                                 >  2:793184/1‑150 (MQ=255)
      aCAACGCGCCGGGGCCGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTAAGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGGCGAGTATACGCGGGACAGGCTGCACCGTGAGACgg                                                                                                                                                <  2:255538/150‑1 (MQ=18)
       cAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGt                                                                                                                                               >  2:672724/1‑150 (MQ=255)
         aCGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGCCGGTgc                                                                                                                                             >  1:1003953/1‑150 (MQ=255)
          cGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTgcg                                                                                                                                            >  1:911197/1‑150 (MQ=255)
          cGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTgcg                                                                                                                                            <  2:93814/150‑1 (MQ=255)
          cGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCGGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgca                                                                                                                                            <  2:235716/150‑2 (MQ=255)
           gAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcag                                                                                                                                           <  1:1032268/150‑3 (MQ=255)
            cgCAGGGGCTGGCAACCCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGGGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCa                                                                                                                                          <  1:1198129/149‑1 (MQ=255)
            cgCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCa                                                                                                                                          <  2:787047/149‑1 (MQ=255)
            cgCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGGGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcagc                                                                                                                                          <  1:655727/149‑4 (MQ=18)
            cgCAGGGGCGGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGGGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCa                                                                                                                                          <  1:491405/149‑1 (MQ=255)
            aGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCa                                                                                                                                          <  2:1176615/150‑1 (MQ=255)
            aGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCa                                                                                                                                          <  1:759100/150‑1 (MQ=255)
             gCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcagcc                                                                                                                                         <  2:941266/150‑5 (MQ=255)
              cAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGGGGAGTGGCTGACGTAGGGCAGCGGTTACCGGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGc                                                                                                                                        <  2:1643464/150‑1 (MQ=255)
                    cTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGTCGCTGGTCGAGTATACGCGGGACAGGCTGCCCCGGGAGACGGTGCGCAGCTTCGgc                                                                                                                                  >  1:153735/1‑150 (MQ=255)
                    cTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCCGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgc                                                                                                                                  >  2:123251/1‑150 (MQ=255)
                      ggCAAACAGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGCCGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgccc                                                                                                                                >  2:704295/1‑148 (MQ=255)
                       gCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgcagc                                                                                                                               >  2:1041307/1‑150 (MQ=255)
                       gCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGATGGGAGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCCCCGTGAGCCGGTGCGCAGCTGCGgcagc                                                                                                                               >  2:1083705/1‑150 (MQ=18)
                         aaaCCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCAt                                                                                                                             <  1:1422980/150‑1 (MQ=255)
                             cgcgTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                                                                                         >  1:951837/1‑150 (MQ=255)
                              gcgTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGCCGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcag                                                                                                                        >  1:768167/1‑150 (MQ=255)
                                     gCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGGCGAGTATACGCGGGACAGGCTGCACCGTGAGCCGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTaa                                                                                                                 >  1:1103182/1‑150 (MQ=255)
                                        ggACAGCCTGCCGCCGGTGGAGTGGCTGAGGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGc                                                                                                              >  2:539066/1‑150 (MQ=255)
                                        ggACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGc                                                                                                              >  2:670388/1‑150 (MQ=255)
                                         gACAGCCGGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGcc                                                                                                             <  1:1344021/150‑1 (MQ=255)
                                           cAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGc                                                                                                           <  2:5075/150‑1 (MQ=255)
                                           cAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGc                                                                                                           <  2:874378/150‑1 (MQ=255)
                                             gCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAt                                                                                                         >  1:774709/1‑150 (MQ=255)
                                              ccggccgccGGGGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATa                                                                                                        <  1:928802/147‑1 (MQ=255)
                                               cTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATAc                                                                                                       <  2:1579438/150‑1 (MQ=255)
                                                 gccgccGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGa                                                                                                     <  2:273595/150‑1 (MQ=255)
                                                 gccgccGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGa                                                                                                     >  2:903159/1‑150 (MQ=255)
                                                  ccgccgGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGaa                                                                                                    <  2:1056535/150‑1 (MQ=255)
                                                     ccgGTGGAGTGGCTGACGTATGGCAGCGGTAACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGAATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCATCTTAGGCAGCATGGCAGGCAGTACTGACGCATAAGAACTg                                                                                                 >  2:1194071/1‑150 (MQ=18)
                                                       ggTGGAGTGGCTGACGTATGGGAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGTCAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGAc                                                                                               <  1:1823333/150‑1 (MQ=255)
                                                          ggAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAg                                                                                            >  2:59922/1‑150 (MQ=255)
                                                           gAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGc                                                                                           <  1:1815242/150‑1 (MQ=255)
                                                             gTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGcac                                                                                         <  2:362598/150‑1 (MQ=255)
                                                             gTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGcac                                                                                         <  1:672724/150‑1 (MQ=255)
                                                               ggCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACAt                                                                                       >  2:1493586/1‑150 (MQ=255)
                                                                gCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATa                                                                                      >  2:599837/1‑150 (MQ=255)
                                                                  tGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATaca                                                                                    <  2:458101/150‑1 (MQ=255)
                                                                  ggACGTATGGCAGCGGTTACCTGGCGGGAAGGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATaca                                                                                    <  2:464828/149‑1 (MQ=255)
                                                                     cGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACAccc                                                                                 >  2:237099/1‑150 (MQ=255)
                                                                       tATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGCGACGCCCCTGGTCGAGTATACGCGGGACAGGCTGCAACGTGAGACGCTGCGCAGCTGCTGCAGCATGTCAGGCAGTAATTCCGCATACGAACTGACCAGCACATACAccccc                                                                               >  2:334482/1‑150 (MQ=18)
                                                                         tGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGc                                                                             <  2:554904/150‑1 (MQ=255)
                                                                          ggCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCCCATACACCCCCGCa                                                                            >  1:83999/1‑150 (MQ=255)
                                                                            cAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAgg                                                                          >  2:327423/1‑150 (MQ=255)
                                                                              gCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACTGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAgggc                                                                        >  1:480194/1‑148 (MQ=255)
                                                                                ggTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAg                                                                      <  2:547850/150‑1 (MQ=255)
                                                                                   tACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTa                                                                   >  2:417048/1‑150 (MQ=255)
                                                                                      cTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACag                                                                <  1:665608/150‑1 (MQ=255)
                                                                                          cGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGcc                                                            <  1:213014/150‑1 (MQ=255)
                                                                                             gAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACCCCCCCGCAGGCCAGTTAGAGAGCCAGc                                                         >  1:1708820/1‑150 (MQ=255)
                                                                                                tGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCAcc                                                      <  1:1660085/150‑1 (MQ=255)
                                                                                                 gAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGCCCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCAccg                                                     >  2:865888/1‑149 (MQ=255)
                                                                                                 gAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCAcca                                                     >  1:135092/1‑149 (MQ=255)
                                                                                                  aaGGTGGGCGGGACTGCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCTCATACGAACTGACCAGCACAAACTCCCCGGCAGGCCAGTTACAGAGCCAGCACCTg                                                    <  1:344506/150‑1 (MQ=17)
                                                                                                   aGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGa                                                   >  1:1378036/1‑150 (MQ=255)
                                                                                                    gCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGaa                                                  <  1:1028096/150‑1 (MQ=255)
                                                                                                    gCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGaa                                                  >  2:1301285/1‑150 (MQ=255)
                                                                                                     cTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACCGAGCCAGCACCTGAAc                                                 >  2:82283/1‑150 (MQ=255)
                                                                                                      tGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACa                                                >  2:1734230/1‑150 (MQ=255)
                                                                                                         gCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGcc                                             >  2:1216291/1‑150 (MQ=255)
                                                                                                          cGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCt                                            <  2:1783242/150‑1 (MQ=255)
                                                                                                          cGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCt                                            <  1:656304/150‑1 (MQ=255)
                                                                                                            ggACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTgg                                          <  2:1428366/150‑1 (MQ=255)
                                                                                                            ggACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTgg                                          >  2:623041/1‑150 (MQ=255)
                                                                                                               cgccgcTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGtat                                       >  2:284516/1‑150 (MQ=255)
                                                                                                                 ccgcTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTatat                                     >  2:1219551/1‑150 (MQ=255)
                                                                                                                     tGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGAcc                                 <  2:31685/150‑1 (MQ=255)
                                                                                                                      ggTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCg                                >  2:218401/1‑150 (MQ=255)
                                                                                                                       gTCGAGTATACGCGGGCCAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACCCCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGGATATGACCGt                               >  1:1233456/1‑150 (MQ=255)
                                                                                                                                aCGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACggg                      >  1:14449/1‑150 (MQ=255)
                                                                                                                                 cgcgGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGt                     >  2:1344119/1‑150 (MQ=255)
                                                                                                                                     ggACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAg                 >  2:1260170/1‑150 (MQ=255)
                                                                                                                                     ggACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAg                 <  2:396372/150‑1 (MQ=255)
                                                                                                                                      gACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGt                <  1:160663/150‑1 (MQ=255)
                                                                                                                                       aCAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTg               <  1:599837/150‑1 (MQ=255)
                                                                                                                                       aCAGGCGGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTg               <  2:1761134/150‑1 (MQ=255)
                                                                                                                                          ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACa            >  2:1747164/5‑150 (MQ=255)
                                                                                                                                          ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACa            >  1:1216076/5‑150 (MQ=255)
                                                                                                                                           tgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACaa           >  1:255069/4‑150 (MQ=255)
                                                                                                                                             tGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACg         <  2:420122/150‑1 (MQ=39)
                                                                                                                                              gCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACgg        >  1:841141/1‑150 (MQ=39)
                                                                                                                                                aCCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCg      <  2:482923/150‑1 (MQ=39)
                                                                                                                                                 ccGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGCACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGa     <  1:65980/150‑1 (MQ=255)
                                                                                                                                                 ccGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATAAGACCGTGACTACGGGTGGAGTGACAACGGCgc     >  1:898381/1‑149 (MQ=255)
                                                                                                                                                  cGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGAc    >  1:891953/1‑150 (MQ=39)
                                                                                                                                                  cGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGAc    <  1:623041/150‑1 (MQ=39)
                                                                                                                                                   gTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACACCGGCGAcc   <  1:1083705/150‑1 (MQ=255)
                                                                                                                                                   gTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGAcc   >  1:523353/1‑150 (MQ=39)
                                                                                                                                                    tGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCt  <  1:410934/150‑1 (MQ=39)
                                                                                                                                                    ggagaCGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCt  <  1:1493586/149‑1 (MQ=39)
                                                                                                                                                     |                                                                                                                                                    
ACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCT  >  NC_000913/525590‑525887

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