Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,189 A→G 6.0% A121A (GCA→GCG yhcG → DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,1890AG 6.0% 859.0 / 30.3 268A121A (GCA→GCGyhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base A (142/107);  new base G (7/9);  total (150/118)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.52e-01

Only 100 of 269 total aligned reads displayed.

AAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTC  >  NC_000913/3368040‑3368337
                                                                                                                                                     |                                                                                                                                                    
aaGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCa                                                                                                                                                      <  2:3331145/150‑1 (MQ=255)
  ggATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGa                                                                                                                                                    <  1:4215193/150‑1 (MQ=255)
   gATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCTGCAGACGGTGTCTGCGGaa                                                                                                                                                   <  1:3219598/150‑1 (MQ=255)
   gATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGaa                                                                                                                                                   <  2:1161115/150‑1 (MQ=255)
      ttATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCTGCAGACGGTGTCTGCGGAATTc                                                                                                                                                <  1:1034318/150‑1 (MQ=255)
        ttGTCTTAGGTATAAGCGTGGGTTCTCTGCAAAACACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGTAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTccc                                                                                                                                              <  2:2915889/149‑1 (MQ=255)
         tGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCa                                                                                                                                             <  2:1941835/150‑1 (MQ=255)
          gTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCAt                                                                                                                                            <  1:4321457/150‑1 (MQ=255)
            cTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCg                                                                                                                                          >  2:4840578/1‑150 (MQ=255)
            cTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTACAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCTGAATTCCCATCg                                                                                                                                          <  1:4540393/150‑1 (MQ=255)
             tAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCgg                                                                                                                                         >  2:1544123/1‑150 (MQ=255)
                 gTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGCGGCTTTTTTACATATTTTTTCAACATGTTTAAATTCACCAGACATTGTCTGGCTAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCaaa                                                                                                                                     <  1:4039977/150‑1 (MQ=255)
                 gTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCaaa                                                                                                                                     >  1:5139943/1‑150 (MQ=255)
                      aGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGAttt                                                                                                                                <  2:4633171/150‑1 (MQ=255)
                      aGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGAttt                                                                                                                                <  2:4596077/150‑1 (MQ=255)
                       gCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTg                                                                                                                               >  2:5001009/1‑150 (MQ=255)
                       gCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTg                                                                                                                               >  1:3041245/1‑150 (MQ=255)
                       gCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATAGGGAAATCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAACATTTg                                                                                                                               >  1:4201938/1‑150 (MQ=255)
                       gCGTGGGCTCCCTGCAACAAACGTCCGCCCCCTGAGGCATTTTTACCTCTTTTTTCAACAGGGGGACAGTCACCCGACACTGTCGGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTg                                                                                                                               <  2:2348859/150‑1 (MQ=255)
                        cGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGATTTCCCATCGGCAAAGATTTgg                                                                                                                              <  1:3795993/150‑1 (MQ=255)
                        cGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTgg                                                                                                                              <  2:1221126/150‑1 (MQ=255)
                        cGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGATTTTTTACCTCTTTTTTCAACATTTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTgg                                                                                                                              <  1:1145137/150‑1 (MQ=255)
                          tGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCa                                                                                                                            <  1:3273685/150‑1 (MQ=255)
                          ggggTTCTCTGCAAAAAACTTACGCCAAATGAGGCTGTTTTACCTCTTTTTTAAACATGTTGAAATTCACCAGACAATGTCTGGCGACTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGcc                                                                                                                            >  1:2760495/2‑149 (MQ=255)
                            ggTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGa                                                                                                                          >  2:4439830/1‑150 (MQ=255)
                              ttCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGCCg                                                                                                                        >  2:4042940/1‑150 (MQ=255)
                              ttCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACg                                                                                                                        >  2:4577140/1‑150 (MQ=255)
                                ctctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAcagc                                                                                                                      <  1:2499122/150‑4 (MQ=255)
                                ctctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAcagc                                                                                                                      <  2:2686336/150‑4 (MQ=255)
                                 tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTg                                                                                                                     >  1:2856860/1‑150 (MQ=255)
                                  ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCCGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt                                                                                                                    >  1:1362455/1‑150 (MQ=255)
                                  ctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt                                                                                                                    <  1:4429352/150‑1 (MQ=255)
                                    gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                  <  1:2350415/150‑1 (MQ=255)
                                      aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                >  1:3134870/1‑150 (MQ=255)
                                      aaaaaaCTTACGCCAAATGAGGATTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                <  2:1967131/150‑1 (MQ=255)
                                       aaaaaaTTACGCCAAATGAGGCTTTTTTACCTATTTTTTCAACTTGTTGAAATTCACCAGACATTGTCTGGCTAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                               <  1:3576194/150‑1 (MQ=255)
                                         aaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTTAAATTCAACAGACAATGTCTGGCGAATTAACACCATTTGGAATTCCTCAGACACCGTCTGCAGATTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTTCCaa                                                                                                             <  2:2364863/150‑1 (MQ=255)
                                          aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCaaa                                                                                                            <  1:752214/150‑1 (MQ=255)
                                           aCTTACCCCAAATGAGGCTTTTTTACCTCTTTTTTAAACATGTTGAAATACACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCATAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAAt                                                                                                           <  1:1721505/150‑1 (MQ=255)
                                            cTTACGACAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCAAATc                                                                                                          <  1:3353835/150‑1 (MQ=255)
                                                 gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGGAGTCGCTGTGTGCGAAATCtttgc                                                                                                     >  1:1098190/1‑148 (MQ=255)
                                                 gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                     >  2:1751502/1‑150 (MQ=255)
                                                 gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                     >  2:4693063/1‑150 (MQ=255)
                                                                  tACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAcc                                                                                    >  2:3561961/1‑150 (MQ=255)
                                                                     ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                 >  2:4487452/1‑150 (MQ=255)
                                                                     ctctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                 >  1:2977167/1‑150 (MQ=255)
                                                                      tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACg                                                                                >  1:966404/1‑150 (MQ=255)
                                                                      gtgaatttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACg                                                                                >  2:267282/6‑150 (MQ=255)
                                                                       ctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCGGACTTCCCATCGGCCAAGATTTGGCAGACGCTGGCTGCCAAAGCTTTCCCCTTGCCCAGGGCCACCTACGt                                                                               >  1:4935986/1‑150 (MQ=255)
                                                                              cAACATGTTGAAATTCACCAGACAATGTCTGGCGATTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTa                                                                        <  2:1713878/150‑1 (MQ=255)
                                                                              cAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTa                                                                        <  1:3848689/150‑1 (MQ=255)
                                                                               aaCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTAc                                                                       >  1:4486306/1‑150 (MQ=255)
                                                                               aaCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTAc                                                                       >  2:2117799/1‑150 (MQ=255)
                                                                                aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACt                                                                      >  2:788842/1‑150 (MQ=255)
                                                                                    gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                  >  1:3620082/1‑150 (MQ=255)
                                                                                    gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                  <  1:1178448/150‑1 (MQ=255)
                                                                                    gTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTcc                                                                  <  1:152426/150‑1 (MQ=255)
                                                                                     ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCAACGTCCGTTTACTTTCCg                                                                 >  1:5113770/1‑150 (MQ=255)
                                                                                     ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATATTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                 >  1:4789151/1‑150 (MQ=255)
                                                                                     tggAAATTCACCAGACAATGGCTGGCGAATTACCACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                 <  1:2159878/148‑1 (MQ=255)
                                                                                      tGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAGTCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGt                                                                >  2:1051704/1‑150 (MQ=255)
                                                                                         aaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaa                                                             <  2:4812662/150‑1 (MQ=255)
                                                                                         aaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaa                                                             <  1:4263461/150‑1 (MQ=255)
                                                                                           ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCCGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                           >  1:3723985/1‑150 (MQ=255)
                                                                                           ttCGCCAGACATTGTCTGGTGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                           >  1:693598/1‑150 (MQ=255)
                                                                                            tCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaaa                                                          >  1:2025209/1‑150 (MQ=255)
                                                                                            tCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaaa                                                          >  1:1624394/1‑150 (MQ=255)
                                                                                              aCCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCCGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTGTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGGTACAAACg                                                        >  2:4942871/1‑150 (MQ=255)
                                                                                              aCCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACg                                                        >  2:1465705/1‑150 (MQ=255)
                                                                                               ccAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                       <  2:868449/150‑1 (MQ=255)
                                                                                               ccAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                       <  2:2361583/150‑1 (MQ=255)
                                                                                               ccAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACTTTCCGTTAAAAACGc                                                       >  1:2313310/1‑150 (MQ=255)
                                                                                                  gACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGa                                                    >  2:294821/1‑150 (MQ=255)
                                                                                                   aCAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGAc                                                   >  2:5074511/1‑150 (MQ=255)
                                                                                                     aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGc                                                 >  1:4768444/1‑150 (MQ=255)
                                                                                                           tGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGAAGCTGTGTGGCCAAACTTTCCCCATGCCCAGGTCCACCAACGTCCGGTTAATTTACGGTAAAAACCCTGACGCCCGccg                                           >  1:1303972/1‑148 (MQ=255)
                                                                                                           tGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAAGCTGACGCCCGCAg                                           >  1:1215534/1‑150 (MQ=255)
                                                                                                             gCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCt                                         >  2:1353788/1‑150 (MQ=255)
                                                                                                                 aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                     <  2:2778051/150‑1 (MQ=255)
                                                                                                                       acCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAACAACGCTGACGCCCGCAGCTTGTATGaaaa                               >  2:4387787/1‑150 (MQ=255)
                                                                                                                        cccttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAACGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                              <  2:3049138/147‑1 (MQ=255)
                                                                                                                         cATTGGGAATTCCGCAGACACAGTCTGCAGAATTCCCATCGGCAAAGATTTTGCAGACGCTGTCTGCCAAATCTTTCCCCATGCCCCGGTCCACCAACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCGTTTTTGAAAAag                             >  1:661640/1‑150 (MQ=255)
                                                                                                                          aTTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                            <  1:19927/150‑1 (MQ=255)
                                                                                                                           ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAagag                           >  2:4111002/1‑150 (MQ=255)
                                                                                                                           ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTGTATGAAAAtgag                           >  2:3944202/1‑150 (MQ=255)
                                                                                                                                 aTTCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTc                     >  2:562635/1‑150 (MQ=255)
                                                                                                                                  ttCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGTTTTGGCAGACGCGGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCGACGTCCGGTTACTTTGCGTTAAAAAAGCAGACGCCCCCCGCGGTTTTGAAAAAGAGACGCTcc                    >  1:179764/1‑150 (MQ=255)
                                                                                                                                  ttCTGCAGACGGTGTCTGCGGAATTCCCATCGGCAAAGATTTTGCAGACGCTGTCTGCCAAATCATTCCCCGTGCCCCGGTCCACCAACCTCCGTTTACTTTACGTTAAAAACGCGGACGCCCGCAGCTTTTATGAAAAAGTGTCGCTcc                    >  1:45948/1‑150 (MQ=255)
                                                                                                                                      gCAGACACCGTCTGCAGAATTCCCATCGGCAATGATTTTGCAGAAGCGGTGTGCCAAAACATTCCCCGTGCCCAGGTCCACCAACCTCCGGTTACGTTCCGTTAAAAAAGCTGACCCCCCCCGCGTTTATGAAAAAGAGTCGcgccgcag                >  1:2587377/1‑148 (MQ=255)
                                                                                                                                       cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt               <  1:4708637/150‑1 (MQ=255)
                                                                                                                                         gACACCGTCTGCAGAATTCCCATCGGCAACGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTgg             <  1:2387664/150‑1 (MQ=255)
                                                                                                                                           cacCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGCCGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTGTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCt           <  1:2734667/150‑1 (MQ=255)
                                                                                                                                             ccGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTgg         <  1:4090447/150‑1 (MQ=255)
                                                                                                                                              cGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt        <  1:4536587/150‑1 (MQ=255)
                                                                                                                                              cGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt        <  2:3798965/150‑1 (MQ=255)
                                                                                                                                              cGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt        <  2:2502111/150‑1 (MQ=255)
                                                                                                                                               gTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTc       >  2:1127347/1‑150 (MQ=255)
                                                                                                                                                  tGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt    >  2:4630977/1‑150 (MQ=255)
                                                                                                                                                    cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc  >  2:2994194/1‑150 (MQ=255)
                                                                                                                                                    cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGGTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCGGTTc  >  2:417749/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGGATTTATGTCTAAGGTATAAGCGTGGGTTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTC  >  NC_000913/3368040‑3368337

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