Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,526,598 G→A 5.7% A118T (GCT→ACT)  yncG → putative glutathione S‑transferase YncG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,526,5980GA 5.7% 222.9 / 2.6 70A118T (GCT→ACT) yncGputative glutathione S‑transferase YncG
Reads supporting (aligned to +/- strand):  ref base G (30/36);  new base A (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGCGATTGCATTGATGGTTCTCGATCGTCGCCCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCTTTGCACTATGC  >  NC_000913/1526450‑1526742
                                                                                                                                                    |                                                                                                                                                
ggcgATTGCATTGATGGTTCTCGATCGTCGCCCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGc                                                                                                                                                 >  1:1270604/1‑150 (MQ=255)
   gATTGCATTGATGGTTCTCGATCGTCGCCCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTcc                                                                                                                                              <  2:1088918/150‑1 (MQ=255)
   gATTGCATTGATGGTTCTCGATCGTCGCCCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTcc                                                                                                                                              <  1:921000/150‑1 (MQ=255)
          ttGATGGTTCTCGATCGTCGCCCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCAATTGGTCAGGCTGGTTGCTAAAGTCTATCCAACATTGACATTCGCCGATTACCCAGAGCGTTGGGCTCCTGACGCa                                                                                                                                       >  2:1081185/1‑150 (MQ=255)
                  tctcGATCGTCGCGCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCAGTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACa                                                                                                                               <  2:336212/150‑1 (MQ=255)
                    tcGATCGTCGCCCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGt                                                                                                                             >  1:917390/1‑150 (MQ=255)
                      gATCGTCGCCCGGACCTTGTGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTa                                                                                                                           >  2:1170945/1‑150 (MQ=255)
                        tcgtcgCCCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTaaa                                                                                                                         <  2:444830/150‑1 (MQ=255)
                           tcgCCCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGTTTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAaagaa                                                                                                                      <  2:1295062/150‑1 (MQ=255)
                               ccGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAAt                                                                                                                  >  1:1283203/1‑150 (MQ=255)
                                cGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTAGCCTGAGCGTTGGGCTCCTGACGCACCTGTTCAGTTAAAGAAGAATg                                                                                                                 <  2:636469/150‑1 (MQ=255)
                                 ggACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGt                                                                                                                <  1:39968/150‑1 (MQ=255)
                                  gACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGtt                                                                                                               <  2:13907/150‑1 (MQ=255)
                                   aCCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTa                                                                                                              <  1:175318/150‑1 (MQ=255)
                                    ccTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTAt                                                                                                             >  1:655293/1‑150 (MQ=255)
                                      ttGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTg                                                                                                           <  2:1270604/150‑1 (MQ=255)
                                        gcgcCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGaa                                                                                                         <  2:584966/150‑1 (MQ=255)
                                            cACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTTGTCTGGCTGTTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATc                                                                                                     <  1:546627/150‑1 (MQ=255)
                                              ccGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGa                                                                                                   <  2:782962/150‑1 (MQ=255)
                                                ggTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGaaa                                                                                                 <  2:394347/150‑1 (MQ=255)
                                                    gggCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCg                                                                                             <  1:1039162/150‑1 (MQ=255)
                                                       cgcgcCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCtt                                                                                          <  1:1332967/150‑1 (MQ=255)
                                                       cgcgcCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGACTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCtt                                                                                          >  2:165743/1‑150 (MQ=255)
                                                       cgcgcCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGACTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCtt                                                                                          >  2:381271/1‑150 (MQ=255)
                                                       cgcgcCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTACGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCtt                                                                                          >  1:730804/1‑150 (MQ=255)
                                                        gcgcCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCttt                                                                                         <  1:1331595/150‑1 (MQ=255)
                                                                 gTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCAACTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGt                                                                                >  1:782161/1‑150 (MQ=255)
                                                                  tCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGtt                                                                               <  1:159352/150‑1 (MQ=255)
                                                                   cAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTg                                                                              <  1:879151/150‑1 (MQ=255)
                                                                       tGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAAtt                                                                          >  1:230454/1‑150 (MQ=255)
                                                                        gTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTc                                                                         <  1:686992/150‑1 (MQ=255)
                                                                           tCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCa                                                                      <  2:917390/150‑1 (MQ=255)
                                                                           tCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCa                                                                      <  2:913461/150‑1 (MQ=255)
                                                                           tCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCa                                                                      <  1:976686/150‑1 (MQ=255)
                                                                             aaCGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGc                                                                    >  1:658984/1‑150 (MQ=255)
                                                                              aCGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCt                                                                   >  1:102117/1‑150 (MQ=255)
                                                                               cGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTc                                                                  >  1:1186676/1‑150 (MQ=255)
                                                                                       ggTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGa                                                          >  1:679925/1‑150 (MQ=255)
                                                                                       ggTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGa                                                          >  1:1123442/1‑150 (MQ=255)
                                                                                          ctggctggTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAAcc                                                       >  2:152857/1‑150 (MQ=255)
                                                                                            ggctggTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTcc                                                                                                                                              <  2:272310/61‑1 (MQ=255)
                                                                                            ggctggTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTcc                                                                                                                                              >  1:272310/1‑61 (MQ=255)
                                                                                                 gTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCg                                                >  2:88931/1‑150 (MQ=255)
                                                                                                    gCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCa                                                                      <  2:582575/125‑1 (MQ=255)
                                                                                                    gCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCa                                                                      >  1:582575/1‑125 (MQ=255)
                                                                                                     cTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTg                                            >  1:1205692/1‑150 (MQ=255)
                                                                                                      tAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTgg                                           <  1:152857/150‑1 (MQ=255)
                                                                                                       aaTGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACATGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGt                                          >  2:456724/1‑150 (MQ=255)
                                                                                                        aTGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCGCACTGCTGAACCTTATGCGTTTGGTg                                         <  2:1348781/150‑1 (MQ=255)
                                                                                                            cTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCa                                     >  1:962953/1‑150 (MQ=255)
                                                                                                                ccAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTa                                 >  2:252180/1‑150 (MQ=255)
                                                                                                                ccAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTa                                 <  2:590102/150‑1 (MQ=255)
                                                                                                                 cAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTaa                                >  2:29136/1‑150 (MQ=255)
                                                                                                                   aCATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACg                              <  2:132543/150‑1 (MQ=255)
                                                                                                                     aTTCACTTTCTCCTATTACCATGAGCGTTGGGCTACTGACGCACCTGAACAGTTAAAGAAGAATTTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTTTGCGTTTGGTGAGCAACTAACGCt                            <  2:1266873/150‑1 (MQ=255)
                                                                                                                      ttCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCtg                           >  1:165940/1‑150 (MQ=255)
                                                                                                                       tCACTTTCGCCGATTACCCTGAGCGTTGGACTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCtgg                          <  1:165743/150‑1 (MQ=255)
                                                                                                                       tCACTTTCGCCGATTACCCTGAGCGTTGGACTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCtgg                          <  1:381271/150‑1 (MQ=255)
                                                                                                                         aCTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCtggtg                        <  1:1227141/150‑1 (MQ=255)
                                                                                                                                ccGATTACCCTGAGCGTTTGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGtt                 <  2:962953/150‑1 (MQ=255)
                                                                                                                                ccGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGtt                 <  2:1283203/150‑1 (MQ=255)
                                                                                                                                      aCCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCttt           >  1:805031/1‑150 (MQ=255)
                                                                                                                                      aCCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCttt           >  1:965666/1‑150 (MQ=255)
                                                                                                                                       cccTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCTTTg          >  1:2217/1‑150 (MQ=255)
                                                                                                                                       cccTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCTTTg          >  1:565400/1‑150 (MQ=255)
                                                                                                                                         cTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCTTTGCa        <  2:102117/150‑1 (MQ=255)
                                                                                                                                         cTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCTTTGCa        <  2:1205692/150‑1 (MQ=255)
                                                                                                                                          tGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCTTTGCAc       <  1:1170945/150‑1 (MQ=255)
                                                                                                                                             gCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCTTTGCACTAt    >  1:1045578/1‑150 (MQ=255)
                                                                                                                                               gTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCTTTGCACTATgc  >  2:535099/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGCGATTGCATTGATGGTTCTCGATCGTCGCCCGGACCTTGCGCCACCGGTTGGGCGCGCCGAACGTCAGCTGTTTCAACGGCTATTGGTCTGGCTGGTTGCTAATGTCTATCCAACATTCACTTTCGCCGATTACCCTGAGCGTTGGGCTCCTGACGCACCTGAACAGTTAAAGAAGAATGTTATTGAATATCGAAAATCGCTTTATATCTGGTTGAATTCGCAGCTCACTGCTGAACCTTATGCGTTTGGTGAGCAACTAACGCTGGTGGACTGTTATCTTTGCACTATGC  >  NC_000913/1526450‑1526742

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