Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 868,422 (C)8→10 coding (903/1872 nt) gsiA → glutathione ABC transporter ATP binding subunit GsiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913868,4141.C99.2% 494.2 / ‑3.0 125coding (895/1872 nt)gsiAglutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (55/69);  total (56/69)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
*NC_000913868,4142.C82.4% 422.2 / 87.1 125coding (895/1872 nt)gsiAglutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base . (13/9);  new base C (43/60);  total (56/69)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 184 total aligned reads displayed.

CCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑CCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTG  >  NC_000913/868314‑868506
                                                                                                     ||                                                                                            
 cGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑c                                                                                             >  1:1476848/1‑101 (MQ=255)
  gTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cc                                                                                            >  2:1476849/1‑101 (MQ=255)
   tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccc                                                                                           >  2:1476850/1‑101 (MQ=255)
   tGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccc                                                                                           <  1:1476616/101‑1 (MQ=255)
    gcgcTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACGTCCAGCGAAACAGG‑‑cccc                                                                                          <  1:1476464/101‑1 (MQ=255)
    gcgcTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccc                                                                                          >  1:1476853/1‑101 (MQ=255)
     cgcTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccc                                                                                         >  2:1476855/1‑101 (MQ=255)
     cgcTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccc                                                                                         >  1:1476854/1‑101 (MQ=255)
      gcTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccc                                                                                        >  1:1476858/1‑101 (MQ=255)
        tGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccc                                                                                      >  2:1476859/1‑101 (MQ=255)
        tGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccc                                                                                      <  2:1476433/101‑1 (MQ=255)
          ttAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccc                                                                                    <  1:1476542/101‑3 (MQ=255)
           tAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccca                                                                                   >  1:1476863/1‑98 (MQ=255)
           tAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccca                                                                                    >  2:1476862/1‑98 (MQ=255)
            aGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccccat                                                                                  >  1:1476867/1‑97 (MQ=255)
             gctgctGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccccatc                                                                                 <  2:1476487/101‑6 (MQ=255)
             gctgctGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccccatc                                                                                 <  1:1476427/101‑6 (MQ=255)
             gctgctGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccccatc                                                                                 <  2:1476575/101‑6 (MQ=255)
             gctgctGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccatcg                                                                                 <  1:1476516/101‑6 (MQ=255)
             gctgctGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccc                                                                                     >  1:1476869/1‑96 (MQ=255)
              ctgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGa                                                                                 <  1:1476480/101‑1 (MQ=255)
              ctgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccc                                                                                      >  2:1476870/1‑95 (MQ=255)
               tgctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccc                                                                                     >  1:1476872/1‑94 (MQ=255)
                gctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccacccccc                                                                                     >  1:1476873/1‑93 (MQ=255)
                 ctgTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCa                                                                              <  2:1476630/101‑1 (MQ=255)
                    ttCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGa                                                                            >  2:1476876/1‑101 (MQ=255)
                    ttCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccc                                                                                       >  2:1476877/1‑88 (MQ=255)
                       cGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGaaaa                                                                         <  2:1476502/101‑1 (MQ=255)
                       cGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccc                                                                                    >  2:1476881/1‑86 (MQ=255)
                        gCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGaaa                                                                          >  2:1476883/1‑99 (MQ=255)
                         cAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACg                                                                       <  2:1476497/101‑1 (MQ=255)
                         cAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACg                                                                       <  2:1476410/101‑1 (MQ=255)
                         cAACTTAGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACg                                                                       >  1:1476885/1‑101 (MQ=255)
                          aaCTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGc                                                                               >  1:1476887/1‑91 (MQ=255)
                             ttGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccc                                                                                    >  2:1476890/1‑80 (MQ=255)
                              tGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACggtggt                                                                  <  2:1476223/101‑1 (MQ=255)
                              tGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACggtggt                                                                  <  2:1476588/101‑1 (MQ=255)
                              tGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACggtggt                                                                  <  2:1476592/101‑1 (MQ=255)
                                gTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGa                                                               <  2:1476546/101‑1 (MQ=255)
                                  gCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccc                                                                                    <  2:1476505/77‑3 (MQ=255)
                                   cGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATgg                                                            >  2:1476895/1‑101 (MQ=255)
                                   cGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccc                                                                                    >  2:1476897/1‑74 (MQ=255)
                                    gATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATgg                                                            <  2:1476714/101‑1 (MQ=255)
                                      tGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCg                                                          >  2:1476902/1‑101 (MQ=255)
                                        aaaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGAt                                                              <  2:1476334/95‑1 (MQ=255)
                                        aaaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGaa                                                        >  1:1476905/1‑101 (MQ=255)
                                        aaaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAc                                                       >  1:1476903/1‑101 (MQ=255)
                                         aaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGaaaa                                                                         >  1:1476908/1‑83 (MQ=255)
                                         aaGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAc                                                       <  1:1476448/101‑1 (MQ=255)
                                          aGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGa                                                         >  1:1476909/1‑98 (MQ=255)
                                           gggTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCt                                                     <  1:1476690/101‑1 (MQ=255)
                                           gggTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccccc                                                                                   >  2:1476910/1‑66 (MQ=255)
                                             gTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGt                                                   <  1:1476699/101‑1 (MQ=255)
                                              ttAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtt                                                  >  1:1476912/1‑101 (MQ=255)
                                                aGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACggtgg                                                                   >  1:1476914/1‑82 (MQ=255)
                                                  aTTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTg                                                    <  1:1476353/95‑1 (MQ=255)
                                                  aTTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGa                                                         >  2:1476919/1‑89 (MQ=255)
                                                  aTTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccc                                                                                    >  1:1476920/1‑59 (MQ=255)
                                                   ttATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACg                                             <  1:1476622/101‑1 (MQ=255)
                                                    tATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGa                                            <  1:1476561/101‑1 (MQ=255)
                                                     aTCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCt                                                     <  1:1476507/91‑1 (MQ=255)
                                                     aTCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAg                                           >  1:1476923/1‑101 (MQ=255)
                                                       ccccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACg                                             >  2:1476924/1‑96 (MQ=255)
                                                         ccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGaaa                                                                          >  1:1476931/1‑66 (MQ=255)
                                                         ccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGt                                      >  2:1476929/1‑101 (MQ=255)
                                                         ccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccc                                                                                    >  2:1476934/1‑52 (MQ=255)
                                                         ccGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccc                                                                                      >  1:1476932/1‑52 (MQ=255)
                                                           gACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTa                                     >  2:1476936/1‑101 (MQ=255)
                                                             cGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAAc                                                       <  2:1476668/81‑1 (MQ=255)
                                                              gTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑cccccccccc                                                                                    >  1:1476939/1‑47 (MQ=255)
                                                                ttCCCGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCtt                                <  1:1476559/101‑1 (MQ=255)
                                                                ttCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccc                                                                                     >  2:1476941/1‑45 (MQ=255)
                                                                  cccGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGc                                        >  2:1476945/1‑91 (MQ=255)
                                                                   ccGTTGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACggtg                                                                    >  1:1476948/1‑62 (MQ=255)
                                                                    cGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccc                                                                                     >  2:1476949/1‑41 (MQ=255)
                                                                       tGATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCAccc                         <  1:1476569/101‑1 (MQ=255)
                                                                        gATATCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTAc                                              >  1:1476952/1‑79 (MQ=255)
                                                                         atatCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGttt                                                 <  1:1476512/75‑1 (MQ=255)
                                                                         atatCGCTTGAACATCCAGCGAAACAGG‑‑ccccccccc                                                                                     >  1:1476956/1‑36 (MQ=255)
                                                                          tatCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGtt                      <  1:1476656/101‑1 (MQ=255)
                                                                           atCGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCg                                                          <  2:1476820/64‑1 (MQ=255)
                                                                            tcccTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTa                                     <  2:1476700/81‑1 (MQ=255)
                                                                             cGCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtt                                                  <  1:1476339/70‑1 (MQ=255)
                                                                              gCTTGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTccc                  >  2:1476959/1‑101 (MQ=255)
                                                                               cTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCtt                >  1:1476960/1‑101 (MQ=255)
                                                                               cTTGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCtt                <  1:1476536/101‑1 (MQ=255)
                                                                                ttGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCtt                >  1:1476963/1‑101 (MQ=255)
                                                                                ttGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGttt                                                 <  2:1476550/67‑1 (MQ=255)
                                                                                ttGAACATCCAGCGAAACAGGCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCttt               >  1:1476962/1‑101 (MQ=255)
                                                                                 tGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCttt               <  2:1476610/101‑1 (MQ=255)
                                                                                 tGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCttt               <  1:1476605/101‑1 (MQ=255)
                                                                                 tGAACATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCttt               <  1:1476600/101‑1 (MQ=255)
                                                                                    aCATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGtttt                                                <  1:1476631/65‑1 (MQ=255)
                                                                                    aCATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgcg            >  1:1476967/1‑101 (MQ=255)
                                                                                     cATCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTgcgc           <  1:1476727/101‑1 (MQ=255)
                                                                                      aTCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCa          <  2:1476688/101‑1 (MQ=255)
                                                                                      aTCCAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCa          <  2:1476639/101‑1 (MQ=255)
                                                                                         cAGCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGtggtgg                                                                >  1:1476969/1‑44 (MQ=37)
                                                                                           gCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCg                                                          >  2:1476970/1‑48 (MQ=37)
                                                                                           gCGAAACAGGCCCCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGt     <  1:1476549/101‑1 (MQ=255)
                                                                                                     ||                                                                                            
CCGTGCGCTGTTAGCTGCTGTTCCGCAACTTGGTGCGATGAAAGGGTTAGATTATCCCCGACGTTTCCCGTTGATATCGCTTGAACATCCAGCGAAACAGG‑‑CCCCCCCCATCGAGCAGAAAACGGTGGTGGATGGCGAACCTGTTTTACGAGTGCGTAATCTTGTCACCCGTTTCCCTTTGCGCAGCGGTTTG  >  NC_000913/868314‑868506

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