Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,113,474 T→A 16.9% L61* (TTG→TAG) ‡ uspD → stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,113,4740TA16.9% 525.9 / 90.1 218L61* (TTG→TAG) uspDstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base T (92/88);  new base A (22/15);  total (115/103)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 232 total aligned reads displayed.

GACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113375‑4113573
                                                                                                    |                                                                                                   
    tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                 <  2:717977/100‑1 (MQ=255)
    tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                 >  1:1125545/1‑100 (MQ=255)
    tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                 >  1:434476/1‑100 (MQ=255)
    tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                 >  1:3877462/1‑100 (MQ=255)
     aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAt‑tgttg                                                                                                 <  2:1351995/99‑1 (MQ=255)
       tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTgaag                                                                                              <  1:733612/100‑1 (MQ=255)
         aCGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTgaagaa                                                                                            <  1:863671/100‑1 (MQ=255)
         aCGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAC‑AATTgaagaa                                                                                            <  2:3302235/100‑1 (MQ=255)
          cGCTCACCTGACTTTAAGTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAAt                                                                                           <  2:3012454/100‑1 (MQ=255)
            cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAC‑AATTGAAGAATat                                                                                         <  2:3595860/100‑2 (MQ=255)
            cTCACCTGACATTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAC‑AATTGAAGAATat                                                                                         <  2:3959634/100‑2 (MQ=255)
             tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCAGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAg                                                                                        <  2:2923078/100‑1 (MQ=255)
              caacTGACGTTAATTCATAGTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACGTCCTTTAAACAGAAGATATGCTTCAAC‑AATTCAAGAATa                                                                                          <  1:2495134/94‑1 (MQ=255)
              cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAC‑AATTGAAGAATatct                                                                                       <  2:1415794/100‑4 (MQ=38)
               aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTc                                                                                      >  2:3319474/1‑100 (MQ=255)
                ccTGACTTTAATTCATATTGATGATGCCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAC‑AATTGAAGAATatcttc                                                                                     <  1:1884523/100‑6 (MQ=37)
                ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCg                                                                                     <  1:611450/100‑1 (MQ=255)
                ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCg                                                                                     <  2:3029020/100‑1 (MQ=255)
                 cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCgg                                                                                    <  2:1017226/100‑1 (MQ=255)
                 cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCgg                                                                                    <  2:3898129/100‑1 (MQ=255)
                 cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCgg                                                                                    >  2:1993827/1‑100 (MQ=255)
                 cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAC‑AATTGAAGAATatcttct                                                                                    <  2:2333843/100‑7 (MQ=38)
                  tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCgg                                                                                    <  1:147037/99‑1 (MQ=255)
                  tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGa                                                                                   >  1:816380/1‑100 (MQ=255)
                  tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGa                                                                                   <  2:192750/100‑1 (MQ=255)
                   gACGTTAATTCATATTGATGATGGCTTGAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGAt                                                                                  <  2:2121987/100‑1 (MQ=255)
                   gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGAt                                                                                  <  2:3801670/100‑1 (MQ=255)
                     cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATaa                                                                                >  2:4194421/1‑100 (MQ=255)
                     cGTTAATTCATAATGAAGATGCCTTAAGCGAGTTGTATCCGGAAATCTACTTCCCTGCAACAGAAGATATACTTCAAT‑TGTTGAAGAATAAGTCgg                                                                                    <  1:1520341/96‑1 (MQ=255)
                       ttAATTCATATTGATGATGGTTTAAGCTAGTTGTATTCTTGTATTTACTTTTCTGCAACAGAAGATATTCTTTAAC‑AATTGAAGAATatct                                                                                       <  1:3841335/91‑4 (MQ=255)
                         aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAg                                                                            >  2:347394/1‑100 (MQ=255)
                         aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑AGTTGAAGAATAAGTCGGATAACAAg                                                                            >  2:604228/1‑100 (MQ=255)
                               tatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGtat                                                                      >  2:54108/1‑100 (MQ=255)
                               tatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGtat                                                                      >  1:4113267/1‑100 (MQ=255)
                                atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                     <  1:1942226/100‑1 (MQ=255)
                                 ttGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATaa                                                                    >  1:1543853/1‑100 (MQ=255)
                                  tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGGTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                     <  1:192963/98‑1 (MQ=255)
                                  tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                   <  1:882077/100‑1 (MQ=255)
                                   gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                  >  2:655176/1‑100 (MQ=255)
                                     tgatgGCTTAAGCGAGTTGTATCCGGGGATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                                <  1:306043/100‑1 (MQ=255)
                                          gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                           >  1:2778934/1‑100 (MQ=255)
                                           cTTGAGCGAGTTGTGTTCGAGTATCCAGGTCCCGGCTACAGAAGATATTCTGCAAT‑TGTTGATGAGGAAGTCGGATAACAAGCGGTATAAAATGACGaa                                                           >  1:2564403/1‑99 (MQ=255)
                                           cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaa                                                          <  1:474638/100‑1 (MQ=255)
                                            ttAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaaa                                                         <  2:1735172/100‑1 (MQ=255)
                                              aaGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAat                                                       <  2:4177458/100‑1 (MQ=255)
                                                gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                     >  2:3754071/1‑100 (MQ=255)
                                                  gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGGATAATATTc                                                   <  1:3443081/100‑1 (MQ=255)
                                                   aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                  <  1:2701704/100‑1 (MQ=255)
                                                      tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                               >  1:993171/1‑100 (MQ=255)
                                                       gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                              <  2:1655623/100‑1 (MQ=255)
                                                        tATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGc                                             <  1:2485373/100‑1 (MQ=255)
                                                        tATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGc                                             <  2:1217185/100‑1 (MQ=255)
                                                         gtCCGGATATCTACTTCGCTGCTGCAGATCATATTTTTCGATGGGCTGAA‑CATCAGCGGGGTAGCTAGCTGCATTAACTGACGAAAAATATTCAATAGcc                                            >  1:1676334/2‑100 (MQ=255)
                                                         aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                            >  2:1561961/1‑100 (MQ=255)
                                                          tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                           >  2:1492023/1‑100 (MQ=255)
                                                           gcGGGTATCTACTTCCCTGCACCCGAAGAGATTGTCCAAT‑ATTTTTTCAATCATTCGGATAACAAGCTGTAACAACTGACGAGAAGTATTCAAt                                                >  1:3694856/2‑94 (MQ=255)
                                                             gggTATCTAGTTCCATGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                        <  1:353229/100‑1 (MQ=255)
                                                                tATCTACTTCCCTGCAACAGAAGATATTCTGCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATAGCCGAAGACa                                     >  1:2896981/1‑100 (MQ=255)
                                                                   cTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaaa                                  >  2:187175/1‑100 (MQ=255)
                                                                    tACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                 <  1:1100687/100‑1 (MQ=255)
                                                                     aCTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGACAGAAAACAAAACt                                >  1:3814695/1‑100 (MQ=255)
                                                                       ttCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                              <  1:2122986/100‑1 (MQ=255)
                                                                        tCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                             >  1:3134699/1‑100 (MQ=255)
                                                                        tCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGACGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaaaat                                >  1:4274824/1‑95 (MQ=255)
                                                                         cccTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGATAAATATTCAATGGCCGATGACAAAACTGcct                            <  2:2964407/100‑3 (MQ=255)
                                                                         cccTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                            >  2:3180778/1‑100 (MQ=255)
                                                                         cccTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                            <  1:2175993/100‑1 (MQ=255)
                                                                          ccTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                           >  2:2399647/1‑100 (MQ=255)
                                                                           cTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                          >  1:2327610/1‑100 (MQ=255)
                                                                             gCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTg                        >  2:271578/1‑100 (MQ=255)
                                                                                aCAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                     <  1:168804/100‑1 (MQ=255)
                                                                                  ccgacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                   >  1:2029479/7‑100 (MQ=38)
                                                                                  agaagaTATTCTTGAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                   <  1:1069173/100‑1 (MQ=255)
                                                                                  agaagaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                   >  1:694021/1‑100 (MQ=255)
                                                                                  agaagaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                   <  1:2229394/100‑1 (MQ=255)
                                                                                  agaagaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                   <  2:2650462/100‑1 (MQ=255)
                                                                                   gaagaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  <  1:2584328/100‑1 (MQ=255)
                                                                                    gacttatTCTTCAAC‑AATTGAAGAATCAGTCCGATACCAACCTGTATAAACTGACGAAAAATCTCCAACGCCCGAAGCCAACCCTCCGTATTGAAcgc                   >  1:2408224/5‑98 (MQ=16)
                                                                                    gacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  >  1:2043207/5‑99 (MQ=38)
                                                                                    gacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  >  2:1095057/5‑99 (MQ=38)
                                                                                    gacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 >  1:3234989/5‑100 (MQ=38)
                                                                                    gacttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 >  1:1217439/5‑100 (MQ=38)
                                                                                     agaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                >  2:851494/1‑100 (MQ=255)
                                                                                     agaTATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                >  2:3525399/1‑100 (MQ=255)
                                                                                     acttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                >  2:2432515/4‑100 (MQ=38)
                                                                                     acttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                >  2:4042574/4‑100 (MQ=38)
                                                                                     acttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                >  1:1487762/4‑100 (MQ=38)
                                                                                      cttatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag               >  1:3664195/3‑100 (MQ=255)
                                                                                        tatTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAAGCGGAGaa             >  1:4070321/1‑100 (MQ=255)
                                                                                        tatTCTTCAAC‑AATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa             >  2:4256140/1‑100 (MQ=255)
                                                                                          ttcttcAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAAt           >  2:1813318/1‑100 (MQ=255)
                                                                                            cttcAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg          >  2:476754/1‑99 (MQ=255)
                                                                                             ttcAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc        >  2:2734225/1‑100 (MQ=255)
                                                                                              tcAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg       <  1:2820547/100‑1 (MQ=255)
                                                                                               cAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg      <  1:4194391/100‑1 (MQ=255)
                                                                                                aaT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa     <  1:3410365/100‑1 (MQ=255)
                                                                                                aaT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa     <  2:3822083/100‑1 (MQ=255)
                                                                                                aaT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa     <  1:3754047/100‑1 (MQ=255)
                                                                                                 aT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa    >  1:4244273/1‑100 (MQ=255)
                                                                                                    tgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAac  <  1:2458642/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
GACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAT‑TGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113375‑4113573

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: ‑