Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,672,424 T→C L163P (CTC→CCC)  b1601 → putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,672,4240TC98.7% 797.0 / ‑0.5 226L163P (CTC→CCC) b1601putative transport protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (125/98);  total (127/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

Only 100 of 237 total aligned reads displayed.

GATGACGTTGCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCTCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTG  >  NC_000913/1672335‑1672513
                                                                                         |                                                                                         
 aTGACGTTGCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGaa                                                                                       <  2:846823/93‑1 (MQ=255)
 aTGACGTTGCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGaa                                                                                       >  1:5088824/1‑93 (MQ=255)
 aTGACGTTGCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGaa                                                                                       >  2:39855/1‑93 (MQ=255)
 aTGACGTTGCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGaa                                                                                       <  1:2258789/93‑1 (MQ=255)
 aTGACGTTGCTCACCAACTTATTGACGCAGTTAGCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGaa                                                                                       <  2:5938957/93‑1 (MQ=255)
  tGACGTTGCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAg                                                                                      >  2:1836776/1‑93 (MQ=255)
      gTTGCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGcc                                                                                  >  1:2986618/1‑93 (MQ=255)
        tGCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCcac                                                                                <  1:3187658/93‑1 (MQ=255)
         gCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCcaca                                                                               >  2:1453882/1‑93 (MQ=255)
          cTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACaa                                                                              <  1:3438051/93‑1 (MQ=255)
           tCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAAt                                                                             >  1:2345763/1‑93 (MQ=255)
               cAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGcc                                                                         >  2:1974486/1‑93 (MQ=255)
                aaCTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGccc                                                                        >  1:2855917/1‑93 (MQ=255)
                    tATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGaa                                                                    >  2:111850/1‑93 (MQ=255)
                    tATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACAGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGaa                                                                    >  1:2385714/1‑93 (MQ=255)
                      ttGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGaaaa                                                                  >  1:5841923/1‑93 (MQ=255)
                      ttGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGaaaa                                                                  <  1:2318666/93‑1 (MQ=255)
                       tGACGCAGTTATCTAATTCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAAt                                                                 >  2:4647190/1‑93 (MQ=255)
                       tGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAAt                                                                 <  1:5144815/93‑1 (MQ=255)
                       tGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAAt                                                                 <  2:244871/93‑1 (MQ=255)
                       tGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAAt                                                                 <  1:1557783/93‑1 (MQ=255)
                          cGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTc                                                              <  2:3151773/93‑1 (MQ=255)
                           gCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCa                                                             >  1:6090253/1‑93 (MQ=255)
                           gCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCa                                                             <  2:4916234/93‑1 (MQ=255)
                           gCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCa                                                             >  2:5298011/1‑93 (MQ=255)
                           gCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCa                                                             >  2:1497345/1‑93 (MQ=255)
                            cAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAg                                                            <  2:5820017/93‑1 (MQ=255)
                             aGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGc                                                           >  2:3122573/1‑93 (MQ=255)
                              gTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCa                                                          >  1:416129/1‑93 (MQ=255)
                              gTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCa                                                          >  2:501257/1‑93 (MQ=255)
                              gTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCa                                                          >  1:1443950/1‑93 (MQ=255)
                               ttATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACTATTGCCCGGAAAATTTCAGCaa                                                         >  1:3756502/1‑93 (MQ=255)
                                tATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCaaa                                                        <  1:4361846/93‑1 (MQ=255)
                                 atttAATGCCATGTCATCCATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAACTTCAGCcaaa                                                       >  1:4558549/4‑92 (MQ=255)
                                 aTCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAat                                                       <  1:6110368/93‑1 (MQ=255)
                                 aTCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAat                                                       >  2:6102406/1‑93 (MQ=255)
                                  tCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAatg                                                      >  1:1122142/1‑93 (MQ=255)
                                  tCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAatg                                                      >  1:2894075/1‑93 (MQ=255)
                                   cTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAatga                                                     <  2:1718379/93‑1 (MQ=255)
                                   cTAATGCCACGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAatga                                                     >  1:3870492/1‑93 (MQ=255)
                                     aaTGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCTCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAatgatg                                                   <  1:2193635/93‑1 (MQ=255)
                                     aaTGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAatgatg                                                   <  2:2391135/93‑1 (MQ=255)
                                     aaTGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAatgatg                                                   <  2:715296/93‑1 (MQ=255)
                                      aTGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATgg                                                  <  1:637456/93‑1 (MQ=255)
                                           aTGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGACACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGt                                             >  1:838317/1‑93 (MQ=255)
                                            tGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTc                                            <  2:2282589/93‑1 (MQ=255)
                                             gTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTcc                                           <  1:4666428/93‑1 (MQ=255)
                                             gTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTcc                                           <  2:4069993/93‑1 (MQ=255)
                                              tcatcaATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCg                                          >  2:3333655/1‑93 (MQ=255)
                                              tcatcaATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCg                                          >  1:1650486/1‑93 (MQ=255)
                                              tcatcaATATTTTTATTGCTGCTGACGGTGCTGTTTATGCGGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCg                                          <  1:5815715/93‑1 (MQ=255)
                                                atcaATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGt                                        <  1:5404475/93‑1 (MQ=255)
                                                 tcaATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGtt                                       <  2:1347698/93‑1 (MQ=255)
                                                 tcaATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGtt                                       >  1:4023114/1‑93 (MQ=255)
                                                 tcaATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGtt                                       <  2:5873340/93‑1 (MQ=255)
                                                  caATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTg                                      >  1:4448933/1‑93 (MQ=255)
                                                  caATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCAGGTTg                                      >  2:3271539/1‑93 (MQ=255)
                                                    atatTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGaa                                    >  1:1215043/1‑93 (MQ=255)
                                                    atatTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGaa                                    <  2:440465/93‑1 (MQ=255)
                                                    atatTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGaa                                    <  1:1855844/93‑1 (MQ=255)
                                                    atatTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGaa                                    <  1:5243014/93‑1 (MQ=255)
                                                     tatTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAg                                   >  2:925599/1‑93 (MQ=255)
                                                     tatTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAg                                   >  2:4195072/1‑93 (MQ=255)
                                                     tatTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAg                                   >  2:2805718/1‑93 (MQ=255)
                                                     tatTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAg                                   >  2:4673685/1‑93 (MQ=255)
                                                      atTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAgg                                  >  1:3331710/1‑93 (MQ=255)
                                                       tttttATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCGGCAAATGATGGCGCGTCCGGTTGAAggg                                 >  2:455774/1‑93 (MQ=255)
                                                       tttttATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAggg                                 >  1:4820034/1‑93 (MQ=255)
                                                        ttttATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGa                                >  2:2977676/1‑93 (MQ=255)
                                                        ttttATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGa                                >  1:678144/1‑93 (MQ=255)
                                                          ttATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATg                              <  1:3635906/93‑1 (MQ=255)
                                                           tATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATgg                             >  2:4409010/1‑93 (MQ=255)
                                                           tATTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATgg                             <  1:92324/93‑1 (MQ=255)
                                                            aTTGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATggc                            >  1:4862287/1‑93 (MQ=255)
                                                             ttGCTGCTGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATggcg                           >  1:1314967/1‑93 (MQ=255)
                                                               gctgctGACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATggcggc                         >  2:1265416/1‑93 (MQ=255)
                                                                 tgctgACGGTGCTGTTTATGCTGCTCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGa                       >  2:2198502/1‑93 (MQ=255)
                                                                 tgctgACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGa                       >  2:3474099/1‑93 (MQ=255)
                                                                 tgctgACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGa                       >  2:1151608/1‑93 (MQ=255)
                                                                 tgctgACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGa                       >  1:3724736/1‑93 (MQ=255)
                                                                   ctgACGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGAtt                     <  1:4247384/93‑1 (MQ=255)
                                                                      aCGGTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCaa                  <  1:1456653/93‑1 (MQ=255)
                                                                        ggTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACg                <  2:6149442/93‑1 (MQ=255)
                                                                         gTGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGt               >  2:3622036/1‑93 (MQ=255)
                                                                          tGCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTg              >  1:4992629/1‑93 (MQ=255)
                                                                           gCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGc             <  1:1042513/93‑1 (MQ=255)
                                                                           gCTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGc             <  1:3467371/93‑1 (MQ=255)
                                                                            cTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCg            <  1:3053612/93‑1 (MQ=255)
                                                                            cTGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGAGGGCGGCGATTCAACGTGCg            >  2:3439582/1‑93 (MQ=255)
                                                                             tGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGa           >  1:5879197/1‑93 (MQ=255)
                                                                             tGTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGa           >  2:5163074/1‑93 (MQ=255)
                                                                              gTTTATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGAt          >  2:5376341/1‑93 (MQ=255)
                                                                               tttATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGGGCGAtt         >  2:1396044/1‑93 (MQ=255)
                                                                                ttATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTg        >  1:5646014/1‑93 (MQ=255)
                                                                                ttATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTg        <  1:3494893/93‑1 (MQ=255)
                                                                                 tATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGa       <  1:2059041/93‑1 (MQ=255)
                                                                                 tATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGa       >  1:5691950/1‑93 (MQ=255)
                                                                                 tATGCTGCCCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGa       <  2:2528714/93‑1 (MQ=255)
                                                                                    gctgccCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAg    <  1:2665571/93‑1 (MQ=255)
                                                                                      tgccCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAgtg  >  1:5304928/1‑93 (MQ=255)
                                                                                         |                                                                                         
GATGACGTTGCTCACCAACTTATTGACGCAGTTATCTAATGCCATGTCATCAATATTTTTATTGCTGCTGACGGTGCTGTTTATGCTGCTCGAAGTGCCACAATTGCCCGGAAAATTTCAGCAAATGATGGCGCGTCCGGTTGAAGGGATGGCGGCGATTCAACGTGCGATTGACAGTG  >  NC_000913/1672335‑1672513

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