Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,980,954 A→G L128L (TTA→TTG aslB → putative arylsulfatase regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,980,9540AG99.4% 1109.3 / NA 359L128L (TTA→TTGaslBputative arylsulfatase regulator
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base G (177/179);  minor base T (0/2);  total (177/182)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 377 total aligned reads displayed.

GAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAC  >  NC_000913/3980864‑3981043
                                                                                          |                                                                                         
gAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTc                                                                                         <  1:6140493/93‑1 (MQ=255)
 aaGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCg                                                                                        >  1:1859890/1‑93 (MQ=255)
   gATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCt                                                                                      >  2:2833315/1‑93 (MQ=255)
    aTAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTg                                                                                     <  2:7094537/93‑1 (MQ=255)
     tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTTGCTgg                                                                                    >  2:5606706/1‑93 (MQ=255)
     tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                                                                    <  2:4155238/93‑1 (MQ=255)
     tAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTgg                                                                                    >  1:3665964/1‑93 (MQ=255)
      aaGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGa                                                                                   <  2:743646/93‑1 (MQ=255)
        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                 >  2:787270/1‑93 (MQ=255)
        gTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATg                                                                                 <  1:6480019/93‑1 (MQ=255)
         tAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATgg                                                                                >  2:615808/1‑93 (MQ=255)
           aCAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGcc                                                                              <  1:3440797/93‑1 (MQ=255)
             aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTTTTGGTTTGTCGCTGGATGGCCCg                                                                            >  1:7625627/1‑93 (MQ=255)
             aGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCg                                                                            >  2:595605/1‑93 (MQ=255)
              gCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGc                                                                           >  2:1409622/1‑93 (MQ=255)
               cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                                                          <  1:2482155/93‑1 (MQ=255)
               cTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGcc                                                                          >  2:2003117/1‑93 (MQ=255)
                ttCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCt                                                                         >  2:6413410/1‑93 (MQ=255)
                 tCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTg                                                                        <  1:328110/93‑1 (MQ=255)
                  ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                                                       <  2:2509305/93‑1 (MQ=255)
                  ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                                                       <  1:4311925/93‑1 (MQ=255)
                  ccAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                                                       <  1:573413/93‑1 (MQ=255)
                  acAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTga                                                                       <  2:1877976/92‑1 (MQ=255)
                   cAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgag                                                                      >  1:4230882/1‑93 (MQ=255)
                    aGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTgaga                                                                     >  1:1003109/1‑93 (MQ=255)
                     gACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGAt                                                                    >  1:263340/1‑93 (MQ=255)
                       cTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTTGCGGAGCATCATTTTTTTTTTGGTTTTTTGCTGGGTGGTGCGCCCGAGATcc                                                                  >  2:7662040/1‑93 (MQ=255)
                       cTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGGGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGTCAGCCTGAGATcc                                                                  >  2:5972399/1‑93 (MQ=255)
                       cTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATcc                                                                  >  1:6380423/1‑93 (MQ=255)
                       cTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATcc                                                                  <  2:1521733/93‑1 (MQ=255)
                       cTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATcc                                                                  <  2:1135359/93‑1 (MQ=255)
                        tAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCca                                                                 >  1:2403895/1‑93 (MQ=255)
                         aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcac                                                                >  2:3767280/1‑93 (MQ=255)
                         aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcac                                                                >  2:6791019/1‑93 (MQ=255)
                         aaCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATCGCCCCCCTGAGATTcac                                                                >  2:4682476/1‑93 (MQ=255)
                          aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcaca                                                               <  1:5817460/93‑1 (MQ=255)
                          aCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCcaca                                                               <  2:7353080/93‑1 (MQ=255)
                            ggCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaat                                                             <  1:7009829/93‑1 (MQ=255)
                             gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                                                            >  2:4043495/1‑93 (MQ=255)
                             gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                                                            <  2:5337658/93‑1 (MQ=255)
                             gCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatc                                                            <  2:7934256/93‑1 (MQ=255)
                               gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                                                          <  1:3963960/93‑1 (MQ=255)
                               gTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaa                                                          <  1:5868783/93‑1 (MQ=255)
                                tgctgcTGGATGACGAATGGTGCGCGTCTCTCGCGGCGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCAcaatcaat                                                         <  1:4069930/93‑1 (MQ=255)
                                  ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                                                       <  1:2108612/93‑1 (MQ=255)
                                  ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                                                       <  2:5463791/93‑1 (MQ=255)
                                  ctgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAatat                                                       <  1:5456574/93‑1 (MQ=255)
                                   tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      >  2:3551756/1‑93 (MQ=255)
                                   tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      <  2:3891322/93‑1 (MQ=255)
                                   tgctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATc                                                      <  1:4542589/93‑1 (MQ=255)
                                     ctgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgc                                                    <  1:3755528/93‑1 (MQ=255)
                                      tgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgcg                                                   >  1:5194073/1‑93 (MQ=255)
                                      tgGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATcgcg                                                   >  2:2574221/1‑93 (MQ=255)
                                          tGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGGTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGAc                                               >  2:5044991/1‑93 (MQ=255)
                                              gAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTaaa                                           >  1:2005567/1‑93 (MQ=255)
                                                aTGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAgg                                         >  2:3439910/1‑93 (MQ=255)
                                                    tgcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGc                                     >  1:5948264/5‑93 (MQ=255)
                                                     gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCa                                    <  2:2700771/90‑1 (MQ=255)
                                                     gcgcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCa                                    <  1:1198702/90‑1 (MQ=255)
                                                       gcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGa                                  >  1:6684932/2‑93 (MQ=255)
                                                       gcGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGa                                  <  2:1698669/92‑1 (MQ=255)
                                                        cGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAc                                 <  2:4463149/93‑1 (MQ=255)
                                                         gTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAcc                                <  2:1261607/93‑1 (MQ=255)
                                                         gTTTCTCGCGGAGCATCATTTTCTTGATGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGAcc                                >  2:5240316/1‑93 (MQ=255)
                                                           ttCTCGCGGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCa                              <  1:5050069/93‑1 (MQ=255)
                                                               cgcgGAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCATGACTAAAGGTGGCAGGCCCACGCa                          <  1:4462439/93‑1 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATat                       <  2:4792285/93‑2 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                       >  1:3191755/1‑93 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                       <  2:6226696/93‑1 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                       >  1:564489/1‑93 (MQ=255)
                                                                  ggAGCATCATTTTCTTGTTGGTTTGTCGCTGGACGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATaa                       >  1:7184469/1‑93 (MQ=255)
                                                                   gAGCATCATGTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAg                      <  1:2253625/93‑1 (MQ=255)
                                                                     gCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCt                    <  2:6660644/93‑1 (MQ=255)
                                                                     gCATCATTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCt                    <  2:3355052/93‑1 (MQ=255)
                                                                       atcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTgg                  >  1:681565/1‑93 (MQ=255)
                                                                        tcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGt                 >  2:4937574/1‑93 (MQ=255)
                                                                        tcatTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGt                 <  1:3532798/93‑1 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               <  2:109560/93‑1 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               <  1:306563/93‑1 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               <  1:7489531/93‑1 (MQ=255)
                                                                          atTTTCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGAGCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGa               <  1:5528269/93‑1 (MQ=255)
                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt              <  2:7115999/93‑1 (MQ=255)
                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt              >  2:3864366/1‑93 (MQ=255)
                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt              >  1:4783930/1‑93 (MQ=255)
                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt              <  1:914829/93‑1 (MQ=255)
                                                                           ttttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGAt              >  1:7264280/1‑93 (MQ=255)
                                                                            tttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATg             <  1:1070398/93‑1 (MQ=255)
                                                                             ttCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGc            <  1:5234166/93‑1 (MQ=255)
                                                                              tCTTGTTGGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCg           <  1:1301884/93‑1 (MQ=255)
                                                                                 tgttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGc        >  2:7278092/1‑93 (MQ=255)
                                                                                  gttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGcc       <  2:3200982/93‑1 (MQ=255)
                                                                                   ttgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGccc      >  1:2814424/1‑93 (MQ=255)
                                                                                    tgGTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCt     <  2:2977110/93‑1 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGCCCCACGCATAAGCTGGTGATGCGTGCCCTGa   <  1:4584257/93‑1 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa   >  2:1363535/1‑93 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa   <  2:2352670/93‑1 (MQ=255)
                                                                                      gTTTGTCGCTGGATGGCCCGCCTGAGATCCACAAGCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGa   <  2:6943369/93‑1 (MQ=255)
                                                                                       tttGTCGCTGGATGGCCCGCCTGGGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAc  >  2:5863673/1‑93 (MQ=255)
                                                                                       tttGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAc  >  1:7935011/1‑93 (MQ=255)
                                                                                       tttGTCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAc  >  2:4369717/1‑93 (MQ=255)
                                                                                          |                                                                                         
GAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGGCCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGCCTGAGATCCACAATCAATATCGCGTGACTAAAGGTGGCAGACCCACGCATAAGCTGGTGATGCGTGCCCTGAC  >  NC_000913/3980864‑3981043

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