Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,106,773 A→G V120V (GTA→GTG sbp → periplasmic sulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,106,7730AG99.4% 552.3 / ‑5.6 170V120V (GTA→GTGsbpperiplasmic sulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (83/86);  total (83/87)
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.99e-01

Only 100 of 173 total aligned reads displayed.

ATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTG  >  NC_000913/4106684‑4106861
                                                                                         |                                                                                        
 ttGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTa                                                                                      <  2:1271382/93‑1 (MQ=255)
 ttGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTa                                                                                      <  2:3181060/93‑1 (MQ=255)
  tGCGGAACGCGGTCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTaa                                                                                     <  1:1099963/93‑1 (MQ=255)
  tGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTaa                                                                                     <  2:2872624/93‑1 (MQ=255)
  tGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTaa                                                                                     <  1:3439820/93‑1 (MQ=255)
   gCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAg                                                                                    <  2:1757009/93‑1 (MQ=255)
    cGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAgg                                                                                   <  2:2681820/93‑1 (MQ=255)
    cGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAgg                                                                                   <  2:2784995/93‑1 (MQ=255)
      gAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGa                                                                                 >  2:2859615/1‑93 (MQ=255)
       aaCGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGaa                                                                                >  1:2130985/1‑93 (MQ=255)
          gcgggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATc                                                                             <  2:1845260/93‑1 (MQ=255)
           cgggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATcc                                                                            <  1:1005149/93‑1 (MQ=255)
           cgggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATcc                                                                            >  1:3773614/1‑93 (MQ=255)
            gggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCg                                                                           >  2:367626/1‑93 (MQ=255)
            gggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCg                                                                           <  2:3458251/93‑1 (MQ=255)
             ggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGa                                                                          >  2:3988774/1‑93 (MQ=255)
             ggcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGa                                                                          >  1:579805/1‑93 (MQ=255)
              gcggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGaa                                                                         <  1:3230173/93‑1 (MQ=255)
                ggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACCCTTCCCCCATTGTTTTCCTGGGGCGTAAGGGAAATCCGCAGc                                                                       >  1:1424718/1‑93 (MQ=255)
                ggATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGGTTTCCTGGTGCGTAAGGGAAATCCGAAGc                                                                       >  1:4219585/1‑93 (MQ=255)
                  aTTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAg                                                                     <  2:1408395/93‑1 (MQ=255)
                  aTTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAg                                                                     >  1:3774926/1‑93 (MQ=255)
                  aTTGATAAAGAGTGGATCAAAAGTCTGCCGGATAACTCCGCACCGTACACTTACACCATTGTTTTCCTGGTGCGTAAGGGCAATCCGAAGCAg                                                                     >  2:2325002/1‑93 (MQ=255)
                   ttGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGa                                                                    >  2:2422976/1‑93 (MQ=255)
                   ttGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGa                                                                    >  2:2756018/1‑93 (MQ=255)
                   ttGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGa                                                                    >  1:4127607/1‑93 (MQ=255)
                    tGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGAt                                                                   <  2:3421615/93‑1 (MQ=255)
                    tGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGAt                                                                   <  1:1600161/93‑1 (MQ=255)
                      aTAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATcc                                                                 >  2:3573240/1‑93 (MQ=255)
                       tAATGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCa                                                                <  2:3955823/93‑1 (MQ=255)
                       tAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCa                                                                <  1:110704/93‑1 (MQ=255)
                         aaGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTATGGGAAACCCCGAGCAGATCCATg                                                              >  2:3386828/1‑93 (MQ=255)
                           gagTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGAc                                                            >  1:3303328/1‑93 (MQ=255)
                             gTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTg                                                          >  2:2068573/1‑93 (MQ=255)
                             gTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTg                                                          <  1:2088548/93‑1 (MQ=255)
                               gaaTCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGa                                                        <  2:3441548/91‑1 (MQ=255)
                                 aTCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAAc                                                      <  1:3685184/93‑1 (MQ=255)
                                 aTCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAAc                                                      <  1:769029/93‑1 (MQ=255)
                                  tCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACg                                                     >  2:4112464/1‑93 (MQ=255)
                                  tCAAACGTCAGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCAATGACTGGAACg                                                     <  2:4209898/93‑1 (MQ=255)
                                   cAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGa                                                    >  1:4416650/1‑93 (MQ=255)
                                   cAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGa                                                    >  1:929922/1‑93 (MQ=255)
                                    aaaCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCATGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGAt                                                   >  1:1418076/1‑93 (MQ=255)
                                     aaCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGGCTTGTACGATc                                                  >  1:639821/1‑93 (MQ=38)
                                     aaCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATc                                                  >  1:3354319/1‑93 (MQ=255)
                                     aaCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATc                                                  >  2:1003603/1‑93 (MQ=255)
                                     aaCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGCAAATTCGAAGCAGATCCATGACTGGAACGATc                                                  >  1:907383/1‑93 (MQ=255)
                                      aCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCt                                                 <  2:1223703/93‑1 (MQ=255)
                                      aCGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCt                                                 <  1:446333/93‑1 (MQ=255)
                                       cGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTg                                                >  2:3618433/1‑93 (MQ=255)
                                        gTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGa                                               >  1:4385424/1‑93 (MQ=255)
                                         tCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGAt                                              >  2:1488896/1‑93 (MQ=255)
                                         tCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGAt                                              >  2:4106495/1‑93 (MQ=255)
                                         cgtGCCGGCTGACTCCGCACCGTACACTTCCACCATTGTTTTCCGGGTGCGTAAGGGAAATCCGAAGCAGATCCATGCCTGGAACGATCTGAt                                              <  1:145618/91‑1 (MQ=255)
                                           tGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTa                                            >  1:1148565/1‑93 (MQ=255)
                                              cGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAc                                         >  2:130097/1‑93 (MQ=255)
                                               ggATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAcc                                        >  2:1846747/1‑93 (MQ=255)
                                               ggATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAAcc                                        >  1:4014593/1‑93 (MQ=255)
                                                gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCg                                       <  2:4464804/93‑1 (MQ=255)
                                                gATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCg                                       >  2:906480/1‑93 (MQ=255)
                                                 aTAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCgg                                      <  1:530064/93‑1 (MQ=255)
                                                  tAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCggg                                     >  1:2408286/1‑93 (MQ=255)
                                                   aaCTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGgt                                    <  1:4046041/93‑1 (MQ=255)
                                                    aCTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGGCTGGAACGATCTGATTAAACCGGgtg                                   <  1:654214/93‑1 (MQ=255)
                                                    aCTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGgtg                                   >  1:3072445/1‑93 (MQ=255)
                                                     cTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCATATCCATGACTGGAACGATCTGATTAAACCGGgtgt                                  <  2:253555/93‑1 (MQ=255)
                                                     cTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGgtgt                                  <  2:1465277/93‑1 (MQ=255)
                                                     cTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGgtgt                                  <  1:2491647/93‑1 (MQ=255)
                                                      tCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGtt                                 >  2:554121/1‑93 (MQ=255)
                                                      tCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGtt                                 >  1:798776/1‑93 (MQ=255)
                                                      tCCGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGtt                                 <  1:4623536/93‑1 (MQ=255)
                                                      cccGCCCCGTACACTGCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGtt                                 <  1:4623537/92‑1 (MQ=38)
                                                        cGCACCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTc                               <  1:1648180/93‑1 (MQ=255)
                                                           aCCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCggg                            >  2:1650166/1‑92 (MQ=255)
                                                           aCCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGt                            >  2:1990107/1‑93 (MQ=255)
                                                           aCCGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGCAGTTGCggg                            >  1:1670081/1‑92 (MQ=37)
                                                            ccGTACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTg                           >  1:2912006/1‑93 (MQ=255)
                                                             cGTACACTTCCCCCAATGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGa                          <  2:3679760/93‑1 (MQ=255)
                                                               tACCCTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATc                        <  2:4047510/93‑1 (MQ=255)
                                                               tACACTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATc                        <  2:4562625/93‑1 (MQ=255)
                                                                acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                       <  1:3040004/93‑1 (MQ=255)
                                                                acacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCa                       >  2:4419490/1‑93 (MQ=255)
                                                                 cacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCAc                      <  1:2191347/93‑1 (MQ=255)
                                                                 cacTTCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCAc                      <  2:449690/93‑1 (MQ=255)
                                                                    ttCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGcc                   <  1:106919/93‑1 (MQ=255)
                                                                     tCCACCATTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGACt                  >  1:2808705/1‑93 (MQ=255)
                                                                      ccaccaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTa                 >  1:1810977/1‑93 (MQ=255)
                                                                        accaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAAt               <  2:2221703/93‑1 (MQ=255)
                                                                         ccaTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATc              >  1:3946088/1‑93 (MQ=255)
                                                                          caTTGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATcc             <  2:3733291/93‑1 (MQ=255)
                                                                            ttGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGa           >  2:4377213/1‑93 (MQ=255)
                                                                            ttGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGa           <  2:3082465/93‑1 (MQ=255)
                                                                            ttGTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGa           <  1:448329/93‑1 (MQ=255)
                                                                              gTTTTCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaa         <  2:3880283/93‑1 (MQ=255)
                                                                               ttttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGaaaa        <  1:2695649/93‑1 (MQ=255)
                                                                                tttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAg       <  2:751038/93‑1 (MQ=255)
                                                                                 ttCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGc      <  1:1184819/93‑1 (MQ=255)
                                                                                  tCCTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGct     <  2:680430/93‑1 (MQ=255)
                                                                                   ccTGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGctc    >  2:2139562/1‑93 (MQ=255)
                                                                                     tGGTGCGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCtg  <  1:2176766/93‑1 (MQ=255)
                                                                                         |                                                                                        
ATTGCGGAACGCGGGCGGATTGATAAAGAGTGGATCAAACGTCTGCCGGATAACTCCGCACCGTACACTTCCACCATTGTTTTCCTGGTACGTAAGGGAAATCCGAAGCAGATCCATGACTGGAACGATCTGATTAAACCGGGTGTTTCGGTGATCACGCCTAATCCGAAAAGCTCTG  >  NC_000913/4106684‑4106861

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: