Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 935,584 A→T 14.5% P787P (CCA→CCT ftsK → DNA translocase at septal ring sorting daughter chromsomes

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913935,5840AT14.5% 104.0 / 4.6 63P787P (CCA→CCTftsKDNA translocase at septal ring sorting daughter chromsomes
Reads supporting (aligned to +/- strand):  ref base A (34/17);  new base T (2/6);  total (39/24)
Fisher's exact test for biased strand distribution p-value = 4.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAC  >  NC_000913/935437‑935729
                                                                                                                                                   |                                                                                                                                                 
tGTACATCATCCTCATCATGCGTTTGCACGGGAGGTGGAATATCAGGAGCCTCACCAACCATTTCCGCCGCAGCATCAGTTTCAGCAGCCACAACTGCCGGCTGCGCCGCATGAACAATATCAGCAGCCGCAACAACCGGTTGCGCcaca                                                                                                                                                 <  2:351212/150‑1 (MQ=12)
  tACAGCAGCCGCAACACCCGCTTGCACCGCAGCAGCAATATCAGCACCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCAcagc                                                                                                                                               <  1:447342/150‑1 (MQ=255)
     agcagcCGCCCCCTCCGTTTGCTCCGCAGCAGGAAGATCTTCAGCCGCAACAACCCGTTCCGCCGCAGCATCAGATTCAGCAGCCACAACAGCCGGTTGCGACGCAGCAACAATATCAGCAGCCGCAACAACCGGTTGCGCCACagcagc                                                                                                                                            <  1:1353959/150‑1 (MQ=18)
     agcagcCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACagcagc                                                                                                                                            <  2:1526544/150‑1 (MQ=255)
        atccGCAACAACCTGTTGCACCGCAGCAGCAATATCAGCAGCAGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAat                                                                                                                                         <  2:910284/148‑1 (MQ=255)
           cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACACCCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACCGCCGGTTGCGCCGCATCCCCACTATCCGCAGCGGCAGCACCCGGGGGCgccgcagccgcagtggc                                                                                                                                      >  2:1762517/1‑145 (MQ=17)
           cGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGGATCAGCAGCCTCAACAGCAGGTTGAGCCGCATCCCCAGACACCGCAGCAGCAACAACCGGTGGCGCCGCTGCCGCaaattc                                                                                                                                      >  2:971184/1‑146 (MQ=17)
             caacaaCCGGTTGCACCGCAGCAGCAATATCAGAAGCCGCAACAACCAGTTCCGCCGCAGCCGCAATATCCGCAGCCACAACCGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCCGCAAGATcac                                                                                                                                    >  1:1486979/1‑149 (MQ=21)
                     ggTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACCGCCGGTTGCGCCTCAGCCACAATATCAGCAGCAGCAACAACCGGATGCGCCTCAGCAGAACTATCAGCAGCCGca                                                                                                                            >  1:1898256/1‑150 (MQ=25)
                     ggTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGca                                                                                                                            <  1:961754/150‑1 (MQ=255)
                         gCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTGCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCACAATACGCAGCATCAACAACCGATTGCGCCACAGCCGCAACTTCATCAGAAGcaacac                                                                                                                        >  2:773992/1‑149 (MQ=17)
                            ccGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCACTATCAGCAGCCACAACAGCCGGTTGCGCCGCAACCGCAATATCAGCAGCCGCGACAACGGGATGAGCGACATCAGCAGTATAAGCAGCCGcacgcccgg                                                                                                                     >  2:1500560/1‑143 (MQ=18)
                                  cagcaATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCACTATCCGCAGCAGCAACCGCCGGTGGCGCCACAGCGGCAAAGACGGCCGGCGCAACAACCGGATgcg                                                                                                               >  1:534914/1‑150 (MQ=18)
                                     caATATCAGCAGCCGCAACAACCCGGTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGcgccg                                                                                                            <  1:1166374/150‑1 (MQ=255)
                                            agcagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCACTATCAGCAGCATCAACAGCCGTTGGGGCCGCATCGGCAATGGCAGCAGCAGCAACAACCGGATACGCCgctgcagg                                                                                                     >  2:806839/1‑149 (MQ=18)
                                            agcagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCCCAATATCCGCAGCATCAACAACCGGTTGAGCCACATCAGCAATATCAGCAGCAGCAACAACCGAATGCGCCgcagccga                                                                                                     >  2:1588177/1‑147 (MQ=25)
                                            agcagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCCGCAGCCGCAACCACCGGTTGCGCCACAGCGGGACTATCCGCAGCCGCAACAACCGCAGGCGCCgcagccgc                                                                                                     >  1:1362245/1‑150 (MQ=25)
                                            agcagcCGCAACAACCAGATCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcagcagc                                                                                                     <  2:922118/150‑1 (MQ=255)
                                             gcagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCgcagcagca                                                                                                    <  1:885209/150‑1 (MQ=255)
                                              cagcCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCACAATCCGCAGCAGCAACAGCCGGTGCCGCCGCAGCCGCCAGATCAGCGGCGGCCACACCCGGATCCGCCgcggccgccg                                                                                                   >  1:1372854/1‑148 (MQ=16)
                                                       caacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCCGCCGCAACAACCGGATGCGCCGCAGCCGCAGTATcagaag                                                                                          >  1:247461/1‑150 (MQ=255)
                                                        aacCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCACCCACAATATCAGCAGCCGCAACAGCAGGTTGCGCCACAGCAGCAATATCCGCAGCCGCATCAACCGGTTCCGCCGCAGCCGCAGGATcagcacc                                                                                         >  1:1701813/1‑148 (MQ=21)
                                                         ccagaatCCGCCGCAGCCGCAGGATCAGCAGCCTCAGCAGCCGACTGCCCGGCAGGCACATGTTCAGCAGCCGCAACAACCGGTTGCTCCTCAGCAGCAATTTCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGcc                                                                                        <  2:83587/144‑1 (MQ=17)
                                                         acCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCAATATCCGCAGCCGCAACCACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGcc                                                                                        >  1:2025410/1‑150 (MQ=255)
                                                          ccAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACACCCGGTTGCGCCACAGCAGCAATATCAGCCGACGCAACCACCGGTTGCGCCGCAGCAGCAGTATCAGCAGccc                                                                                       >  1:1757325/1‑149 (MQ=255)
                                                             gTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCGCAGCCACAATATCCCCAGCAGCAACAACAGTTTGGGCGACAGCACGAAAAACAGCAGCGGCATCAACCGCTTCCGCCGCAGCCGAAGTATCAGCAGCGacaa                                                                                    >  2:627446/1‑150 (MQ=21)
                                                                      cagccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGAGCCACAGCCGCAATATCAGCAGCCGCAACCACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                                           >  1:1829040/1‑150 (MQ=255)
                                                                       agccgcagTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCGCACTATCCGCAGCCGCAACCACCGGTTGCGCCTCAGCAGCAATATCCGCCGCAGCAACACCCGGTTGCGCCGCAGCTGCAGGAGCAGCAGCGACAgccgccggagg                                                                          >  1:347234/1‑145 (MQ=255)
                                                                                ataagcagcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCTCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGGAGTATCAGCAGCAACAACCGCCAGTTGCGCcacacc                                                                 >  2:1034676/4‑148 (MQ=21)
                                                                                  ccgcagcCCCAACACCCGGTCGCGCCGCAGCCACATTAGCAGCCCCCGCAACTACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGGTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCg                                                               <  1:1500561/148‑1 (MQ=21)
                                                                                  cagcagCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACATCAGCAATATCAGCAGCAGCAACAACCGGTTGCGCCGCAGCAGCAGTAGCAGCAGCAACAACAGCCCGATGCGCCACAGCCg                                                               >  2:1181081/1‑150 (MQ=21)
                                                                                   agcagcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCCGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTTTCCGCAGCCACTACAGCCAGTGCCGCCACACCCGc                                                              >  2:1882985/1‑150 (MQ=21)
                                                                                    gcagcCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                                             <  1:1568583/150‑1 (MQ=255)
                                                                                              caGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTg                                                   >  2:1894085/1‑150 (MQ=255)
                                                                                                     ttgatctgcagcaACATTTTCAGCAGCCTCATCAACCTGTTGCTCCTCAGCATCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATc                                            <  2:1965603/143‑1 (MQ=18)
                                                                                                     ttGCGCCGCAGCAACATTATCAGCAGCCGCAACAACCGGTTGCGCGACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCTGTTGCGCCACAACCGCAGGATACCCTGCTTCATc                                            <  2:1862447/150‑1 (MQ=255)
                                                                                                       gcgcCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACCGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCCGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATccc                                          >  1:1582123/1‑149 (MQ=255)
                                                                                                         gccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCCCAACAGCCAGTTGCGCCACAACCGCCGGATACCCTGCTTCATCCGcg                                        >  1:32978/1‑149 (MQ=255)
                                                                                                          ccgcAGGAAGAATATCAGCAGCCGCAACTAGTGTCGGCTACTGATCATCAATATCAGCAGCCGCAACAACCGGTTGCTCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTg                                       <  1:137502/150‑1 (MQ=21)
                                                                                                           cgcAGCCACAAGCACAGCAGCCGCAACAACCGCCTGCTCCCCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt                                      <  1:87385/150‑1 (MQ=21)
                                                                                                            gcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGCTACCCTGCTTCATCCGCTggt                                     >  2:873849/1‑148 (MQ=255)
                                                                                                                 cacaATATCAGCAGCCGCAACAACCGGTTGCTCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCGTCCGCTGTTGATGc                                <  2:1481328/150‑1 (MQ=255)
                                                                                                                   caATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGGTGTTGAGGCgg                              >  2:52278/1‑149 (MQ=255)
                                                                                                                    aaTATCAGCAGCCGCAACAACCGGTTGCGCCTCAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa                             <  2:1398733/150‑1 (MQ=255)
                                                                                                                      tatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAAt                           <  2:1172370/150‑1 (MQ=255)
                                                                                                                       atCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATg                          >  1:1129595/1‑150 (MQ=255)
                                                                                                                        tCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATgg                         <  1:2008914/150‑1 (MQ=255)
                                                                                                                          agcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCg                       >  2:710561/1‑150 (MQ=255)
                                                                                                                             agcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGGTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACa                    <  2:1520126/150‑1 (MQ=255)
                                                                                                                              gcCTCAACAACCTGTTCCTACTGAGCATCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAg                   <  2:1982704/150‑1 (MQ=25)
                                                                                                                               ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCCGCAGTAGCAGCAGCCACAACCGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGc                  >  1:1302462/1‑150 (MQ=255)
                                                                                                                                 gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccg                >  2:527810/1‑150 (MQ=255)
                                                                                                                                 gCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccg                >  1:1176504/1‑150 (MQ=255)
                                                                                                                                   aacaacCCGACGCGCCACAGAAGAAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtc              <  2:568707/150‑1 (MQ=37)
                                                                                                                                    acaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACCACCGCCATTTGCGCCACAACCGCAGGATACCCTGCTTCATCGGCGGTTGATGAGTAATGGCGACAGccgtca             >  1:1260078/1‑149 (MQ=25)
                                                                                                                                    acaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtcc             >  2:1729853/1‑150 (MQ=255)
                                                                                                                                     caacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACCGCCAGTTGCGCCACGACCGCCGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGccgtccg            >  1:1532522/1‑150 (MQ=255)
                                                                                                                                       acCGGTTGAGAAACAGCAGCAATATCAGCAGCCTCAAAAACCGTTTGCGCCGCAGCAGCAGTATCAGAAGCCAAAACAGCCAGTTGCGACACAACCGCAGGATACACTGCTTCATACGATGTTTATGCGTAATGGCGACAGACGTCAGtt          >  2:1149924/1‑150 (MQ=255)
                                                                                                                                         cGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGCCCGCCGTCCGTTgg        >  2:1268839/1‑149 (MQ=255)
                                                                                                                                          ggTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa       >  2:1799509/1‑150 (MQ=255)
                                                                                                                                          ggTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa       >  1:500796/1‑150 (MQ=255)
                                                                                                                                          ggTTGCGCCACAGCAGCAATAACAGCAGCCGCAACAACCTCTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCa       <  1:131423/150‑1 (MQ=255)
                                                                                                                                           ggtGCTCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCTCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCAt      <  1:1226321/148‑1 (MQ=255)
                                                                                                                                               cgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc  >  1:1561042/1‑150 (MQ=255)
                                                                                                                                               cgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCATCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAc  >  2:1104161/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGTACAGCAGCCGCAACAACCGGTTGCACCGCAGCAGCAATATCAGCAGCCGCAACAACCAGTTCCGCCGCAGCCGCAGTATCAGCAGCCACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATGGCGACAGCCGTCCGTTGCATAAAC  >  NC_000913/935437‑935729

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