Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,113,473 T→C 20.6% L61L (TTG→CTG) ‡ uspD → stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,113,4730TC20.6% 317.2 / 77.3 159L61L (TTG→CTG) uspDstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (62/64);  new base C (12/21);  total (74/85)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 166 total aligned reads displayed.

AGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAA  >  NC_000913/4113374‑4113571
                                                                                                   |                                                                                                  
gatcttAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt                                                                                                    <  2:3628141/95‑1 (MQ=255)
 gACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAtt                                                                                                   >  1:1778209/1‑100 (MQ=255)
   cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                 <  2:2705566/100‑1 (MQ=255)
    aTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATATACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                >  2:2866695/1‑100 (MQ=255)
     tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                               >  2:2235861/1‑100 (MQ=255)
       aTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaa                                                                                             <  2:2043791/100‑1 (MQ=255)
           cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                                                                                         >  2:321590/1‑100 (MQ=255)
            gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                        >  1:3430705/1‑100 (MQ=255)
            gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                        <  1:1520416/100‑1 (MQ=255)
            gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                        <  1:3914878/100‑1 (MQ=255)
            gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATa                                                                                        <  1:562742/100‑1 (MQ=255)
             cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATat                                                                                       <  1:3442273/100‑2 (MQ=255)
                aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatctt                                                                                    <  1:2456231/100‑5 (MQ=38)
                 ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg                                                                                   >  2:3172336/1‑100 (MQ=255)
                 ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                   <  1:1291306/100‑6 (MQ=38)
                  cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                  <  1:1604790/100‑1 (MQ=255)
                  cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                  >  1:1117377/1‑100 (MQ=255)
                  cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                  <  1:2160175/100‑7 (MQ=38)
                  cTGACGGTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACGTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                  <  1:1765902/100‑7 (MQ=37)
                   tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                  <  2:2122027/99‑7 (MQ=38)
                    gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGAt                                                                                <  2:3428016/100‑1 (MQ=255)
                    gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctgt                                                                                <  1:3592998/100‑9 (MQ=37)
                    gACGTTAATTCATATTGATGATGGCTGAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctgt                                                                                <  2:163070/100‑9 (MQ=37)
                     aCGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATa                                                                               <  1:3399102/100‑1 (MQ=255)
                      cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                              >  1:2909824/1‑100 (MQ=255)
                        ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                            <  2:418141/100‑1 (MQ=255)
                        ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                            <  2:569351/100‑1 (MQ=255)
                         tAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACaa                                                                           <  2:1448191/100‑1 (MQ=255)
                          aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACagg                                                                          >  2:204548/1‑98 (MQ=255)
                          aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAg                                                                          <  1:321590/100‑1 (MQ=255)
                             tCATATTGATGATGTCTTAAGCGATTTGTATCCGGGGATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCtt                                                                       <  2:4966/100‑2 (MQ=255)
                             tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                       >  2:2788075/1‑100 (MQ=255)
                               atatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGta                                                                     <  1:3677633/100‑1 (MQ=255)
                                 atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                   <  2:2197877/100‑1 (MQ=255)
                                 atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                   <  1:3910630/100‑1 (MQ=255)
                                 agtgatgaTGGCTAAAGCGAGTTGTATCCGGATATCTACTTCGCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                   >  1:2663548/3‑100 (MQ=255)
                                   tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                 <  1:734118/100‑1 (MQ=255)
                                         tgGCTTAAGCGAGTTGTATCCTGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACg                                                           <  1:926688/100‑1 (MQ=255)
                                         tgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACg                                                           >  2:2380248/1‑100 (MQ=255)
                                           gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                         >  1:317230/1‑100 (MQ=255)
                                           gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                         <  1:689099/100‑1 (MQ=255)
                                           gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAACAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                         >  2:1503633/1‑100 (MQ=255)
                                            cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaa                                                        <  1:480132/100‑1 (MQ=255)
                                              tAAGCGATTTGTATCCGGGTATCTACTTCCCTTCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTTACGaaa                                                        >  2:3464998/1‑98 (MQ=255)
                                               aaGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAat                                                     <  1:3615036/100‑1 (MQ=255)
                                                  cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                  <  2:1182279/100‑1 (MQ=255)
                                                   gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                 <  2:911334/100‑1 (MQ=255)
                                                    aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGAGGAAAAATATTCa                                                >  1:1018400/1‑100 (MQ=255)
                                                       tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                             >  2:1685651/1‑100 (MQ=255)
                                                        gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                            >  1:2831887/1‑100 (MQ=255)
                                                         tATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGc                                           >  2:410790/1‑100 (MQ=255)
                                                          aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                          >  2:3918049/1‑100 (MQ=255)
                                                           tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                         >  2:3094158/1‑100 (MQ=255)
                                                           tCCAGGTATCCACTAACCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAAAAAAAAAAAATATTCAATGGc                                           <  2:1053063/98‑1 (MQ=255)
                                                              gggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCTaa                                       <  1:1429715/99‑1 (MQ=255)
                                                              gggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                      <  1:2632210/100‑1 (MQ=255)
                                                              gggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                      <  2:2469965/100‑1 (MQ=255)
                                                                    cTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaaa                                >  2:758282/1‑100 (MQ=255)
                                                                    cTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaa                                 >  2:1288759/1‑99 (MQ=255)
                                                                     tactacCCTGCAACAGAAGATATTCTTCAATTGTTGAAAAATAAGTCGGATAACAAGCTGTAAAAAAAAAAAAAAAATATTCAATGGCCGAAGACAAAAc                               <  2:2341510/100‑1 (MQ=255)
                                                                     tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                               >  2:661474/1‑100 (MQ=255)
                                                                     tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                               >  2:490897/1‑100 (MQ=255)
                                                                      aCTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACt                              >  2:2933289/1‑100 (MQ=255)
                                                                       nttCCCTGCAACAGAAGATATTCTCCGATTGTTGAAGAATAACGCTGATAACAGGCTGTAGAAACTGACAAAAAATATGCAAGGACCGAATACAAAAc                               >  2:3292325/2‑98 (MQ=255)
                                                                        ttCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                            >  1:917139/1‑100 (MQ=255)
                                                                        ttCCCTGCAACAGAAGATATTCATCAATTGTCGAAGAGTAAGTCGGACAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                            <  2:2663531/100‑1 (MQ=255)
                                                                         tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                           <  2:1228946/100‑1 (MQ=255)
                                                                         tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                           <  1:1720950/100‑1 (MQ=255)
                                                                         tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                           >  1:3245069/1‑100 (MQ=255)
                                                                          cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                          >  2:1884855/1‑100 (MQ=255)
                                                                           ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                         <  1:410790/100‑1 (MQ=255)
                                                                            cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                        <  1:1120686/100‑1 (MQ=255)
                                                                                aaCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTg                      <  2:3087817/98‑1 (MQ=255)
                                                                                aaCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                    >  1:2380755/1‑100 (MQ=255)
                                                                                 aCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                   <  1:1834888/100‑1 (MQ=255)
                                                                                 aCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTAACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                    >  2:2531579/1‑99 (MQ=255)
                                                                                  cAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcg                  <  2:3685579/100‑1 (MQ=255)
                                                                                   agaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                 <  1:469807/100‑1 (MQ=255)
                                                                                   agaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                 <  1:1999935/100‑1 (MQ=255)
                                                                                    cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                >  2:443349/6‑100 (MQ=38)
                                                                                    cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                >  1:1671041/6‑100 (MQ=38)
                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg               >  1:2343200/5‑100 (MQ=38)
                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg               >  2:3373817/5‑100 (MQ=38)
                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg               >  1:1200941/5‑100 (MQ=38)
                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAAAATTCAATGACCGAAGACAAAACTGCGTATTGAACGCgg               >  2:1768296/5‑100 (MQ=38)
                                                                                     gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTAGAAACTGACGAAAAATATTCAATGGTCGAAGAAAAAACCGCGTATTAAAagc                 >  2:385884/5‑98 (MQ=25)
                                                                                     aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg               <  1:2309095/100‑1 (MQ=255)
                                                                                     aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg               <  1:876826/100‑1 (MQ=255)
                                                                                      agaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga              <  2:810324/100‑1 (MQ=255)
                                                                                       gaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag             >  1:1937014/1‑100 (MQ=255)
                                                                                       cttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag             >  2:1033371/3‑100 (MQ=255)
                                                                                         tatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa           >  1:1023416/1‑100 (MQ=255)
                                                                                          atTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaaa          >  1:1980714/1‑100 (MQ=255)
                                                                                            tcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg        >  1:2679548/1‑100 (MQ=255)
                                                                                              ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc      >  1:3208652/1‑100 (MQ=255)
                                                                                               tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg     >  1:2281435/1‑100 (MQ=255)
                                                                                                cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAATACAAAACTGCGTATTGAACGCGGAGAAATGCCgg    >  2:361356/1‑100 (MQ=255)
                                                                                                cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg    >  2:2855733/1‑100 (MQ=255)
                                                                                                cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg    >  2:1089311/1‑100 (MQ=255)
                                                                                                 aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa   <  2:3430682/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
AGACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAA  >  NC_000913/4113374‑4113571

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑