Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 435,390 T→A 15.0% L85Q (CTG→CAG)  nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913435,3900TA15.0% 389.3 / 57.6 147L85Q (CTG→CAG) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (72/52);  new base A (4/18);  total (76/71)
Fisher's exact test for biased strand distribution p-value = 8.43e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

Only 100 of 147 total aligned reads displayed.

GGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAA  >  NC_000913/435292‑435487
                                                                                                  |                                                                                                 
ggCCGGGGTGGGGACTAATGCCGCATGCCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTg                                                                                                  <  1:2238509/100‑1 (MQ=255)
  ccGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTgcg                                                                                                <  2:2139895/100‑1 (MQ=255)
   cGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTgcgc                                                                                               <  1:2606021/100‑1 (MQ=255)
     gggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCAt                                                                                             >  2:160586/1‑100 (MQ=255)
          gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                        <  1:2367770/100‑1 (MQ=255)
          gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                        >  2:2337349/1‑100 (MQ=255)
           cGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCt                                                                                       >  1:1456838/1‑100 (MQ=255)
           cGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATgcgc                                                                                        <  2:2718958/99‑1 (MQ=255)
             aCTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGt                                                                                     <  2:2003047/100‑1 (MQ=255)
              cTAATACCGCAAGTCGCGAGGAACTGATGAAGCGATACCTGTCTCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATgcgcag                                                                                      <  2:1173214/98‑3 (MQ=255)
                aaTACCGCATACCTCGGCGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTAct                                                                                  <  1:1700917/100‑2 (MQ=255)
                aaTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACg                                                                                  >  1:565603/1‑100 (MQ=255)
                aaTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTAct                                                                                  <  1:1443073/100‑2 (MQ=255)
                 aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActg                                                                                 <  2:2551547/100‑3 (MQ=255)
                  tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGaa                                                                                >  2:849261/1‑100 (MQ=255)
                  tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                <  2:1678706/100‑4 (MQ=255)
                  tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                <  2:1378007/100‑4 (MQ=255)
                   aCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAAc                                                                               <  2:2070687/100‑1 (MQ=255)
                   aCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAAc                                                                               >  2:415189/1‑100 (MQ=255)
                   aCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                <  2:2718229/99‑4 (MQ=255)
                    ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACt                                                                              >  1:638854/1‑100 (MQ=255)
                    ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCACTGTACGAACt                                                                              >  2:1604515/1‑100 (MQ=255)
                    ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                              <  2:2512444/100‑6 (MQ=38)
                    ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                              <  2:567780/100‑6 (MQ=38)
                    ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                              <  2:1517853/100‑6 (MQ=38)
                       cATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTc                                                                           <  1:581747/100‑1 (MQ=255)
                         tACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctttttct                                                                         <  2:1384641/100‑11 (MQ=37)
                         tACCTCGACGGACCGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctttttc                                                                          <  1:2264984/99‑10 (MQ=37)
                          aCCGCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAACAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaa                                                                        <  2:1345683/100‑1 (MQ=255)
                           ccTCGACGGACTCATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                       <  1:163652/100‑1 (MQ=255)
                            cTCGACGGACTGATGAAGCCATACGTATCCCGCCTGCTGGAAGAACTGAGACAGGAAGAAAAAGCAGTACTGCGCATTGCGCTGTACGCACTGTCTaaa                                                                       >  1:2835954/1‑99 (MQ=255)
                               gACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTa                                                                   <  2:1794446/100‑1 (MQ=255)
                                aCGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAg                                                                  <  1:1177994/100‑1 (MQ=255)
                                 cGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                 <  2:748939/100‑1 (MQ=255)
                                 cGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                 <  1:554651/100‑1 (MQ=255)
                                   gACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGa                                                               <  1:234124/100‑1 (MQ=255)
                                    aCTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAt                                                              <  2:2047341/100‑1 (MQ=255)
                                     cTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACCGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtg                                                             <  1:1858728/100‑1 (MQ=255)
                                     cTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtg                                                             <  1:1040089/100‑1 (MQ=255)
                                       gatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgtg                                                           >  1:1586621/1‑100 (MQ=255)
                                       gatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgtg                                                           <  1:1305808/100‑1 (MQ=255)
                                       gatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgtg                                                           >  2:2571476/1‑100 (MQ=255)
                                       gatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgtg                                                           >  2:1398536/1‑100 (MQ=255)
                                         tgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgtg                                                           >  2:1666437/1‑98 (MQ=255)
                                         tgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGcc                                                         >  1:854747/1‑100 (MQ=255)
                                         tgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGcc                                                         >  2:1761321/1‑100 (MQ=255)
                                         tgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGcc                                                         >  1:1929986/1‑100 (MQ=255)
                                            aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCAt                                                       >  1:359432/1‑99 (MQ=255)
                                            aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATa                                                      <  2:2128782/100‑1 (MQ=255)
                                             gCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATAc                                                     <  2:1394814/100‑1 (MQ=255)
                                               cATACCTGTCCCGCATGCTGGAAAAACTGGGACAGGTAGAAAAAGCAGTACGGCGTATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACa                                                    <  1:1312425/99‑1 (MQ=255)
                                                   ccTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTg                                               >  2:496656/1‑100 (MQ=255)
                                                   ccTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTg                                               <  1:2774913/100‑1 (MQ=255)
                                                   ccTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTg                                               >  1:1517004/1‑100 (MQ=255)
                                                     tGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGc                                             <  2:1573285/100‑1 (MQ=255)
                                                      gTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                            >  1:1608121/1‑100 (MQ=255)
                                                      gTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                            >  2:1760306/1‑100 (MQ=255)
                                                      gTCCCGCCTGCTGGAAGAACTGGGAAAGGAAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                            <  1:291388/100‑1 (MQ=255)
                                                        cccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCAt                                          >  2:58119/1‑100 (MQ=255)
                                                        cccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCAt                                          >  2:2658261/1‑100 (MQ=255)
                                                        cccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCAt                                          >  2:1049795/1‑100 (MQ=255)
                                                        cccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCAt                                          <  2:705128/100‑1 (MQ=255)
                                                         ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCAtt                                         >  2:2522649/1‑100 (MQ=255)
                                                           gCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTaa                                       >  2:1010470/1‑100 (MQ=255)
                                                              tgctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa                                    <  2:680313/100‑1 (MQ=255)
                                                              tgctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa                                    <  1:255038/100‑1 (MQ=255)
                                                                ctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                  >  1:1155147/1‑100 (MQ=255)
                                                                 tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGc                                 <  2:2909286/100‑1 (MQ=255)
                                                                  ggAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCg                                >  2:1744927/1‑100 (MQ=255)
                                                                  ggAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCg                                >  1:1006538/1‑100 (MQ=255)
                                                                    aagaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt                              >  1:1237999/1‑100 (MQ=255)
                                                                      gaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCg                            <  2:171545/100‑1 (MQ=255)
                                                                       aaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGa                           <  2:565604/100‑1 (MQ=255)
                                                                        aCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGaa                          >  1:2748419/1‑100 (MQ=255)
                                                                        aCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGaa                          >  2:70910/1‑100 (MQ=255)
                                                                         cTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAc                         <  2:2687657/100‑1 (MQ=255)
                                                                          tGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACt                        <  1:1668211/100‑1 (MQ=255)
                                                                           gggACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTg                       >  2:2016578/1‑100 (MQ=255)
                                                                               cAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGa                   >  2:2066455/1‑100 (MQ=255)
                                                                               cAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGa                   >  2:2597721/1‑100 (MQ=255)
                                                                               cAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGAGATCggaagagcac                   >  2:1861855/1‑91 (MQ=255)
                                                                                 ggTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                 >  1:2440005/1‑100 (MQ=255)
                                                                                 ggTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                 >  2:473216/1‑100 (MQ=255)
                                                                                  gTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                 >  1:153712/1‑99 (MQ=255)
                                                                                  gTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAAt                <  2:638855/100‑1 (MQ=255)
                                                                                   tAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATc               >  1:1982565/1‑100 (MQ=255)
                                                                                     agacagttcgTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGt             >  2:2865936/10‑100 (MQ=37)
                                                                                      aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGtt            >  1:999128/1‑100 (MQ=255)
                                                                                      aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGtt            >  1:2146076/1‑100 (MQ=255)
                                                                                       aaaaGCAATACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc           >  2:2587469/1‑100 (MQ=255)
                                                                                        cagttcgTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCg          >  1:1204790/7‑100 (MQ=38)
                                                                                        cagttcgTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCg          >  2:696497/7‑100 (MQ=38)
                                                                                        aaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCg          <  2:150799/100‑1 (MQ=255)
                                                                                         aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCgg         <  1:1046468/100‑1 (MQ=255)
                                                                                          aGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgc        <  2:2878578/100‑1 (MQ=255)
                                                                                          aGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgc        <  2:245836/100‑1 (MQ=255)
                                                                                              gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag    >  2:433869/1‑100 (MQ=255)
                                                                                                aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa  <  2:1877277/100‑1 (MQ=255)
                                                                                                aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa  <  2:1830434/100‑1 (MQ=255)
                                                                                                aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa  <  1:2343356/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
GGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAA  >  NC_000913/435292‑435487

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