Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,147,021 T→C 59.0% M174V (ATG→GTG)  alkA ← 3‑methyl‑adenine DNA glycosylase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,147,0210TC59.0% 164.1 / 300.4 239M174V (ATG→GTG) alkA3‑methyl‑adenine DNA glycosylase II
Reads supporting (aligned to +/- strand):  ref base T (45/53);  new base C (69/72);  total (114/125)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 241 total aligned reads displayed.

AAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGT  >  NC_000913/2146872‑2147169
                                                                                                                                                     |                                                                                                                                                    
  aataatTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATAc                                                                                                                                                    >  1:39521/1‑150 (MQ=255)
  aataatTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATAc                                                                                                                                                    >  2:2528534/1‑150 (MQ=255)
   ataatTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATAcc                                                                                                                                                   <  2:1473454/150‑1 (MQ=255)
       gtCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAAc                                                                                                                                               <  2:3467600/149‑1 (MQ=255)
           ccGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCtt                                                                                                                                           >  1:1351341/1‑150 (MQ=255)
           ccGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCtt                                                                                                                                           >  1:106233/1‑150 (MQ=255)
            cGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCttt                                                                                                                                          <  1:2072483/150‑1 (MQ=255)
              tCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTa                                                                                                                                        <  1:815443/150‑1 (MQ=255)
              tCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTa                                                                                                                                        <  1:248308/150‑1 (MQ=255)
                cAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAAt                                                                                                                                      <  2:1885139/150‑1 (MQ=255)
                     cccGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCTGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTg                                                                                                                                 <  1:448197/150‑1 (MQ=255)
                     cccGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTg                                                                                                                                 <  1:2595541/150‑1 (MQ=255)
                       cGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCg                                                                                                                               <  1:279392/150‑1 (MQ=255)
                         aTACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCggg                                                                                                                             >  2:3694491/1‑150 (MQ=255)
                              cGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTctgc                                                                                                                        >  1:1778892/1‑150 (MQ=255)
                              cGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTctgc                                                                                                                        >  1:322780/1‑150 (MQ=255)
                              cGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgc                                                                                                                        <  1:725578/150‑1 (MQ=255)
                                gAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTctgctg                                                                                                                      >  1:971329/1‑150 (MQ=255)
                                  aaaGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgctgct                                                                                                                    <  2:2809171/150‑1 (MQ=255)
                                  aaaGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTctgctgct                                                                                                                    <  2:2896991/150‑1 (MQ=255)
                                    aGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCtgctgctgc                                                                                                                  <  1:3351239/150‑1 (MQ=255)
                                     gTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGcc                                                                                                                 >  2:2267004/1‑150 (MQ=255)
                                        tGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGTTATGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCACGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGc                                                                                                              <  1:552492/150‑1 (MQ=255)
                                          cAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCg                                                                                                            >  1:2217115/1‑150 (MQ=255)
                                           aGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGc                                                                                                           <  1:431346/150‑1 (MQ=255)
                                            gCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCt                                                                                                          >  2:2352852/1‑150 (MQ=255)
                                                 ttCATCGCCTGCTCCACAGCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGc                                                                                                     <  2:3434169/150‑1 (MQ=255)
                                                    aTCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTc                                                                                                  <  2:539867/150‑1 (MQ=255)
                                                    aTCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTc                                                                                                  <  2:2699058/150‑1 (MQ=255)
                                                     tCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCg                                                                                                 >  1:2304330/1‑150 (MQ=255)
                                                     tCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCg                                                                                                 >  1:3480095/1‑150 (MQ=255)
                                                     tCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCg                                                                                                 >  1:2864765/1‑150 (MQ=255)
                                                     tCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCg                                                                                                 >  2:295989/1‑150 (MQ=255)
                                                     tCGCCTGCTCCACATCGCCCGGTAGGGTCATTGGTAAGGTGCCCTCCAGCGCCGCAGTTGCCAGATGAATCAGCGCCTCTGCCCGGTTCAACGGCACACCTAACGCTTGTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCg                                                                                                 >  2:1328251/1‑150 (MQ=255)
                                                      cGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCgg                                                                                                <  1:2272177/150‑1 (MQ=255)
                                                         cTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGaa                                                                                             <  2:1415724/150‑1 (MQ=255)
                                                         cTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGaa                                                                                             <  1:3555931/150‑1 (MQ=255)
                                                            cTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCa                                                                                          >  2:895249/1‑150 (MQ=255)
                                                            cTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCa                                                                                          >  2:3274697/1‑150 (MQ=255)
                                                             tCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAg                                                                                         <  1:897250/150‑1 (MQ=255)
                                                              ccACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGa                                                                                        <  1:1541055/150‑1 (MQ=255)
                                                               cacaTCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGat                                                                                       >  2:4191027/1‑150 (MQ=255)
                                                                acaTCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGata                                                                                      <  1:4171761/150‑1 (MQ=255)
                                                                    cGCCCGGTATTGTCATTGGTAAGGTGCCCTCCCGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACt                                                                                  <  2:605445/150‑1 (MQ=255)
                                                                           tATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGaaa                                                                           >  2:601487/1‑150 (MQ=255)
                                                                            aTTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGaaaa                                                                          <  1:614658/150‑1 (MQ=255)
                                                                            aTTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGaaaa                                                                          <  2:1090260/150‑1 (MQ=255)
                                                                             ttGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAat                                                                         >  1:841559/1‑150 (MQ=255)
                                                                             ttGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAat                                                                         <  2:2367932/150‑1 (MQ=255)
                                                                             ttGTCATTGGTAAGGTGCCCACCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAat                                                                         <  2:1178081/150‑1 (MQ=255)
                                                                                tCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcat                                                                      >  1:1679589/1‑150 (MQ=255)
                                                                                tCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcat                                                                      >  1:3130273/1‑150 (MQ=255)
                                                                                tCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcat                                                                      >  1:2729277/1‑150 (MQ=255)
                                                                                 gaTTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGTATCAGCGCCTCTGCCCGTTTCTACGGCATACCTAACGCTTTTATTGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGTAAatcatc                                                                     >  2:3654281/2‑150 (MQ=255)
                                                                                 cATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcatc                                                                     <  2:3346848/150‑1 (MQ=255)
                                                                                 cATTGGTAAGGTGCCCTCCAGCGCCGCATGTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAAGGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAatcatc                                                                     <  2:2310012/150‑1 (MQ=255)
                                                                                  aTTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATcc                                                                    >  2:361999/1‑150 (MQ=255)
                                                                                  aTTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATcc                                                                    >  1:1139549/1‑150 (MQ=255)
                                                                                     ggTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGc                                                                 <  2:3541790/150‑1 (MQ=255)
                                                                                          ggTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTc                                                            >  2:1842019/1‑150 (MQ=255)
                                                                                                cTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATa                                                      >  1:1149217/1‑150 (MQ=255)
                                                                                                   cAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAg                                                   <  1:2477760/150‑1 (MQ=255)
                                                                                                     gcgcCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCt                                                 >  2:60654/1‑150 (MQ=255)
                                                                                                        ccgcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCtgtg                                              <  2:4058809/150‑1 (MQ=255)
                                                                                                         cgcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGc                                             <  1:2961932/150‑1 (MQ=255)
                                                                                                         cgcATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGc                                             >  1:2394516/1‑150 (MQ=255)
                                                                                                             ttgggcAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCAct                                         >  2:93527/6‑150 (MQ=255)
                                                                                                                 ccAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTgg                                     >  2:1511169/1‑150 (MQ=255)
                                                                                                                  cAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCGGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGc                                    >  2:3016635/1‑150 (MQ=255)
                                                                                                                  cAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGc                                    >  2:2176835/1‑150 (MQ=255)
                                                                                                                  cAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGc                                    >  1:4221295/1‑150 (MQ=255)
                                                                                                                   aGATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCg                                   >  1:1291511/1‑150 (MQ=255)
                                                                                                                    gATGAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCgg                                  >  2:822349/1‑150 (MQ=255)
                                                                                                                       gAATCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATaa                                                     >  1:1256312/1‑128 (MQ=255)
                                                                                                                        aaTCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCaa                              <  2:4179856/150‑1 (MQ=255)
                                                                                                                         aTCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAt                             >  1:2487697/1‑150 (MQ=255)
                                                                                                                         aTCAGCGCCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGACATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAt                             >  1:3699215/1‑150 (MQ=255)
                                                                                                                           catcgcCTCTGACCGTTTCAAAGGCATACCTAAAGCTTTTAATGACTGCGGGTCTGATGCTGCCAGCCGCTGAGTCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAttt                           <  2:2027868/147‑1 (MQ=255)
                                                                                                                            aGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAtttt                          <  2:1149079/150‑1 (MQ=255)
                                                                                                                            aGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAtttt                          <  1:3579540/150‑1 (MQ=255)
                                                                                                                            aGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAtttt                          <  2:3947974/150‑1 (MQ=255)
                                                                                                                            aGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAtttt                          >  1:3797792/1‑150 (MQ=255)
                                                                                                                            aGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAAtttt                          <  2:1149227/150‑1 (MQ=255)
                                                                                                                             gcgcCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTg                         >  1:64753/1‑150 (MQ=255)
                                                                                                                             gcgcCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGCGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTg                         >  2:726653/1‑150 (MQ=255)
                                                                                                                              cgcCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgc                        <  2:2511655/150‑1 (MQ=255)
                                                                                                                               gcCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgcc                       <  1:4159998/150‑1 (MQ=255)
                                                                                                                               gcCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgcc                       >  2:1758582/1‑150 (MQ=255)
                                                                                                                               gcCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgcc                       >  1:400822/1‑150 (MQ=255)
                                                                                                                               gcCTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgcc                       <  1:2689200/150‑1 (MQ=255)
                                                                                                                                ccTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgccg                      <  2:2217152/150‑1 (MQ=255)
                                                                                                                                ccTCTGCCCGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTgccg                      >  1:3871207/1‑150 (MQ=255)
                                                                                                                                    tGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCAt                  >  1:3294827/1‑150 (MQ=255)
                                                                                                                                      cccGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGCTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCg                >  1:945271/1‑150 (MQ=255)
                                                                                                                                      cccGTTTCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCg                >  2:1023821/1‑150 (MQ=255)
                                                                                                                                            tCAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGc          >  2:408864/1‑150 (MQ=255)
                                                                                                                                             cAACGGCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCt         <  1:699241/150‑1 (MQ=255)
                                                                                                                                                cGGCATACCTACCGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCAc      <  1:1462340/150‑1 (MQ=255)
                                                                                                                                                 ggCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATGTTGCCGCCATCGCCACGCTCAcc     >  2:2443616/1‑150 (MQ=255)
                                                                                                                                                  gCACACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCa    >  1:2850087/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCAGATATACTCCGGAAAATCATCCAGCCGTTCGCCATAAAGCTGTGCCACTCTGGCGGTCAATTTTGCCGCCATCGCCACGCTCACCAGT  >  NC_000913/2146872‑2147169

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