Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 892,459 G→C 7.4% D165H (GAC→CAC)  rimK → ribosomal protein S6 modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913892,4590GC 7.4% 334.2 / 16.6 108D165H (GAC→CAC) rimKribosomal protein S6 modification protein
Reads supporting (aligned to +/- strand):  ref base G (46/54);  new base C (3/5);  total (49/59)
Fisher's exact test for biased strand distribution p-value = 7.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 111 total aligned reads displayed.

CTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTG  >  NC_000913/892310‑892608
                                                                                                                                                      |                                                                                                                                                     
ctgTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTg                                                                                                                                                       >  2:2053085/1‑150 (MQ=255)
  gTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGAc                                                                                                                                                     <  1:1202041/150‑1 (MQ=255)
   gcACGGGCATTGCGCAATCGCCGGAGGATACCAGCGATTTACTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTATGATGAAGGCACGCAGGGAAATGGCGTCGTGCTGTCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACg                                                                                                                                                    <  1:1866640/149‑1 (MQ=255)
    ctcGGGCATTGCGCATTCGGCGGATGATACCAGCGATTTAATCGACCTGGTCGGTGGTTCGCCGCTGGTGGTCATGTTGGCTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGc                                                                                                                                                   <  2:133907/148‑1 (MQ=255)
    cACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGc                                                                                                                                                   >  2:308544/1‑150 (MQ=255)
     aCGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCt                                                                                                                                                  >  1:1273793/1‑150 (MQ=255)
        ggCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGTTCGGTGGTGCGCCGGTGGTGGTCAAGTTGGTTGAAGGCACTCAGGGAATTGGCGTCGTGCTGGCTGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTc                                                                                                                                               <  1:356395/150‑1 (MQ=255)
        ggCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTc                                                                                                                                               <  1:1109993/150‑1 (MQ=255)
        ggCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTc                                                                                                                                               >  1:429438/1‑150 (MQ=255)
           aTTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCcg                                                                                                                                             <  1:1703348/150‑1 (MQ=38)
           aTTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCTGGCGGCGGTAAGCGTGATTGACGCTTTCcgc                                                                                                                                            >  1:144778/1‑150 (MQ=255)
           aTTGCGCATTCGCAGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCcgc                                                                                                                                            <  2:876051/150‑1 (MQ=255)
                   ttCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGaa                                                                                                                                    >  1:386934/1‑150 (MQ=255)
                    tCGCCGGATGATACCAGCGATTTACTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAc                                                                                                                                   <  1:1553009/150‑1 (MQ=255)
                    tCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGACCGCGGAAAGCGTCAATCACGCTTTCCGCCGCCTGacg                                                                                                                                   <  1:1381547/150‑3 (MQ=38)
                     cGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAAcg                                                                                                                                  >  2:1663013/1‑150 (MQ=255)
                      gCCGGATGATACCAGCGATTAAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCCGCCGCCTGacgc                                                                                                                                  <  1:1035636/150‑4 (MQ=25)
                       ccGGCTGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCCGCCGCCTGacgcg                                                                                                                                 <  1:330927/150‑5 (MQ=25)
                       ccGGATGATACCAGCGATTTTATCGTCATGGTCGGTGGTGCGCCGCTGGTTGTCAAGTTGGTGGAAGGCACGCAGGGAATTGGCGGCGTGCTGGCGGAGACGCGTTCAGACCGCGGAAAGCGTCAATCACGCTTTCCGCCGCCTGacgcg                                                                                                                                 <  1:139266/150‑5 (MQ=21)
                               tACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCt                                                                                                                        <  2:1253882/150‑1 (MQ=255)
                                aCCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTg                                                                                                                       >  1:1699694/1‑150 (MQ=255)
                                aCCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTg                                                                                                                       >  2:1469146/1‑150 (MQ=255)
                                 ccAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTgg                                                                                                                      <  2:2042130/150‑1 (MQ=255)
                                 ccAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTgg                                                                                                                      >  1:1364723/1‑150 (MQ=255)
                                    ccGATTTAATCGTCATGGGCGGTGGGGCGCCGCTGGTGGGCAAGTTGGTTGAAGGCACGGAGGGAATTGGCGTCGTGCGGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTGTCCGCGGTCTGAACGCTCATATTCTGGTGc                                                                                                                   <  2:388700/149‑1 (MQ=255)
                                     cGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCa                                                                                                                  >  1:272457/1‑150 (MQ=255)
                                      gATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGGCGTGCGGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAg                                                                                                                 <  2:1823267/150‑1 (MQ=255)
                                      gATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGGTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAg                                                                                                                 <  1:536773/150‑1 (MQ=255)
                                      gATTTAATCGACATGGGCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAg                                                                                                                 <  2:447150/150‑1 (MQ=255)
                                           aaTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAata                                                                                                            >  2:1801305/1‑150 (MQ=255)
                                            aTCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAatat                                                                                                           >  1:764123/1‑150 (MQ=255)
                                            aTCGACATGGTCGGTGGGGCGCCGCTGGTGGTCAAGTTGGGTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGGGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAatat                                                                                                           <  2:1522852/150‑1 (MQ=255)
                                              cGACATGGTCGGTGGTGCGCCGCTGGTGGTCATGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACTCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAatatat                                                                                                         <  2:1441473/150‑1 (MQ=255)
                                              cGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAatatat                                                                                                         <  1:2082106/150‑1 (MQ=255)
                                               gataTGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATc                                                                                                        >  1:989000/4‑150 (MQ=255)
                                               gACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATc                                                                                                        >  1:1377565/1‑150 (MQ=255)
                                                aCATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCa                                                                                                       <  2:549492/150‑1 (MQ=255)
                                                 cATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaa                                                                                                      <  1:56493/150‑1 (MQ=255)
                                                 cATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCaa                                                                                                      >  1:520404/1‑150 (MQ=255)
                                                   tGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAag                                                                                                    >  2:311015/1‑150 (MQ=255)
                                                    ggTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAaga                                                                                                   <  1:255046/150‑1 (MQ=255)
                                                    ggTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCAGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAaga                                                                                                   >  1:739353/1‑150 (MQ=255)
                                                            gtgCGCCGCTGGTGGTCAAGATGTGTGAAGGCACGCAGGGAATTGGCGTCGTGATGGCGTAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAg                                                                                           <  2:1953995/150‑1 (MQ=255)
                                                                 ccgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGc                                                                                      >  2:941251/1‑150 (MQ=255)
                                                                  cgcTGGTGGTCAAGTTGGTTGAAGGCCCGCAGGGATTTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGAAGGAATATATCAAAGAGGCGCAAGGGTGCg                                                                                     <  1:830802/150‑1 (MQ=255)
                                                                  cgcTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGGCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCg                                                                                     <  2:1480399/150‑1 (MQ=255)
                                                                       gtggtCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATc                                                                                >  1:59990/1‑150 (MQ=255)
                                                                        tggtCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATcc                                                                               >  1:115356/1‑150 (MQ=255)
                                                                        gggtCAAGTTGGTTGAAGGCACGCTGGGAATTGGCGTCGTGCTGGCTGAGTCGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTATGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATcc                                                                               <  2:923270/149‑1 (MQ=255)
                                                                        gggtCAAGTTGGGTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGTCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATcc                                                                               <  1:575098/149‑1 (MQ=255)
                                                                          gtCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGc                                                                             <  1:1246164/150‑1 (MQ=255)
                                                                          gtCAAGTTGGTTGAAGGCACGCAGGGAATTGGAGTCGTGCTGGCGGAGACGCGTCAGCCGGCGGAAAGCGTGATTGACGCTTTCAGCGGTCAGAACGCGCATATTCTGGTGCAGGAATATATCCAAGAGGCGCAAGGGTGCGATATCCGc                                                                             >  1:137673/1‑150 (MQ=255)
                                                                           tCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCt                                                                            >  1:1902530/1‑150 (MQ=255)
                                                                            cAAGTTGGTTTAAGGCACGCAGGGAATTGGCGTCGTGCTGGCTGAGACGCGGCAGGCGGCTGTAAGCGTGATTGACGCTTTCCTCGGTCTGTACGCGCATATTCTGGTGCTGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTg                                                                           <  1:721064/150‑1 (MQ=255)
                                                                               gttggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCt                                                                        <  2:745621/150‑1 (MQ=255)
                                                                                 tggttgAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGCCGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATAAATCAATGAGGCGCAAGGGTGCGATATCCGCTGTCTgg                                                                      >  2:234330/1‑150 (MQ=255)
                                                                                  ggttgAAGGCACGCAGGGAATTGGCGGCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgt                                                                     <  1:1474019/150‑1 (MQ=255)
                                                                                  ggtggAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgt                                                                     <  1:554310/150‑1 (MQ=255)
                                                                                   gttgAAGGCACGCAGGGAATTGCCGTCGTGCTGGCGTAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgtt                                                                    >  1:183310/1‑150 (MQ=255)
                                                                                    gggAAGGCACGCAGGGCAGTGGCGTCGTGCTGGCGGACACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGGCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttg                                                                   <  2:1150064/148‑1 (MQ=255)
                                                                                      gAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGgttgtt                                                                 >  1:27066/1‑150 (MQ=255)
                                                                                        aGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTgg                                                               <  2:1387126/150‑1 (MQ=255)
                                                                                         ggCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGc                                                              >  1:1220622/1‑150 (MQ=255)
                                                                                         ggCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGc                                                              >  1:1113636/1‑150 (MQ=255)
                                                                                         ggCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGc                                                              <  2:1794801/150‑1 (MQ=255)
                                                                                          gCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCg                                                             <  1:1514450/150‑1 (MQ=255)
                                                                                             cGCAGGGAATTGGCGCCGGGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGGCTGAACGCGCATATTCTGGTGCTGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATg                                                          <  2:454205/150‑1 (MQ=255)
                                                                                                 gggAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGAGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGt                                                      >  2:903086/1‑150 (MQ=255)
                                                                                                    aaTTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGt                                                   >  1:1990108/1‑150 (MQ=255)
                                                                                                     agtGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGGGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTc                                                  <  2:2058215/148‑1 (MQ=255)
                                                                                                       tGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGc                                                >  2:1454568/1‑150 (MQ=255)
                                                                                                            tcgtGCTGGCGCACCCGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGGGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGa                                           <  2:209052/150‑1 (MQ=255)
                                                                                                             cgtGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGGGAGTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGAt                                          <  2:1576554/150‑1 (MQ=255)
                                                                                                             cgggCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGAt                                          <  1:1972235/147‑1 (MQ=255)
                                                                                                                gCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGa                                       >  2:770435/1‑150 (MQ=255)
                                                                                                                        agaCGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACggcggg                               <  2:520404/150‑1 (MQ=255)
                                                                                                                           cgcgTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATGTCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGa                            >  1:1092431/1‑150 (MQ=255)
                                                                                                                              gTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAag                         <  2:1699694/150‑1 (MQ=255)
                                                                                                                                   ggcggcGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAggg                     >  1:2002474/1‑150 (MQ=255)
                                                                                                                                   gaccgcggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAggg                     >  2:1492749/5‑150 (MQ=255)
                                                                                                                                   gaccgcggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAggg                     >  2:1499890/5‑150 (MQ=255)
                                                                                                                                    gcggcgGAAAGCGTGATTGACGCTTTCCGCGGTGTGAACGCGCATATTCTGGTGCAGGAATATAGCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGc                    <  1:188833/150‑1 (MQ=255)
                                                                                                                                    gcggcgGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGc                    <  2:386934/150‑1 (MQ=255)
                                                                                                                                      ggcggAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGa                  >  1:1996826/1‑150 (MQ=255)
                                                                                                                                       gcggAAAGCGTCAATCACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAt                 >  1:769874/1‑150 (MQ=255)
                                                                                                                                        cggAAAGCGTGATTGATGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGAtt                >  1:1033724/1‑150 (MQ=255)
                                                                                                                                            aaGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGCGGCGCAAGGGTGCGATATCCCCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGTCGGTCGAAAGTGGGCGCTTTTCg            >  2:357540/1‑150 (MQ=255)
                                                                                                                                              gCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGGGAAAGAGGGCGATTTTCGtt          >  2:175695/1‑150 (MQ=255)
                                                                                                                                              gCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGtt          <  2:817807/150‑1 (MQ=255)
                                                                                                                                              gCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGtt          >  2:937746/1‑150 (MQ=255)
                                                                                                                                               cGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTc         <  2:764123/150‑1 (MQ=255)
                                                                                                                                               cGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTc         <  2:1712452/150‑1 (MQ=255)
                                                                                                                                               cGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGGGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTc         <  1:1773211/150‑1 (MQ=255)
                                                                                                                                                 tGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCa       >  2:557487/1‑150 (MQ=255)
                                                                                                                                                 tGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGATGTCTGGTTGTTGGCGATGAAGTGGTAGCTGCGATTGAACGGCGGTCGAAAGAGGGCTATTTTCGTTCCa       >  2:225861/1‑150 (MQ=255)
                                                                                                                                                   aTTGCCGCGTTCCGCGGTCGGAACGCGCATATACTGGTGCAGGAATATATCACCGCGGCGCACGGGTGCCATATCCGCTGGCTGGTTGTGGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAAt     <  2:144778/150‑1 (MQ=255)
                                                                                                                                                    ttGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAAtt    >  1:1767364/1‑150 (MQ=255)
                                                                                                                                                     tGACGCTTTCCGCGGGCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAAttt   <  1:1820701/150‑1 (MQ=255)
                                                                                                                                                      gACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTg  <  2:535441/150‑1 (MQ=255)
                                                                                                                                                      gACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGGTGTTGGCGATGAAGTGGGCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTg  <  1:2102857/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTGTCACGGGCATTGCGCATTCGCCGGATGATACCAGCGATTTAATCGACATGGTCGGTGGTGCGCCGCTGGTGGTCAAGTTGGTTGAAGGCACGCAGGGAATTGGCGTCGTGCTGGCGGAGACGCGTCAGGCGGCGGAAAGCGTGATTGACGCTTTCCGCGGTCTGAACGCGCATATTCTGGTGCAGGAATATATCAAAGAGGCGCAAGGGTGCGATATCCGCTGTCTGGTTGTTGGCGATGAAGTGGTCGCTGCGATTGAACGGCGGGCGAAAGAGGGCGATTTTCGTTCCAATTTG  >  NC_000913/892310‑892608

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