Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,897,431 T→C H327H (CAT→CAC b1815 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,897,4310TC95.8% 558.7 / 9.9 167H327H (CAT→CACb1815hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base C (71/89);  total (76/91)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 185 total aligned reads displayed.

TATTGGTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCATTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGTTGGCAT  >  NC_000913/1897341‑1897520
                                                                                          |                                                                                         
tATTGGTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCActt                                                                                         <  1:1307130/93‑1 (MQ=255)
tATTGGTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCActt                                                                                         <  1:4038993/93‑1 (MQ=255)
tATTGGTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCActt                                                                                         >  1:3354707/1‑93 (MQ=255)
tATTGGTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCActt                                                                                         >  2:3796416/1‑93 (MQ=255)
tATTGGTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGCATGTGTTTATTCCTATCGCGGAAGAACATCActt                                                                                         <  1:2511701/93‑1 (MQ=255)
  ttGGTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTaa                                                                                       >  1:2935102/1‑93 (MQ=255)
  ttGGTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTaa                                                                                       >  1:36889/1‑93 (MQ=255)
    ggTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGTATGTGTTTATTCCTATCGCGGAAGAACATCACTTAAtt                                                                                     >  2:1414977/1‑93 (MQ=255)
    ggTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAAtt                                                                                     >  2:25712/1‑93 (MQ=255)
      tgtAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATtgt                                                                                   <  2:914722/93‑1 (MQ=255)
       gtAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCATTTAATtgtg                                                                                  >  2:1846890/1‑93 (MQ=255)
       gtAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATtgtg                                                                                  >  1:2263445/1‑93 (MQ=255)
         agagATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGcc                                                                                <  1:3046871/93‑1 (MQ=255)
            gaTTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACt                                                                             <  1:5272877/93‑1 (MQ=255)
            gaTTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACt                                                                             >  1:2784345/1‑93 (MQ=255)
             aTTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTg                                                                            <  1:6023529/93‑1 (MQ=255)
              ttttGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGa                                                                           >  2:4990733/1‑93 (MQ=255)
              ttttGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGa                                                                           >  2:584708/1‑93 (MQ=255)
              ttttGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGa                                                                           >  1:4586237/1‑93 (MQ=255)
              ttttGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTCAATTGTGCCACTGa                                                                           <  1:2208498/93‑1 (MQ=255)
               tttGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGAc                                                                          >  1:5379221/1‑93 (MQ=255)
                 tgctgcGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCATTTAATTGTGCCACTGAccc                                                                        >  1:1148380/1‑93 (MQ=255)
                  gctgcGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCg                                                                       >  1:5173175/1‑93 (MQ=255)
                   ctgcGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGc                                                                      <  2:96379/93‑1 (MQ=255)
                    tgcGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCt                                                                     <  2:1375203/93‑1 (MQ=255)
                     gcgcTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTa                                                                    <  2:596779/93‑1 (MQ=255)
                     gcgcTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTa                                                                    <  2:667625/93‑1 (MQ=255)
                      cgcTGGAACAATCCGCGTCAGGCCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTAt                                                                   <  1:3025130/93‑1 (MQ=255)
                        cTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTAtgt                                                                 >  1:2120508/1‑93 (MQ=255)
                        cTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTAtgt                                                                 <  2:5722476/93‑1 (MQ=255)
                         tGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTAtgtg                                                                <  1:2898557/93‑1 (MQ=255)
                         tGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTAtgtg                                                                <  2:451966/93‑1 (MQ=255)
                          ggAACAATCCGCGTCAGGTCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGTCCCGCTATGTGa                                                               <  2:2395072/93‑1 (MQ=255)
                           gAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGAt                                                              <  2:1484065/93‑1 (MQ=255)
                           gAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGAt                                                              <  2:1074488/93‑1 (MQ=255)
                           gAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGAt                                                              <  1:1905494/93‑1 (MQ=255)
                           gAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGAt                                                              >  1:4905927/1‑93 (MQ=255)
                                aTCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAg                                                         >  1:1029716/1‑93 (MQ=255)
                                 accGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGa                                                        <  2:3939410/92‑1 (MQ=255)
                                   cgcgTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGaaa                                                      <  1:3868727/93‑1 (MQ=255)
                                    gcgTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAAc                                                     <  2:2595925/93‑1 (MQ=255)
                                    gcgTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAAc                                                     <  1:5918383/93‑1 (MQ=255)
                                      gTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCa                                                   >  2:2355408/1‑93 (MQ=255)
                                      gTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCa                                                   <  2:1622625/93‑1 (MQ=255)
                                       tCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCATTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCAt                                                  >  1:4552698/1‑93 (MQ=255)
                                       tCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCAt                                                  >  1:1019781/1‑93 (MQ=255)
                                       tCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCAt                                                  >  1:6123512/1‑93 (MQ=255)
                                        cAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCAtt                                                 <  2:1236899/93‑1 (MQ=255)
                                         aGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTc                                                <  1:1484851/93‑1 (MQ=255)
                                          gggCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCg                                               >  2:2012497/1‑93 (MQ=255)
                                              tGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAg                                           <  1:1921867/93‑1 (MQ=255)
                                               ggATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAgc                                          <  1:4769291/93‑1 (MQ=255)
                                               ggATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAgc                                          >  2:3336330/1‑93 (MQ=255)
                                               ggATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAgc                                          >  1:3299292/1‑93 (MQ=255)
                                                gATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCATTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAgcg                                         <  2:2697261/93‑1 (MQ=255)
                                                   ttCACCGAAAGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCa                                      <  1:875192/93‑1 (MQ=255)
                                                    tCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCAt                                     <  1:5143925/93‑1 (MQ=255)
                                                      aCCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGt                                   <  1:3947038/93‑1 (MQ=255)
                                                      aCCGGATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGt                                   >  2:2056400/1‑93 (MQ=255)
                                                         ggATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGtttt                                <  1:4698968/93‑1 (MQ=255)
                                                          gATGTGTTTATTCCTATCGCGGAAGAACATCATTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTc                               >  2:5129086/1‑93 (MQ=255)
                                                          gATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTc                               <  1:3340110/93‑1 (MQ=255)
                                                          gATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTc                               <  2:2722622/93‑1 (MQ=255)
                                                          gATGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTc                               >  1:3789775/1‑93 (MQ=255)
                                                           aTGTGTTTATTCCTATCGCGGAAGAACATCATTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTcc                              <  1:4172123/93‑1 (MQ=255)
                                                           aTGTGTTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTcc                              <  2:3267755/93‑1 (MQ=255)
                                                            tgtgTTTATTCCTATCGCGGACGAACATAATTTAATTGTGCAACTGACCCGCTATTCGATGTCAGAACCAATTCGCAAGCGAAATGCTTcccc                             >  2:218707/1‑93 (MQ=255)
                                                             gtgtTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCg                            >  2:3371541/1‑93 (MQ=255)
                                                             gtgtTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCg                            >  1:801512/1‑93 (MQ=255)
                                                              tgtTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGa                           <  1:3619104/93‑1 (MQ=255)
                                                               gtTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAAGCATTCGTCAGCGCCATGTTTTCCCGAt                          >  1:3448913/1‑93 (MQ=255)
                                                               gtTTATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGAt                          <  1:1497432/93‑1 (MQ=255)
                                                                 ttATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGa                        <  2:560442/93‑1 (MQ=255)
                                                                 ttATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGa                        >  1:1650664/1‑93 (MQ=255)
                                                                 ttATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGa                        <  1:4927222/93‑1 (MQ=255)
                                                                  tATTCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGag                       <  1:1504347/93‑1 (MQ=255)
                                                                    ttCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGagta                     <  2:2108097/93‑1 (MQ=255)
                                                                    ttCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGagta                     <  1:4907323/93‑1 (MQ=255)
                                                                    ttCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGagta                     <  1:3344573/93‑1 (MQ=255)
                                                                    ttCCTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGagta                     >  1:2336440/1‑93 (MQ=255)
                                                                       cTATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTc                  >  1:5196518/1‑93 (MQ=255)
                                                                        tATCGCGGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGCGTAGTCa                 >  2:5376634/1‑93 (MQ=255)
                                                                         agcgcaGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAg                <  2:164205/91‑1 (MQ=255)
                                                                           cgctGAAGAACATCACTTGACTGGGCCACGGACCCGCGATCTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGtt              <  2:4428056/93‑1 (MQ=255)
                                                                           cgcgGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGtt              <  1:1408627/93‑1 (MQ=255)
                                                                             cgGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTc            >  2:4847140/1‑93 (MQ=255)
                                                                             cgGAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTc            >  2:2794662/1‑93 (MQ=255)
                                                                              ggAAGAACATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCa           >  1:2824693/1‑93 (MQ=255)
                                                                               gaagaaCATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCAt          >  2:2195471/1‑93 (MQ=255)
                                                                                aagaaCATCTCTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATg         <  2:5195238/93‑1 (MQ=255)
                                                                                aagaaCATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATg         >  1:2628538/1‑93 (MQ=255)
                                                                                  gaaCATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGtt       >  2:3907855/1‑93 (MQ=255)
                                                                                  gaaCATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGtt       <  2:3506816/93‑1 (MQ=255)
                                                                                   aaCATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGTTg      >  2:3348878/1‑93 (MQ=255)
                                                                                    aCATCACTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGTTgg     <  2:5707896/93‑1 (MQ=255)
                                                                                     catcacTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGTTGGc    <  1:736699/93‑1 (MQ=255)
                                                                                     catcacTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGTTGGc    >  2:3129483/1‑93 (MQ=255)
                                                                                     catcacTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGTTGGc    >  1:1276381/1‑93 (MQ=255)
                                                                                      atcacTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGTTGGCa   >  2:3759744/1‑93 (MQ=255)
                                                                                       tcacTTCATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGAAGCAAGTTTCATGTTGGCAt  >  1:932672/1‑93 (MQ=255)
                                                                                          |                                                                                         
TATTGGTGTAGAGATTTTGCTGCGCTGGAACAATCCGCGTCAGGGCTGGATTTCACCGGATGTGTTTATTCCTATCGCGGAAGAACATCATTTAATTGTGCCACTGACCCGCTATGTGATGGCAGAAACCATTCGTCAGCGCCATGTTTTCCCGATGAGTAGTCAGTTTCATGTTGGCAT  >  NC_000913/1897341‑1897520

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: