Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,818,046 A→G 9.6% intergenic (‑186/+13) argQ ← / ← argZ tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,818,0460AG 9.6% 288.5 / 17.4 97intergenic (‑186/+13)argQ/argZtRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base A (54/32);  new base G (4/5);  total (59/38)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTA  >  NC_000913/2817897‑2818192
                                                                                                                                                     |                                                                                                                                                  
gCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGa                                                                                                                                                    >  2:2161212/1‑150 (MQ=255)
           cATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAat                                                                                                                                         <  2:901350/150‑1 (MQ=255)
              atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                      <  2:2161415/150‑1 (MQ=255)
              atCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTGGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTCGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATg                                                                                                                                      <  2:383397/150‑1 (MQ=255)
                     aaGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATc                                                                                                                               >  1:607660/1‑150 (MQ=255)
                     aaGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATc                                                                                                                               <  2:2135767/150‑1 (MQ=255)
                     aaGAGATGGTGCATCCGGGAGGATTACTCGCCGTGCGGCTCGCCCTACGGGCCGTTGCTATAGCAACGTTATCATCGCTATTGCGTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATc                                                                                                                               <  2:1616862/150‑1 (MQ=18)
                       gagaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCg                                                                                                                             >  1:1507880/1‑150 (MQ=255)
                       gagaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCg                                                                                                                             >  1:1507836/1‑150 (MQ=255)
                       gagaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCg                                                                                                                             >  1:1479392/1‑150 (MQ=255)
                       gagaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGATGCTAAAGCAACGATATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTAGATACCTTCATTGCTGAAGAATTCGGAGAATATTGTGCATCCg                                                                                                                             >  1:2522195/1‑150 (MQ=17)
                        agaTGGTGCATCCGGGAGGATTACTCGCGTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCTTACCTAGTGCCTGCGGTTAATTCCGCTAGCTCTTCGATACCTTCATTGCTGTAGTTTTCGGAGATTATGGTGCTTCCgg                                                                                                                            >  2:1249366/1‑150 (MQ=17)
                        agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  1:1034743/1‑150 (MQ=255)
                        agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  1:2314650/1‑150 (MQ=255)
                        agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  1:1877708/1‑150 (MQ=255)
                        agaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCgg                                                                                                                            >  2:2473761/1‑150 (MQ=255)
                         gaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCACCGGTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTGCATAGCTGAAGATTTCGGAGAATATGGTGCATCCggg                                                                                                                           >  1:1790773/1‑150 (MQ=255)
                         gaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCggg                                                                                                                           >  1:125828/1‑150 (MQ=255)
                         gaTGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCggg                                                                                                                           >  1:1922603/1‑150 (MQ=255)
                          aTGGTGCACCCGGGAGGCTAACTCGCCCTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGgga                                                                                                                          <  2:1441510/150‑1 (MQ=255)
                             gTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGggagga                                                                                                                       >  2:1560254/1‑150 (MQ=255)
                                 aTCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCg                                                                                                                   >  1:297344/1‑150 (MQ=255)
                                   ccGGGAGGTGTTCTCGCCCTGCGGCTCGCCCTTCGGTCCGTTGCTTTAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGaa                                                                                                                 <  2:560233/150‑1 (MQ=18)
                                   ccGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTCCAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGaa                                                                                                                 <  2:656769/150‑1 (MQ=255)
                                    cGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGACTATGGTGCATCCGGGAGGATTCGAAc                                                                                                                >  1:1851056/1‑150 (MQ=255)
                                     gggAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAAcc                                                                                                               >  1:1558479/1‑150 (MQ=255)
                                              aCTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGc                                                                                                      <  2:1690662/150‑1 (MQ=255)
                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                    >  2:776651/1‑150 (MQ=255)
                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTc                                                                                                    >  2:776745/1‑150 (MQ=255)
                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATCCGAACCTCCGACCGCTc                                                                                                    >  1:1939214/1‑150 (MQ=255)
                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATCCGAACCTCCGACCGCTc                                                                                                    >  1:1939247/1‑150 (MQ=255)
                                                tCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCAGCGGTTAATTCCGCTACCTCTTAGATACCTTCATTGCTGAAGATTTCGGAGAATATTGTGCATCCGGGAGGATCCGAACCTCCGACCGCTc                                                                                                    >  1:659252/1‑150 (MQ=255)
                                                     ttGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGGTc                                                                                               >  1:2530891/1‑150 (MQ=255)
                                                     ttGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATACCTAGTGCCTGCGGTTAATTACGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTc                                                                                               >  1:2531024/1‑150 (MQ=255)
                                                           ttCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGga                                                                                                                                             >  2:2091546/2‑98 (MQ=255)
                                                           ttCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGga                                                                                                                                             <  1:2091546/97‑1 (MQ=255)
                                                             cGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGa                                                                                       <  1:737596/150‑1 (MQ=255)
                                                                 cTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTAGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAc                                                                                   >  1:1018292/1‑150 (MQ=255)
                                                                     gggCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                               >  2:537518/1‑150 (MQ=255)
                                                                       gCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                             <  1:2021512/150‑1 (MQ=255)
                                                                           ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                                         >  2:1203503/1‑150 (MQ=255)
                                                                           ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                                         >  2:1203482/1‑150 (MQ=255)
                                                                           ttGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGc                                                                         >  2:1632641/1‑150 (MQ=255)
                                                                                 aaGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                                   >  1:2012312/1‑150 (MQ=255)
                                                                                 aaGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                                   >  1:138226/1‑150 (MQ=255)
                                                                                        gTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATg                                                            >  1:2167543/1‑150 (MQ=255)
                                                                                            tcatcCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATc                                                        <  1:409864/150‑1 (MQ=255)
                                                                                                     gTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAattt                                               >  1:1091984/1‑150 (MQ=255)
                                                                                                     gTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAagtt                                               >  1:2258606/1‑150 (MQ=255)
                                                                                                           gCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTa                                         >  1:781070/1‑150 (MQ=255)
                                                                                                            cGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTAc                                        >  2:221285/1‑150 (MQ=255)
                                                                                                                  aTTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa                                  <  2:2547131/150‑1 (MQ=255)
                                                                                                                     ccGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGatat                               <  2:1507880/150‑1 (MQ=255)
                                                                                                                     ccGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGatat                               <  2:1507836/150‑1 (MQ=255)
                                                                                                                     ccGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGatat                               <  2:1479392/150‑1 (MQ=255)
                                                                                                                      cGCTACCTCTTCGATACCTTCATTGCTGAAGGTTTCGGAGAATATGGTGCATCCGGGAGGATTCGa                                                                                                                  <  1:322617/66‑1 (MQ=25)
                                                                                                                       gCTAGCTCTTCGATACCATCATTTCTGAAGATTTAGGATAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGAttttt                             <  2:2531024/150‑2 (MQ=11)
                                                                                                                             tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGACACTTACCTTACTGCTGATATTGAGAcac                       >  2:1171853/1‑148 (MQ=11)
                                                                                                                             tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCg                       >  2:2395011/1‑150 (MQ=255)
                                                                                                                             tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCg                       >  2:346433/1‑150 (MQ=255)
                                                                                                                             tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACATCCGACCGCTCGGTTCGTAGCCGAGTACACAATACAGATGAGCTACGGATGCATCGGGAAATTTATGTTACTGCTGATATTGTTACCg                       >  1:2256093/1‑150 (MQ=14)
                                                                                                                             tctTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACAGCCGACCGCTCGGTTCTTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGAATCGGGAAAGTTATTTTACTGCTGATATTGATAGCg                       >  1:1613155/1‑150 (MQ=14)
                                                                                                                               ttCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGt                                                                                                                                    <  1:2020280/39‑1 (MQ=255)
                                                                                                                               ttCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAGTTATTTTACTGCTGATATTGATACCGct                     >  2:2159346/1‑150 (MQ=255)
                                                                                                                                tCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                    <  2:811619/150‑1 (MQ=255)
                                                                                                                                tCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                    <  2:579508/150‑1 (MQ=255)
                                                                                                                                tCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGcta                    <  2:171045/150‑1 (MQ=255)
                                                                                                                                 cGATACCTTCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGctac                   <  1:1973602/150‑1 (MQ=18)
                                                                                                                                    tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCtaataa                >  1:1038987/1‑150 (MQ=255)
                                                                                                                                    tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTaa                >  1:431941/1‑150 (MQ=255)
                                                                                                                                    tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACAGCTACTaa                >  1:1059327/1‑150 (MQ=255)
                                                                                                                                    tACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAGTTAATTTACTGCTGATATTGTTAACGCtaataa                >  1:1059343/1‑150 (MQ=14)
                                                                                                                                      cttcgtttgttgAAGTTTTCGGAGAATATGGTGCATCGGGGAGGATTCGAACCTCCGACCGCTGGGTTCGTAGCCGTGTACTCTATCCAGCTGAGCTACGGATGCATGGGTAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              <  1:1143430/144‑1 (MQ=11)
                                                                                                                                      ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              <  1:473396/150‑1 (MQ=255)
                                                                                                                                      ccTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              <  1:1143429/150‑1 (MQ=255)
                                                                                                                                        ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTCCGGATGCATCGGGAAATATATTTGACGGCTGATATTGATACCGCTACTAAAGcc            >  2:2526085/1‑150 (MQ=11)
                                                                                                                                        ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc            >  2:1025157/1‑150 (MQ=18)
                                                                                                                                        ttCATTGCTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc            >  2:191668/1‑150 (MQ=18)
                                                                                                                                        ttCATTGCTGAAGATTTCGGAGTATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCg                                                       <  2:1370478/107‑1 (MQ=25)
                                                                                                                                        ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCg                                                       >  1:1370478/1‑107 (MQ=35)
                                                                                                                                        ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc            >  1:1600896/1‑150 (MQ=255)
                                                                                                                                         tCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCa           >  1:1402424/1‑150 (MQ=255)
                                                                                                                                            ttGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCatat        <  1:1203482/150‑1 (MQ=255)
                                                                                                                                            ttGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCatat        <  1:1203503/150‑1 (MQ=255)
                                                                                                                                            ttGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCatat        <  1:1632641/150‑1 (MQ=255)
                                                                                                                                              ccTGAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCa      <  1:717246/149‑1 (MQ=18)
                                                                                                                                               cTGAAGATTTCTGTGTATTTGGTGCTTGGGGGATGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa     <  1:1463915/150‑1 (MQ=14)
                                                                                                                                               cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTTATAACGCTAATAAAGCCATATCaa     >  1:933721/1‑150 (MQ=255)
                                                                                                                                               cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa     >  2:138882/1‑150 (MQ=255)
                                                                                                                                               cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTAGAACCACAGACCCCACCCGTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCAGCGCGCAATTTATGTTCCGGCGGATAATGATACCGATACTCAAGCAATATcta     >  1:933781/1‑148 (MQ=9)
                                                                                                                                                tGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGa                                                                                                          <  1:2398504/48‑1 (MQ=255)
                                                                                                                                                tGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAg    <  2:1558479/150‑1 (MQ=255)
                                                                                                                                                 gAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              >  2:1143464/1‑139 (MQ=18)
                                                                                                                                                 gAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAg              <  1:1143464/139‑1 (MQ=18)
                                                                                                                                                 gAAGGTTACGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGt   <  1:2504094/150‑1 (MQ=25)
                                                                                                                                                  aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCt                                                                                                     <  2:1987118/51‑1 (MQ=255)
                                                                                                                                                  aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCt                                                                                                     >  1:1987118/1‑51 (MQ=255)
                                                                                                                                                  aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATAACGCTACAAAAGCCATATCAAGTa  >  2:1777102/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCTACTAAAGCCATATCAATTAAGAGATGGTGCATCCGGGAGGATTACTCGCCTTGCGGCTCGCCCTTCGGGCCGTTGCTAAAGCAACGTTATCATCCCTAGTGCCTGCGGTTAATTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTA  >  NC_000913/2817897‑2818192

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: