Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 435,396 T→A 15.8% I87N (ATT→AAT)  nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913435,3960TA15.8% 372.2 / 58.3 145I87N (ATT→AAT) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (69/53);  new base A (5/18);  total (74/71)
Fisher's exact test for biased strand distribution p-value = 2.74e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 148 total aligned reads displayed.

GGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCA  >  NC_000913/435297‑435495
                                                                                                   |                                                                                                   
gggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCAt                                                                                                     >  2:160586/1‑100 (MQ=255)
    ggCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcg                                                                                                 >  2:110906/1‑100 (MQ=255)
    ggCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATgcg                                                                                                 <  2:375221/100‑1 (MQ=255)
     gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                                <  1:2367770/100‑1 (MQ=255)
     gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                                >  2:2337349/1‑100 (MQ=255)
      cGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCt                                                                                               >  1:1456838/1‑100 (MQ=255)
        aCTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGt                                                                                             <  2:2003047/100‑1 (MQ=255)
         cTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTa                                                                                            <  2:456689/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)
            aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGa                                                                                         >  1:848750/1‑100 (MQ=255)
            aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGa                                                                                         >  1:2356378/1‑100 (MQ=255)
            aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActg                                                                                         <  2:2551547/100‑3 (MQ=255)
             tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGaa                                                                                        >  2:849261/1‑100 (MQ=255)
             tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                        <  2:1378007/100‑4 (MQ=255)
             tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                        <  2:1678706/100‑4 (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:1517853/100‑6 (MQ=38)
                 gCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctttt                                                                                    <  2:2573823/100‑8 (MQ=38)
                    tACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctttttct                                                                                 <  2:1384641/100‑11 (MQ=37)
                    tACCTCGACGGACCGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctttttc                                                                                  <  1:2264984/99‑10 (MQ=37)
                      ccTCGACGGACTCATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                               <  1:163652/100‑1 (MQ=255)
                       cTCGACGGACTGATGAAGCCATACGTATCCCGCCTGCTGGAAGAACTGAGACAGGAAGAAAAAGCAGTACTGCGCATTGCGCTGTACGCACTGTCTaaa                                                                               >  1:2835954/1‑99 (MQ=255)
                         cGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGt                                                                            >  1:2545027/1‑100 (MQ=255)
                          gACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTa                                                                           <  2:1794446/100‑1 (MQ=255)
                           aCGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAg                                                                          <  1:1177994/100‑1 (MQ=255)
                            cGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                         <  2:748939/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                                                                   >  2:2571476/1‑100 (MQ=255)
                                  gatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgtg                                                                   >  1:1586621/1‑100 (MQ=255)
                                  gatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgtg                                                                   <  1:1305808/100‑1 (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)
                                      aaGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCAt                                                               <  2:2805494/100‑1 (MQ=255)
                                       aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCAt                                                               >  1:359432/1‑99 (MQ=255)
                                       aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATa                                                              <  2:2128782/100‑1 (MQ=255)
                                        gCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATAc                                                             <  2:1394814/100‑1 (MQ=255)
                                          cATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACaa                                                           >  2:2630589/1‑100 (MQ=255)
                                          cATACCTGTCCCGCATGCTGGAAAAACTGGGACAGGTAGAAAAAGCAGTACGGCGTATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACa                                                            <  1:1312425/99‑1 (MQ=255)
                                              ccTGTCCCTCCTGCTGGGAGAACTGGCACAGGTAGAAAAAGCAGTACTGCGCAATGCGCTGTACGAACTGTTTGAAAGTAGCTCTGTGCCATAAAAAGTg                                                       >  1:2305002/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:705128/100‑1 (MQ=255)
                                                    ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCAtt                                                 >  2:2522649/1‑100 (MQ=255)
                                                      gCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTaa                                               >  2:1010470/1‑100 (MQ=255)
                                                      gCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTaa                                               >  2:1010410/1‑100 (MQ=255)
                                                      gCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTaa                                               >  1:161793/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)
                                                           ctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                          >  1:2794632/1‑100 (MQ=255)
                                                            tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGc                                         <  2:2909286/100‑1 (MQ=255)
                                                            tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGc                                         <  1:58118/100‑1 (MQ=255)
                                                             ggAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCg                                        >  1:1006538/1‑100 (MQ=255)
                                                             ggAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCg                                        >  2:1744927/1‑100 (MQ=255)
                                                               aagaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt                                      >  1:1237999/1‑100 (MQ=255)
                                                                agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                                     >  2:51577/1‑100 (MQ=255)
                                                                 gaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCg                                    <  2:171545/100‑1 (MQ=255)
                                                                  aaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGa                                   <  2:565604/100‑1 (MQ=255)
                                                                   aCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGaa                                  >  2:70910/1‑100 (MQ=255)
                                                                    cTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAc                                 <  2:2687657/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                         >  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:2146076/1‑100 (MQ=255)
                                                                                  aaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc                   >  2:2103661/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:415187/100‑1 (MQ=255)
                                                                                    aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCgg                 <  1:1046468/100‑1 (MQ=255)
                                                                                    aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCgg                 <  1:17352/100‑1 (MQ=255)
                                                                                     aGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgc                <  2:245836/100‑1 (MQ=255)
                                                                                           aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa          <  1:2343356/100‑1 (MQ=255)
                                                                                           aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa          <  2:1877277/100‑1 (MQ=255)
                                                                                           aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa          <  2:1830434/100‑1 (MQ=255)
                                                                                               cgcATTGCGTTGTACGAACTGTCTATGCGCGGCGATGTGCCATACAAAGTGGCTATTAACGAAGGGATGGAACTTGCGAAATCGTTAGGCGCagaag         <  1:746188/97‑1 (MQ=255)
                                                                                                   ttGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCa  <  2:433334/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCA  >  NC_000913/435297‑435495

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