Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,113,475 G→A 20.8% 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,4750GA20.8% 331.7 / 78.2 159L61L (TTG→TTA) uspDstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base G (62/64);  new base A (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 = 6.98e-01

Only 100 of 167 total aligned reads displayed.

CATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACA  >  NC_000913/4113377‑4113574
                                                                                                  |                                                                                                   
 aTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATATACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                   >  2:2866695/1‑100 (MQ=255)
    aTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaa                                                                                                <  1:3935200/100‑1 (MQ=255)
    aTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaa                                                                                                <  2:2043791/100‑1 (MQ=255)
     tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                               <  2:491178/100‑1 (MQ=255)
        cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAAt                                                                                            <  2:2046082/100‑1 (MQ=255)
         gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                           >  1:3430947/1‑100 (MQ=255)
         gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                           <  1:1520416/100‑1 (MQ=255)
         gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                           >  1:3430705/1‑100 (MQ=255)
         gCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATa                                                                                           <  1:562742/100‑1 (MQ=255)
          cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATat                                                                                          <  1:3181514/100‑2 (MQ=255)
           tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAg                                                                                         >  1:3242231/1‑100 (MQ=255)
             aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatctt                                                                                       <  2:481184/100‑5 (MQ=38)
             aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatctt                                                                                       <  2:1795820/100‑5 (MQ=38)
             aCCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatctt                                                                                       <  1:2456231/100‑5 (MQ=38)
              ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCg                                                                                      >  2:3172336/1‑100 (MQ=255)
              ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                      <  1:1291306/100‑6 (MQ=38)
              ccTGACGTTAATTCATATTGATGATCGCTTAAGCTAGTTGTATACGGGTCTATACTACCCTTCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                      <  2:1932944/100‑6 (MQ=21)
               cTGATGTTATTTTATATTGACGATTATTTTAGCGAATTGTATCCGGGTATATAGTCCACTGCAGCAGAAGATATTCTTCAACAATTGAataatatct                                                                                        <  2:393418/97‑4 (MQ=255)
               cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                     <  1:1604790/100‑1 (MQ=255)
               cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCgg                                                                                     >  1:1117377/1‑100 (MQ=255)
               cTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                     <  2:249975/100‑7 (MQ=38)
               cTGACGGTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACGTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                     <  1:1765902/100‑7 (MQ=37)
                tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctg                                                                                    <  1:3398503/100‑8 (MQ=37)
                 gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctg                                                                                    <  1:1674231/99‑8 (MQ=37)
                   cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTCTCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                                 >  2:3827640/1‑100 (MQ=255)
                   cGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATaa                                                                                 >  1:2909824/1‑100 (MQ=255)
                     ttACTTCAGATTTACGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTGGAAGAATAAGTCGGATAACa                                                                               <  2:1039332/100‑1 (MQ=255)
                     ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                               <  2:418141/100‑1 (MQ=255)
                     ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                               <  2:3651661/100‑1 (MQ=255)
                      tAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACaa                                                                              <  2:1448191/100‑1 (MQ=255)
                       aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAg                                                                             <  1:321590/100‑1 (MQ=255)
                        aTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc                                                                            >  2:2614567/1‑100 (MQ=255)
                        aTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGc                                                                            <  1:231746/100‑1 (MQ=255)
                          tCATATTGATGATGTCTTAAGCGATTTGTATCCGGGGATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCtt                                                                          <  2:4966/100‑2 (MQ=255)
                          tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                          >  2:2788075/1‑100 (MQ=255)
                            atatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGta                                                                        <  2:718959/100‑1 (MQ=255)
                            atatTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGta                                                                        <  1:3677633/100‑1 (MQ=255)
                              atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                      <  2:2197877/100‑1 (MQ=255)
                              atTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGtata                                                                      <  1:3910630/100‑1 (MQ=255)
                                tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                    <  1:734118/100‑1 (MQ=255)
                                   tgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTg                                                                 >  1:900129/1‑100 (MQ=255)
                                      tgGCTTAAGCGAGTTGTATCCTGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACg                                                              <  1:926688/100‑1 (MQ=255)
                                      tgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACg                                                              >  2:2380248/1‑100 (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)
                                              gCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAatat                                                      >  2:2124747/1‑100 (MQ=255)
                                               cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                     >  1:1182288/1‑100 (MQ=255)
                                               cGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATAtt                                                     <  2:1182279/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)
                                                     gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                               >  2:1644902/1‑100 (MQ=255)
                                                     gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAAAATTCAAt                                                 >  2:3191006/1‑98 (MQ=255)
                                                        tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                            >  2:3094158/1‑100 (MQ=255)
                                                          cGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGaa                                          <  2:2798288/100‑1 (MQ=255)
                                                           gggTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAg                                         <  1:2632210/100‑1 (MQ=255)
                                                            ggTATCTACTTCCCTGCAACAGAAGATCTCCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGaa                                          <  1:1430053/98‑1 (MQ=255)
                                                                 cTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaaa                                   >  2:758282/1‑100 (MQ=255)
                                                                 cTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaa                                    >  2:1288759/1‑99 (MQ=255)
                                                                  tgcTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                  >  1:2670999/3‑100 (MQ=255)
                                                                  tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                  >  2:2762407/1‑100 (MQ=255)
                                                                  tACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAAc                                  >  2:2357956/1‑100 (MQ=255)
                                                                   aCTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACt                                 >  2:2933289/1‑100 (MQ=255)
                                                                    nttCCCTGCAACAGAAGATATTCTCCGATTGTTGAAGAATAACGCTGATAACAGGCTGTAGAAACTGACAAAAAATATGCAAGGACCGAATACAAAAc                                  >  2:3292325/2‑98 (MQ=255)
                                                                     ttCCCTGCAACAGAAGATATTCATCAATTGTCGAAGAGTAAGTCGGACAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGc                               <  2:2663531/100‑1 (MQ=255)
                                                                      tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                              <  1:1720950/100‑1 (MQ=255)
                                                                      tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                              <  2:1228946/100‑1 (MQ=255)
                                                                       cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                             >  2:3385105/1‑100 (MQ=255)
                                                                       cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGt                             >  1:2067426/1‑100 (MQ=255)
                                                                        ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                            <  1:410790/100‑1 (MQ=255)
                                                                           gCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                            <  2:905648/97‑1 (MQ=255)
                                                                             aaCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGaa                       >  1:2380755/1‑100 (MQ=255)
                                                                              aCAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                      <  1:1834888/100‑1 (MQ=255)
                                                                               cAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcg                     >  2:934704/1‑100 (MQ=255)
                                                                                agaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgc                    <  1:1999935/100‑1 (MQ=255)
                                                                                 gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                   >  1:680315/1‑100 (MQ=255)
                                                                                  gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  >  1:1200941/5‑100 (MQ=38)
                                                                                  gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  >  2:3373817/5‑100 (MQ=38)
                                                                                  gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTAGAAACTGACGAAAAATATTCAATGGTCGAAGAAAAAACCGCGTATTAAAagc                    >  2:385884/5‑98 (MQ=25)
                                                                                  aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  <  1:876826/100‑1 (MQ=255)
                                                                                  aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  <  2:618510/100‑1 (MQ=255)
                                                                                  aagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                  <  1:2309095/100‑1 (MQ=255)
                                                                                   acttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 >  1:3150776/4‑100 (MQ=38)
                                                                                    ttccgatctcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                 <  1:3586540/91‑1 (MQ=255)
                                                                                    cttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag                >  2:1033371/3‑100 (MQ=255)
                                                                                    cttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag                >  1:1909085/3‑100 (MQ=255)
                                                                                      tatTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa              <  2:23130/100‑1 (MQ=255)
                                                                                      tatTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa              >  2:413412/1‑100 (MQ=255)
                                                                                      tatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaa              >  1:1023416/1‑100 (MQ=255)
                                                                                       atTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGaaa             >  1:1980714/1‑100 (MQ=255)
                                                                                         tcttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATg           >  1:2679548/1‑100 (MQ=255)
                                                                                            tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg        >  1:2281435/1‑100 (MQ=255)
                                                                                             cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgaagatcggaag       >  2:3586512/1‑90 (MQ=255)
                                                                                             cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg       >  2:1089311/1‑100 (MQ=255)
                                                                                              aaTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGa      <  2:3430682/100‑1 (MQ=255)
                                                                                               aTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa     <  2:3067939/100‑1 (MQ=255)
                                                                                                  gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAaca  >  2:1628721/1‑100 (MQ=255)
                                                                                                  gttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAaca  >  1:2461383/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAACA  >  NC_000913/4113377‑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: ‑