Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 523,741 (G)6→7 coding (481/4281 nt) 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_000913523,7351.G97.1% 381.5 / 1.5 105coding (475/4281 nt)rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base . (2/1);  new base G (54/48);  total (56/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

GACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTGCCGCCGTAC  >  NC_000913/523585‑523878
                                                                                                                                                        |                                                                                                                                                
gaCGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCt‑gt                                                                                                                                                   <  1:205209/151‑1 (MQ=255)
   ggACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑ggg                                                                                                                                              <  4:449524/151‑1 (MQ=255)
     aCTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑ggggg                                                                                                                                            <  1:774403/151‑1 (MQ=255)
        gATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑gggggggc                                                                                                                                         <  2:211131/151‑3 (MQ=255)
        gATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑gggggggc                                                                                                                                         <  2:19502/151‑3 (MQ=255)
        gATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑ggggggg                                                                                                                                          <  1:1730349/150‑2 (MQ=255)
        gATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑gggggg                                                                                                                                           <  4:84087/149‑1 (MQ=255)
        gATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑gggggg                                                                                                                                           <  1:982057/149‑1 (MQ=255)
         aTACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑gggggggcg                                                                                                                                        <  3:490023/151‑4 (MQ=255)
         aTACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑gggggggcg                                                                                                                                        <  1:1451879/151‑4 (MQ=255)
           aCTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑gggggggcgct                                                                                                                                      <  4:1049673/151‑6 (MQ=255)
           aCTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑gggggcgcgc                                                                                                                                       <  4:1786056/150‑6 (MQ=255)
              cAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTg                                                                                                                                      <  1:1102258/149‑1 (MQ=255)
               aaCGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTgccg                                                                                                                                   <  1:597217/151‑1 (MQ=255)
               aaCGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTgccg                                                                                                                                   <  4:1073576/151‑1 (MQ=255)
                      aCGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGata                                                                                                                            <  1:549616/151‑1 (MQ=255)
                         gcgggcggCGCATTCACTTTGAGCCGCTGCGGGCGGGGGAGGCGGTGGACAGCCGCAGTGAGGCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGGTGGCGCGGGT‑GTG‑GGGGGGCGCTGCCGCCGGagagcg                                                                                                                         <  3:2601/151‑6 (MQ=255)
                          cgggcggAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGGG‑‑GGGGGGCGCTGCCGCCGGCTATcc                                                                                                                         >  1:105576/1‑150 (MQ=255)
                              cggAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGtt                                                                                                                     >  3:469770/1‑150 (MQ=255)
                                gAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCTGGG‑‑GGGGGGCGCTGCCGCCGGATATCCGGTTa                                                                                                                    >  1:470384/1‑149 (MQ=255)
                                    aTTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCg                                                                                                              >  1:787325/1‑151 (MQ=255)
                                     ttCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGGGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCCGCCGCCGGATATCCGGTTAAGCCCGc                                                                                                             >  3:344537/1‑151 (MQ=255)
                                         cTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATc                                                                                                          >  1:459251/1‑150 (MQ=255)
                                         cTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCt                                                                                                         >  4:1448173/1‑151 (MQ=255)
                                           ttGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCtt                                                                                                        <  3:533700/150‑1 (MQ=255)
                                           ttGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCt                                                                                                         <  1:1693128/149‑1 (MQ=255)
                                           ttGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCttt                                                                                                       <  4:361888/150‑1 (MQ=255)
                                            tGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCttt                                                                                                       <  1:1794764/150‑1 (MQ=255)
                                              aGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTAcc                                                                                                    >  2:58908/1‑151 (MQ=255)
                                              aGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGGGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTAcc                                                                                                    >  1:1149245/1‑151 (MQ=255)
                                               gCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCt                                                                                                   <  1:193646/151‑1 (MQ=255)
                                                ccGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTg                                                                                                  <  1:1817959/151‑1 (MQ=255)
                                                 cGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTgg                                                                                                 >  4:1326406/1‑151 (MQ=255)
                                                    tgctgcCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCg                                                                                               >  1:938029/1‑150 (MQ=255)
                                                    tgctgcCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCg                                                                                               >  3:541153/1‑150 (MQ=255)
                                                      ctgcCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGAcc                                                                                            >  1:1028869/1‑151 (MQ=255)
                                                       tgcCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGAcc                                                                                            >  2:81750/1‑150 (MQ=255)
                                                        gcCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGAcc                                                                                            <  4:770469/149‑1 (MQ=255)
                                                         ccGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGGGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAAc                                                                                         >  1:251592/1‑151 (MQ=255)
                                                           gggggAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAg                                                                                       <  2:447617/151‑1 (MQ=255)
                                                             gggAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAgcg                                                                                     <  1:1420603/151‑1 (MQ=255)
                                                               gAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAgcgc                                                                                    >  1:780582/1‑150 (MQ=255)
                                                                   cGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAgg                                                                               <  3:276790/151‑1 (MQ=255)
                                                                   cGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAgg                                                                               >  1:213176/1‑151 (MQ=255)
                                                                    ggTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAggg                                                                              >  4:452173/1‑151 (MQ=255)
                                                                    ggTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAggg                                                                              <  1:57656/151‑1 (MQ=255)
                                                                     gtgtACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAgggg                                                                             <  2:350423/151‑1 (MQ=255)
                                                                        tACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCg                                                                          >  3:418692/1‑151 (MQ=255)
                                                                         aCAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGGGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCAAAGGGGCCg                                                                          >  3:306646/1‑150 (MQ=255)
                                                                            gCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTGGGGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtgg                                                                       <  3:161835/151‑1 (MQ=255)
                                                                            gCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggt                                                                      <  4:721317/151‑1 (MQ=255)
                                                                             ccGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCgtggtg                                                                     >  1:704196/1‑151 (MQ=255)
                                                                                cAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGAt                                                                  >  4:1806132/1‑151 (MQ=255)
                                                                                  gTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATa                                                                 <  1:1287585/150‑1 (MQ=255)
                                                                                    gAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACt                                                               <  1:844477/150‑1 (MQ=255)
                                                                                    gAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTg                                                              <  4:336611/151‑1 (MQ=255)
                                                                                    gAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTg                                                              >  3:383360/1‑151 (MQ=255)
                                                                                     aGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTgg                                                             >  2:436883/1‑151 (MQ=255)
                                                                                     aGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTgg                                                             >  1:1003001/1‑151 (MQ=255)
                                                                                        cAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTgggg                                                           >  1:1433038/1‑150 (MQ=255)
                                                                                         aaTGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGggtg                                                         >  4:649986/1‑151 (MQ=255)
                                                                                              ggCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTg                                                    <  2:278565/151‑1 (MQ=255)
                                                                                              ggCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTg                                                    <  4:1661847/151‑1 (MQ=255)
                                                                                                cTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTg                                                    >  4:1194012/1‑149 (MQ=255)
                                                                                                 tGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGc                                                 >  1:1124179/1‑151 (MQ=255)
                                                                                                  ggTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcg                                                >  1:128523/1‑151 (MQ=255)
                                                                                                   gTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgg                                               <  4:470558/151‑1 (MQ=255)
                                                                                                   gTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgg                                               <  4:1626463/151‑1 (MQ=255)
                                                                                                   gTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgg                                               >  1:1529872/1‑151 (MQ=255)
                                                                                                       gcgGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgc                                           >  2:15849/1‑151 (MQ=255)
                                                                                                       gcgGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgc                                           >  1:1020749/1‑151 (MQ=255)
                                                                                                         ggtggtAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgccg                                         >  4:1024789/1‑151 (MQ=255)
                                                                                                         ggtggtAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgccg                                         <  1:1610937/151‑1 (MQ=255)
                                                                                                           tggtAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgccggg                                       >  4:1713561/1‑151 (MQ=255)
                                                                                                            ggtAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgccgggt                                      <  2:418706/151‑1 (MQ=255)
                                                                                                              tAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGcgggtgccgggtgc                                    >  3:424536/1‑151 (MQ=255)
                                                                                                                 ggCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCCGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGa                                 <  1:495419/151‑1 (MQ=255)
                                                                                                                  gcagcaCAGGCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTg                                  >  3:151947/1‑149 (MQ=255)
                                                                                                                   cagcaCAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAgg                               >  4:1493040/1‑151 (MQ=255)
                                                                                                                   cagcaCAGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAgg                               >  1:125906/1‑151 (MQ=255)
                                                                                                                         aGCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGAGCTGAGGACGTAc                         >  1:319004/1‑151 (MQ=255)
                                                                                                                          gCCGGACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACt                        <  2:332287/151‑1 (MQ=255)
                                                                                                                             ggACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGcc                     >  4:1719264/1‑151 (MQ=255)
                                                                                                                              gACGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGcc                     >  1:1013733/1‑150 (MQ=255)
                                                                                                                                cGGCCATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAgc                  >  1:994857/1‑151 (MQ=255)
                                                                                                                                   ccATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGcgcc               >  2:412189/1‑151 (MQ=255)
                                                                                                                                   ccATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACAAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGCCTGAGCGGGGGCCGGGTGCTGAGGACGTACTGCCAAcgcc               >  3:302883/1‑151 (MQ=255)
                                                                                                                                    cATACGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGcgccg              >  4:1379076/1‑151 (MQ=255)
                                                                                                                                       aCGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCt            <  4:1188127/150‑1 (MQ=255)
                                                                                                                                       aCGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTg           <  1:97209/151‑1 (MQ=255)
                                                                                                                                        cGCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTgc          <  4:1260185/151‑1 (MQ=255)
                                                                                                                                         gCTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTgcc         <  1:1450372/151‑1 (MQ=255)
                                                                                                                                          cTGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTgccg        >  4:599752/1‑151 (MQ=255)
                                                                                                                                           tGGCGCGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTgccgc       >  4:192689/1‑151 (MQ=255)
                                                                                                                                             gcgcGGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTgccgcct     <  4:539720/151‑2 (MQ=255)
                                                                                                                                              cgcgGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTGCCGCCGt    <  2:139735/151‑1 (MQ=255)
                                                                                                                                               gcgGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTGCCGCCGt    >  4:166085/1‑150 (MQ=255)
                                                                                                                                               gcgGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTGCCGCCGTa   <  1:1715871/151‑1 (MQ=255)
                                                                                                                                                cgGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTGCCGCCGTAc  >  4:1677212/1‑151 (MQ=255)
                                                                                                                                                cgGCT‑GTG‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTGCCGCCGTAc  >  1:1036203/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                
GACGGACTGATACTCAACGACAACGGCGGGCGGAGCATTCACTTTGAGCCGCTGCTGCCGGGGGAGGCGGTGTACAGCCGCAGTGAGTCAATGTGGCTGGTGCGCGGTGGTAAGGCAGCACAGCCGGACGGCCATACGCTGGCGCGGCT‑GT‑‑GGGGGGCGCTGCCGCCGGATATCCGGTTAAGCCCGCATCTTTACCTGGCGACCAACAGCGCACAGGGGCCGTGGTGGATACTGGGGTGGTCTGAGCGGGTGCCGGGTGCTGAGGACGTACTGCCAGCGCCGCTGCCGCCGTAC  >  NC_000913/523585‑523878

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