Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,202,618 T→C I95T (ATC→ACC)  ygiH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,202,6180TC98.2% 523.0 / ‑1.0 165I95T (ATC→ACC) ygiHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (70/92);  total (71/94)
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.37e-01

Only 100 of 173 total aligned reads displayed.

AAGGAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACATCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCACCGGAG  >  NC_000913/3202528‑3202707
                                                                                          |                                                                                         
aaGGAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCt                                                                                         <  1:1084021/93‑1 (MQ=255)
aaGGAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCt                                                                                         <  1:3423009/93‑1 (MQ=255)
aaGGAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCt                                                                                         <  2:3209824/93‑1 (MQ=255)
aaGGAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCt                                                                                         >  2:3938679/1‑93 (MQ=255)
 aGGAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGCCACACCTg                                                                                        <  2:4656809/93‑1 (MQ=255)
 aGGAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTg                                                                                        <  2:1402390/93‑1 (MQ=255)
 aGGAATGGTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTg                                                                                        <  1:313367/93‑1 (MQ=255)
  ggAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTgg                                                                                       >  2:1506421/1‑93 (MQ=255)
  ggAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTgg                                                                                       >  1:4364622/1‑93 (MQ=255)
    aaTGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGcc                                                                                     <  2:4549459/93‑1 (MQ=255)
     aTGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACATCTGGccc                                                                                    >  1:1592075/1‑93 (MQ=255)
      tGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCg                                                                                   >  2:4402396/1‑93 (MQ=255)
       gTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGt                                                                                  <  2:4195083/93‑1 (MQ=255)
         tacccGTCTGGGGCCCGTTTGAATTAGGTGTTATCCCCTTCTGGCTAGTCTCAAGTGCCATCTCTTCCTGTCTTAGACACACCTTGCCCGttt                                                                                <  2:2672226/91‑1 (MQ=255)
          gCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGtttt                                                                               <  2:1019108/93‑1 (MQ=255)
             cGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTttctt                                                                            >  1:3367369/1‑93 (MQ=255)
              gTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTttcttc                                                                           <  1:2399470/93‑1 (MQ=255)
               tCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCg                                                                          <  2:1411945/93‑1 (MQ=255)
                cTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCgg                                                                         <  1:1670355/93‑1 (MQ=255)
                cTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGACACACCTGGCCCGTTTTCTTCGGa                                                                        <  1:821957/93‑1 (MQ=255)
                 tGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGa                                                                        >  1:3323724/1‑93 (MQ=255)
                  ggggCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGAt                                                                       >  1:2566270/1‑93 (MQ=255)
                  ggggCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGAt                                                                       >  1:4621914/1‑93 (MQ=255)
                   gggCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGAtt                                                                      >  2:2005598/1‑93 (MQ=255)
                    ggCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGAttt                                                                     <  1:1926800/93‑1 (MQ=255)
                       gcgTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTAGACACACCTGGCCCGTTTTCTTCGGATTTaaa                                                                  <  2:3778572/93‑1 (MQ=255)
                         gTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAagg                                                                >  1:4235044/1‑93 (MQ=255)
                         gTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAagg                                                                <  2:1816260/93‑1 (MQ=255)
                         gTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAagg                                                                <  1:1983282/93‑1 (MQ=255)
                          tATGAATTAGGTGTCAGCCCCTTCTGTCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAagga                                                               <  2:4067070/93‑1 (MQ=255)
                            tGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAaggagg                                                             >  1:4093125/1‑93 (MQ=255)
                            tGAATTAGGTGTCAGCCCCTTCTGGCTAGACTTAATTGCCATCGCCGCCTGTCTTGGCCACACCTGGCCCGTTTTCTTCGGATTTAAaggagg                                                             >  1:4093126/1‑93 (MQ=255)
                             gAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAggagga                                                            <  1:1291767/93‑1 (MQ=255)
                             gAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCAGCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCGTCGGATTTTGAtgagga                                                            >  2:3145350/1‑93 (MQ=255)
                               aTTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTTGACCCACCTGGCCCGTTTTTTTCTGATTTAAAGGAGGgaa                                                          >  2:3570906/1‑93 (MQ=255)
                                ttAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGaaaa                                                         <  2:2590648/93‑1 (MQ=255)
                                ttAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGaaaa                                                         <  2:501498/93‑1 (MQ=255)
                                ttAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGaaaa                                                         >  1:3759737/1‑93 (MQ=255)
                                ttAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGaaaa                                                         <  2:738378/93‑1 (MQ=255)
                                 tAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAg                                                        <  2:4491908/93‑1 (MQ=255)
                                 tAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCAGCTCCGCCTGTCTTGGACACACCTGGCCCGTGTTCTTCGGATTTAAAGGAGGAAAAg                                                        <  2:1912808/93‑1 (MQ=255)
                                  aGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAgg                                                       >  1:2224851/1‑93 (MQ=255)
                                  aGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAgg                                                       <  1:4050258/93‑1 (MQ=255)
                                  aGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAgg                                                       <  1:1045762/93‑1 (MQ=255)
                                   ggTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGc                                                      <  2:2816246/93‑1 (MQ=255)
                                    gtgtCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCg                                                     >  1:562206/1‑93 (MQ=255)
                                     tgtCAGCCCCTTCTTGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGt                                                    >  1:2660213/1‑93 (MQ=255)
                                     tgtCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGt                                                    <  2:2551125/93‑1 (MQ=255)
                                          gCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTa                                               <  1:2355/93‑1 (MQ=255)
                                           ccccTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACATCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTAc                                              <  2:651852/93‑1 (MQ=255)
                                           ccccTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTAc                                              <  2:2263174/93‑1 (MQ=255)
                                            cccTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTAcc                                             >  1:2841829/1‑93 (MQ=255)
                                            cccTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTAcc                                             <  2:1757175/93‑1 (MQ=255)
                                             ccTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCg                                            >  1:3088106/1‑93 (MQ=255)
                                                 cTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCttt                                        >  1:2670616/1‑93 (MQ=255)
                                                 cTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCttt                                        <  2:3713229/93‑1 (MQ=255)
                                                  tGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCtttt                                       >  1:904017/1‑93 (MQ=255)
                                                    gCTAGTCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTgg                                     <  1:3827805/93‑1 (MQ=255)
                                                    gCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTgg                                     <  1:3827804/93‑1 (MQ=255)
                                                      tAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTg                                   >  2:7706/1‑93 (MQ=255)
                                                        ggCTTAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGcc                                 >  1:2625368/1‑93 (MQ=255)
                                                           ttAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATc                              <  2:4100652/93‑1 (MQ=255)
                                                           ttAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATc                              <  1:496988/93‑1 (MQ=255)
                                                            tAATTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCg                             >  1:3986782/1‑93 (MQ=255)
                                                             aaTTGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGc                            >  2:2339459/1‑93 (MQ=255)
                                                               ttGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCAc                          >  2:2764254/1‑93 (MQ=255)
                                                                tGCCATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCAcc                         >  2:2114993/1‑93 (MQ=255)
                                                                   cATCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCAt                      >  1:2207654/1‑93 (MQ=255)
                                                                    aTCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCAtt                     <  2:2007113/93‑1 (MQ=255)
                                                                    aTCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCAtt                     <  1:1481070/93‑1 (MQ=255)
                                                                     tCGCCGCCTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGTAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTg                    >  1:1298117/1‑93 (MQ=255)
                                                                       gccgccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGc                  <  2:887185/93‑1 (MQ=255)
                                                                       gccgccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGc                  >  1:2558447/1‑93 (MQ=255)
                                                                        ccgccTGTCTTGGACCCACCTGGCCCGTTTTCTTTGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCt                 <  1:1176930/93‑1 (MQ=255)
                                                                         cgccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTg                <  2:2206910/93‑1 (MQ=255)
                                                                         cgccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTg                <  2:2207054/93‑1 (MQ=255)
                                                                          gccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTgg               <  1:1898796/93‑1 (MQ=255)
                                                                          gccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTgg               >  2:3097399/1‑93 (MQ=255)
                                                                           ccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTTAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTggg              >  1:933590/1‑93 (MQ=255)
                                                                           ccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTggg              >  2:4682323/1‑93 (MQ=255)
                                                                           ccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTggg              <  1:4199169/93‑1 (MQ=255)
                                                                           ccTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTggg              <  2:1069126/93‑1 (MQ=255)
                                                                            cTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGa             >  1:1867584/1‑93 (MQ=255)
                                                                            cTGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGa             <  2:4109847/93‑1 (MQ=255)
                                                                             tGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGAt            <  1:3487516/93‑1 (MQ=255)
                                                                             tGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGAt            >  1:4443096/1‑93 (MQ=255)
                                                                             tGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGAt            >  1:4433889/1‑93 (MQ=255)
                                                                             tGTCTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGAt            >  1:4443095/1‑93 (MQ=255)
                                                                             tGTCTTGGACACACCTGGCCCGTTTTCTTAGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGAt            <  2:2100619/93‑1 (MQ=255)
                                                                                tttGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGAtctc         <  2:977581/92‑1 (MQ=255)
                                                                                cTTGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGAtctc         <  2:2442972/93‑1 (MQ=255)
                                                                                 ttGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGAtctcc        >  2:2943941/1‑92 (MQ=255)
                                                                                 ttGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCa        >  1:2529721/1‑93 (MQ=255)
                                                                                  tGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCAc       <  1:1808944/93‑1 (MQ=255)
                                                                                  tGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCAc       <  1:684102/93‑1 (MQ=255)
                                                                                  tGGACACACCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCAc       >  1:321126/1‑93 (MQ=255)
                                                                                    gACACACCTGGCCCTTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCACCg     <  2:1963080/93‑1 (MQ=255)
                                                                                      cacaCCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCACCGGa   >  1:4264170/1‑93 (MQ=255)
                                                                                       acacCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCACCGGAg  >  1:1687890/1‑93 (MQ=255)
                                                                                       acacCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCACCGGAg  <  2:4593042/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAGGAATGTTGCCCGTCTGGGGCGCGTATGAATTAGGTGTCAGCCCCTTCTGGCTAGGCTTAATTGCCATCGCCGCCTGTCTTGGACACATCTGGCCCGTTTTCTTCGGATTTAAAGGAGGAAAAGGCGTTGCTACCGCTTTTGGTGCCATCGCACCCATTGGCTGGGATCTCACCGGAG  >  NC_000913/3202528‑3202707

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