Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,113,475 G→A 18.9% 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,4750GA18.9% 406.8 / 81.7 185L61L (TTG→TTA) uspDstress‑induced protein
Reads supporting (aligned to +/- strand):  ref base G (79/70);  new base A (18/17);  total (98/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-03

Only 100 of 204 total aligned reads displayed.

ACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113376‑4113573
                                                                                                   |                                                                                                  
aCATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttg                                                                                                    <  1:352001/100‑1 (MQ=255)
 cATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgt                                                                                                   >  1:1129418/1‑100 (MQ=255)
  aTAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgtt                                                                                                  <  1:1952550/100‑1 (MQ=255)
   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 >  2:4197795/1‑100 (MQ=255)
   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 <  2:2592211/100‑1 (MQ=255)
   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 <  2:1311708/100‑1 (MQ=255)
   tAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAAttgttg                                                                                                 >  1:3483599/1‑100 (MQ=255)
   tAATGACGCTCACCTGACATTAATTCATATTGAAGATCGCTTAAGCGAGTTGTATACGGGTATATACTTCCCTGCAACAGAAGATATTCTTCACttg                                                                                                    >  1:705508/1‑97 (MQ=255)
    aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTga                                                                                                >  1:2690945/1‑100 (MQ=255)
    aaTGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGCTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTga                                                                                                >  1:22702/1‑100 (MQ=255)
      tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                              <  2:83308/100‑1 (MQ=255)
      tGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTATATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaag                                                                                              >  2:1187729/1‑100 (MQ=255)
       gACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaaga                                                                                             <  2:1726244/100‑1 (MQ=255)
        aCGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCGGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTgaagaa                                                                                            >  1:876148/1‑100 (MQ=255)
         cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                                                                                           >  2:3495389/1‑100 (MQ=255)
         cGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAAt                                                                                           >  2:1373005/1‑100 (MQ=255)
          gCTCACCTGACGTTCATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                          <  2:1982235/100‑1 (MQ=255)
          gCTCACCTGACGTTAAGTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATa                                                                                          <  1:2420147/100‑1 (MQ=255)
           ntCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATaa                                                                                         >  1:4042819/2‑100 (MQ=255)
           cTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATat                                                                                         <  2:3753910/100‑2 (MQ=255)
            tCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatc                                                                                        <  1:1032100/100‑3 (MQ=255)
             cACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGt                                                                                       <  2:3672930/100‑1 (MQ=255)
               ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                     <  1:2052242/100‑6 (MQ=38)
               ccTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttc                                                                                     <  2:906825/100‑6 (MQ=38)
               ccTGACGTCAATTCATATTGATGATGGCTTAAGCGAGTTCTATCTGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCCATTTTTAAAGAATTAgcc                                                                                      >  1:1531662/1‑97 (MQ=255)
                cTGACTTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTCCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttct                                                                                    <  1:1235717/100‑7 (MQ=37)
                 tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctg                                                                                   <  1:3268366/100‑8 (MQ=37)
                 tGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctg                                                                                   <  2:4307651/100‑8 (MQ=37)
                  gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGAt                                                                                  <  1:13572/100‑1 (MQ=255)
                  gACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAACAATTGAAGAATatcttctgc                                                                                  <  1:3306280/100‑9 (MQ=37)
                   aCGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATa                                                                                 <  1:4321667/100‑1 (MQ=255)
                     gTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAAc                                                                               >  2:1054642/1‑100 (MQ=255)
                     gTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAAc                                                                               >  1:3352020/1‑100 (MQ=255)
                      ttAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACa                                                                              <  2:1078219/100‑1 (MQ=255)
                        aaTTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAg                                                                            <  1:3891305/100‑1 (MQ=255)
                           tCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTg                                                                         <  2:2561672/100‑1 (MQ=255)
                           tCATATTGATGACGGCTTAAGCGAGGTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATCGTTGAAGAATAAGTCGGATAACAAGCTg                                                                         <  1:1579149/100‑1 (MQ=255)
                                 tgatgaTGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                   >  1:1160554/1‑100 (MQ=255)
                                  gatgatGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATaaa                                                                   >  2:337114/1‑99 (MQ=255)
                                   atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACt                                                                 <  1:4197764/100‑1 (MQ=255)
                                   atgatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCCGTATAAAc                                                                  <  1:933342/99‑1 (MQ=255)
                                     gatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGGTATTCATCAAATGTTGAAGAATAAGCCGGATAAAAAACTGTATAAACTGa                                                               >  1:1394611/1‑100 (MQ=255)
                                     gatgGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGa                                                               <  1:2745128/100‑1 (MQ=255)
                                         gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                           >  2:4086224/1‑100 (MQ=255)
                                         gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                           >  1:999470/1‑100 (MQ=255)
                                         gCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGaa                                                           >  2:86194/1‑100 (MQ=255)
                                              aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                      >  2:3177165/1‑100 (MQ=255)
                                              aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                      >  2:3473514/1‑100 (MQ=255)
                                              aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                      >  2:2706600/1‑100 (MQ=255)
                                              aGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAata                                                      >  1:3160619/1‑100 (MQ=255)
                                                 gggTTGTATCCGGGTATCTACTTCTCTTTAACAGAATATATTCTTTAATTGTTGAAGAATATGTCGTATAACAAGCTGTATAAACTGACGAAAAATAtt                                                    <  1:2348279/97‑1 (MQ=255)
                                                 gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                   <  2:3501374/100‑1 (MQ=255)
                                                 gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                   >  1:3393236/1‑100 (MQ=255)
                                                 gAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTc                                                   <  2:2391165/100‑1 (MQ=255)
                                                  aGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCa                                                  <  2:1326750/100‑1 (MQ=255)
                                                      gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                              >  2:3584745/1‑100 (MQ=255)
                                                      gTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATgg                                              >  2:2678979/1‑100 (MQ=255)
                                                        aTCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGcc                                            >  2:2946651/1‑100 (MQ=255)
                                                         tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCg                                           >  2:4343047/1‑100 (MQ=255)
                                                         tCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGACCg                                           >  2:790205/1‑100 (MQ=255)
                                                          ccGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGa                                          >  2:707888/1‑100 (MQ=255)
                                                          ccGGGTATCTACTTCCCTGCAACAGAAGAGATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGTCGAAAAATATTCAATGGCCg                                           <  1:1602121/99‑1 (MQ=255)
                                                           cGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGaa                                         <  1:3592938/100‑1 (MQ=255)
                                                                     cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACt                                >  2:2173899/1‑99 (MQ=255)
                                                                     cTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTg                               >  2:1547176/1‑100 (MQ=255)
                                                                      ttCCCTTCAACAGAAGATAGTCTTCAATTGTTGAAGGATAAGTCGGATAATAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACaaaa                                  <  1:3314750/96‑1 (MQ=255)
                                                                       tCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                             >  1:397787/1‑100 (MQ=255)
                                                                        cccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                             <  2:3836100/99‑1 (MQ=255)
                                                                        cagtGCAACAGAAGATATCCTGCAATGGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                             <  2:4215035/96‑1 (MQ=255)
                                                                         ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCg                             <  1:3690450/98‑1 (MQ=255)
                                                                         ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                           <  1:1315711/100‑1 (MQ=255)
                                                                         ccTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTa                           >  2:209907/1‑100 (MQ=255)
                                                                            gCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTg                        <  1:1819411/100‑1 (MQ=255)
                                                                               atccgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                     >  1:4228329/9‑100 (MQ=37)
                                                                               atccgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAc                     >  1:4130807/9‑100 (MQ=37)
                                                                                cAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcg                    <  1:1171832/100‑1 (MQ=255)
                                                                                  gaagaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  >  1:2096953/1‑100 (MQ=255)
                                                                                  cgacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAAcgcg                  >  2:1511108/6‑100 (MQ=38)
                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 >  2:1640541/5‑100 (MQ=38)
                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 >  2:3212166/5‑100 (MQ=38)
                                                                                   gacttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCgg                 >  2:212823/5‑100 (MQ=38)
                                                                                    agaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                <  1:2706585/100‑1 (MQ=255)
                                                                                    acttatTCTTCAACAATTGAAGGCTAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGga                >  1:2567306/4‑100 (MQ=255)
                                                                                     gaTATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag               <  2:3699100/100‑1 (MQ=255)
                                                                                      ttatTCTTCAACAATTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgaga              >  2:250000/2‑100 (MQ=255)
                                                                                      atatTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGgag               <  1:499357/99‑1 (MQ=255)
                                                                                           cttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGc         <  2:1038848/100‑1 (MQ=255)
                                                                                            ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc        <  2:2699866/100‑1 (MQ=255)
                                                                                            ttcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGcc        <  2:22702/100‑1 (MQ=255)
                                                                                             tcAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCg       >  2:1921196/1‑100 (MQ=255)
                                                                                              cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg      >  1:2127554/1‑100 (MQ=255)
                                                                                              cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg      <  2:3604804/100‑1 (MQ=255)
                                                                                              cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg      <  2:3182665/100‑1 (MQ=255)
                                                                                              cAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCgg      <  2:3042297/100‑1 (MQ=255)
                                                                                                aTTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa    >  1:3030349/1‑100 (MQ=255)
                                                                                                aTTGTTGAAGAATAAGCCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa    >  2:606986/1‑100 (MQ=255)
                                                                                                 ttgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaaa   <  2:520045/100‑1 (MQ=255)
                                                                                                 ttgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaaa   >  1:2947366/1‑100 (MQ=255)
                                                                                                 ttgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGaa    >  2:2665298/1‑99 (MQ=255)
                                                                                                  tgttgAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAac  <  2:3062992/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
ACATAATGACGCTCACCTGACGTTAATTCATATTGATGATGGCTTAAGCGAGTTGTATCCGGGTATCTACTTCCCTGCAACAGAAGATATTCTTCAATTGTTGAAGAATAAGTCGGATAACAAGCTGTATAAACTGACGAAAAATATTCAATGGCCGAAGACAAAACTGCGTATTGAACGCGGAGAAATGCCGGAAAC  >  NC_000913/4113376‑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: ‑