Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,113,475 G→A 15.2% L61L (TTG→TTA) ‡ uspD → stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,113,4750GA15.2% 489.2 / 63.4 196L61L (TTG→TTA) uspDstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base G (99/67);  new base A (16/14);  total (115/81)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

Only 100 of 219 total aligned reads displayed.

ACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACA  >  NC_000913/4113376‑4113574
                                                                                                   |                                                                                                   
aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                     <  2:1487085/100‑1 (MQ=255)
aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                     <  1:388731/100‑1 (MQ=255)
 cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                    >  1:1176460/1‑100 (MQ=255)
 cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                    <  1:1891675/100‑1 (MQ=255)
   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                  >  1:963961/1‑100 (MQ=255)
   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                  >  1:2985896/1‑100 (MQ=255)
   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                  <  2:2433925/100‑1 (MQ=255)
   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                  >  2:3660859/1‑100 (MQ=255)
   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                  <  2:1048298/100‑1 (MQ=255)
       gACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaaga                                                                                              <  2:2300595/100‑1 (MQ=255)
       gACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                               <  2:1442215/99‑1 (MQ=255)
        aCGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaagaa                                                                                             <  2:74776/100‑1 (MQ=255)
         cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                                                                                            <  2:1131665/100‑1 (MQ=255)
          gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                           <  1:2037803/100‑1 (MQ=255)
          gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                           <  1:2736007/100‑1 (MQ=255)
           cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATat                                                                                          <  1:1911777/100‑2 (MQ=255)
            tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                         <  1:3635091/100‑1 (MQ=255)
            tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                         <  1:2304630/100‑1 (MQ=255)
             cACCTGATGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATat                                                                                          <  1:179498/98‑2 (MQ=255)
             cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGt                                                                                        <  1:4085124/100‑1 (MQ=255)
             cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGt                                                                                        >  2:1193542/1‑100 (MQ=255)
             cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGt                                                                                        >  1:2556291/1‑100 (MQ=255)
              aCCTGACTTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTc                                                                                       <  1:4041605/100‑1 (MQ=255)
              aCCTGACGTTAATTCATATTGATGATGTCTTAAGCGAGTTGTATCCGGGTATCTACTACCCTGAAACAGAAGATATTCTTCAATTGTTGAAGAATAAGt                                                                                        <  2:4166354/99‑1 (MQ=255)
               ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg                                                                                      >  1:707490/1‑100 (MQ=255)
               ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                      <  1:470023/100‑6 (MQ=38)
                cTGACTTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                      <  2:3773713/99‑6 (MQ=38)
                 tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGa                                                                                    <  1:680726/100‑1 (MQ=255)
                  gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctgt                                                                                   <  1:838467/100‑9 (MQ=37)
                   aCGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATa                                                                                  >  1:2604887/1‑100 (MQ=255)
                   aCGTTAATTCATATTGAGGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTCCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctgtt                                                                                  <  2:610540/100‑10 (MQ=25)
                    cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                                 >  1:3639440/1‑100 (MQ=255)
                    cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATa                                                                                  >  2:2271568/1‑99 (MQ=255)
                     gTTAATTTCTATTTCTCATTGCTTAATAGAGTTATATCCTGATATCTACTTCCCTTCCACAGAAGATATTCTTCAAAAATTCAAGAATatcttc                                                                                      <  2:3748612/94‑6 (MQ=255)
                      ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                               >  2:1897495/1‑100 (MQ=255)
                           tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGGAGAATAAGTGGGATAACAAGCt                                                                           >  2:1377593/1‑99 (MQ=255)
                           tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                          <  1:4150175/100‑1 (MQ=255)
                           tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                          <  2:580118/100‑1 (MQ=255)
                             atatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGta                                                                        <  1:2126025/100‑1 (MQ=255)
                               atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                      >  2:1421848/1‑100 (MQ=255)
                                ttGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaa                                                                     >  1:3489932/1‑100 (MQ=255)
                                ttGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaa                                                                     >  1:142092/1‑100 (MQ=255)
                                  gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGTAGAATAAGTCGGATAACAAGCTGTATAAAc                                                                   >  2:3573673/1‑100 (MQ=255)
                                   atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                  <  2:2012297/100‑1 (MQ=255)
                                    tgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                                 >  2:3861158/1‑100 (MQ=255)
                                     gatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGa                                                                >  1:3265994/1‑100 (MQ=255)
                                       tgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACg                                                              >  2:528439/1‑100 (MQ=255)
                                          cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaa                                                           <  1:1074769/100‑1 (MQ=255)
                                          cTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaaa                                                           >  1:156242/1‑100 (MQ=255)
                                               gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGagatcggaaga                                                      >  1:906365/1‑90 (MQ=255)
                                                cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                     <  1:1193547/100‑1 (MQ=255)
                                                cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                     <  1:1447836/100‑1 (MQ=255)
                                                   gTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTGGGATAACAAGCTGTATAAAGTGACGAAAGATGTCCaa                                                  >  2:3824114/1‑100 (MQ=255)
                                                   gTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCaa                                                  >  1:4054648/1‑100 (MQ=255)
                                                     tGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATg                                                >  1:2653891/1‑100 (MQ=255)
                                                      gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                               >  2:3688461/1‑100 (MQ=255)
                                                        aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                             <  2:30107/100‑1 (MQ=255)
                                                        aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                             >  2:2709189/1‑100 (MQ=255)
                                                         tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                            <  2:3502204/100‑1 (MQ=255)
                                                            gggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                         <  1:3660886/100‑1 (MQ=255)
                                                            gggTATCCACTTCCCTGCAACAGAAGATAATCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                         <  2:2237442/100‑1 (MQ=255)
                                                             ggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGa                                        <  1:1634566/100‑1 (MQ=255)
                                                              gTATCTACTTCCCTACAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGAc                                       >  1:980475/1‑100 (MQ=255)
                                                               tATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACa                                      >  1:636476/1‑100 (MQ=255)
                                                                  gtACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGCATAAGTCGGATAACATGTTGTATAAACTGACGAAAAAGATTCAATGGCCGAAGACa                                      >  2:376569/2‑97 (MQ=255)
                                                                   tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                  >  1:2219001/1‑100 (MQ=255)
                                                                    aCTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACt                                 <  1:3595157/100‑1 (MQ=255)
                                                                     cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                                >  2:2353312/1‑100 (MQ=255)
                                                                       tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                              <  1:2266044/100‑1 (MQ=255)
                                                                        cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                             >  1:1614821/1‑100 (MQ=255)
                                                                        cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                             <  1:2153203/100‑1 (MQ=255)
                                                                          cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                              >  2:3992046/1‑97 (MQ=255)
                                                                          cTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAt                           >  2:1141006/1‑100 (MQ=255)
                                                                           tGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTAtt                          >  1:2859481/1‑100 (MQ=255)
                                                                            gCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTg                         >  2:2830288/1‑100 (MQ=255)
                                                                            gCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTg                         >  1:1318055/1‑100 (MQ=255)
                                                                             ttatccgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGa                        >  2:1274018/11‑100 (MQ=37)
                                                                                  gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcg                     >  1:899021/1‑98 (MQ=255)
                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGCCAAAACTGCGTATTAAACGCgg                  >  2:1394292/5‑100 (MQ=38)
                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  >  1:879286/5‑100 (MQ=38)
                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  >  1:2103716/5‑100 (MQ=39)
                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAATCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  >  1:2925470/5‑100 (MQ=38)
                                                                                    agaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  <  2:4079923/99‑1 (MQ=255)
                                                                                    agaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 <  1:79707/100‑1 (MQ=255)
                                                                                    acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 >  2:3614544/4‑100 (MQ=38)
                                                                                    acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 >  2:1645656/4‑100 (MQ=38)
                                                                                    acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 >  1:3027905/4‑100 (MQ=38)
                                                                                     gaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag                <  2:1327122/100‑1 (MQ=255)
                                                                                       tatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa              >  2:3050925/1‑100 (MQ=255)
                                                                                       tatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGAAAAAACTGCGTATTGAACGCGGAGaa              >  1:3816154/1‑100 (MQ=255)
                                                                                           cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGc          >  2:1794351/1‑100 (MQ=255)
                                                                                            ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc         <  1:1299195/100‑1 (MQ=255)
                                                                                            ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc         <  1:1982606/100‑1 (MQ=255)
                                                                                              cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg       >  2:4143761/1‑100 (MQ=255)
                                                                                              cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg        >  2:3926690/1‑99 (MQ=255)
                                                                                               aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGGAGACAAAACTGCGTATAGAACGCGGAGAAATGc          >  2:883903/1‑96 (MQ=255)
                                                                                               aaTTGGTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAAGGCCgg       <  2:199245/99‑1 (MQ=255)
                                                                                                 ttgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaaa    >  2:849733/1‑100 (MQ=255)
                                                                                                  tgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGGAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAac   >  2:3320668/1‑100 (MQ=255)
                                                                                                   gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAaca  >  1:1129661/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACA  >  NC_000913/4113376‑4113574

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