Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 435,396 T→A 16.1% 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,3960TA16.1% 285.8 / 43.1 116I87N (ATT→AAT) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (54/43);  new base A (7/12);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 2.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

Only 100 of 118 total aligned reads displayed.

GGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCA  >  NC_000913/435297‑435495
                                                                                                   |                                                                                                   
gggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCAt                                                                                                     <  1:140277/100‑1 (MQ=255)
    ggCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcg                                                                                                 <  2:251333/100‑1 (MQ=255)
    ggCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcg                                                                                                 >  1:1785399/1‑100 (MQ=255)
     gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                                >  2:128674/1‑100 (MQ=255)
     gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                                <  2:1034323/100‑1 (MQ=255)
     gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                                <  2:326329/100‑1 (MQ=255)
      cGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCt                                                                                               >  1:2310875/1‑100 (MQ=255)
       gACTAATCCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTg                                                                                              >  2:1362650/1‑100 (MQ=255)
          tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTAc                                                                                           <  1:1062919/100‑1 (MQ=255)
          tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTAc                                                                                           >  1:1424283/1‑100 (MQ=255)
           aaTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACg                                                                                          >  1:1037578/1‑100 (MQ=255)
           aaTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACg                                                                                          <  1:768881/100‑1 (MQ=255)
           aaTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACg                                                                                          >  2:75473/1‑100 (MQ=255)
            aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActg                                                                                         <  1:626739/100‑3 (MQ=255)
            aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActg                                                                                         <  1:284930/100‑3 (MQ=255)
            aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActg                                                                                         <  2:542128/100‑3 (MQ=255)
             tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                        <  2:682455/100‑4 (MQ=255)
             tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                        <  2:353792/100‑4 (MQ=255)
               ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACt                                                                                      >  2:1401781/1‑100 (MQ=255)
               ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                                      <  2:1666523/100‑6 (MQ=38)
               ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                                      <  2:2434037/100‑6 (MQ=38)
               ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                                      <  1:192757/100‑6 (MQ=38)
               ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                                      <  1:1481152/100‑6 (MQ=38)
                cGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTg                                                                                     <  2:1599689/100‑1 (MQ=255)
                 gCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTg                                                                                     <  1:2034944/99‑1 (MQ=255)
                  cATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTTGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTc                                                                                   <  1:1652558/100‑1 (MQ=255)
                      ccTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                               >  1:197472/1‑100 (MQ=255)
                      ccTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                               >  2:1468398/1‑100 (MQ=255)
                      ccTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                               >  1:1012315/1‑100 (MQ=255)
                       cTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAAc                                                                              <  2:890110/100‑1 (MQ=255)
                       cTCGACGGACTGATGAAGCCATACCTGTCCCGCATGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAAc                                                                              >  2:832839/1‑100 (MQ=255)
                        tCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACg                                                                             >  1:2285176/1‑100 (MQ=255)
                        tCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACg                                                                             <  2:2072512/100‑1 (MQ=255)
                            cGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                         <  2:56696/100‑1 (MQ=255)
                             ggACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                         <  2:2425284/99‑1 (MQ=255)
                             ggACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCg                                                                        >  1:653005/1‑100 (MQ=255)
                              gACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGa                                                                       >  2:1121170/1‑100 (MQ=255)
                               aCTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAt                                                                      >  2:717776/1‑100 (MQ=255)
                                cTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtg                                                                     <  2:1817993/100‑1 (MQ=255)
                                cTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtg                                                                     >  1:2338657/1‑100 (MQ=255)
                                 tgatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgt                                                                    >  2:1923492/1‑100 (MQ=255)
                                  gatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgtg                                                                   >  2:239076/1‑100 (MQ=255)
                                    tgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGcc                                                                 >  2:196620/1‑100 (MQ=255)
                                    tgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGcc                                                                 >  2:777582/1‑100 (MQ=255)
                                     gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCa                                                                <  2:1783216/100‑1 (MQ=255)
                                      aaGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCAt                                                               >  1:2271369/1‑100 (MQ=255)
                                       aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATa                                                              <  1:2507535/100‑1 (MQ=255)
                                       aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATa                                                              <  1:162373/100‑1 (MQ=255)
                                           aTACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACaaa                                                          >  2:1815803/1‑100 (MQ=255)
                                               cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTgg                                                      <  1:1143438/100‑1 (MQ=255)
                                                tGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGc                                                     >  1:678968/1‑100 (MQ=255)
                                                 gTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                    <  2:2014214/100‑1 (MQ=255)
                                                  tcacGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACTAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCa                                                   <  1:2009824/97‑1 (MQ=255)
                                                   cccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCAt                                                  >  2:2394532/1‑100 (MQ=255)
                                                    ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCAtt                                                 >  1:1410035/1‑100 (MQ=255)
                                                      gCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTa                                                <  1:1955527/99‑1 (MQ=255)
                                                       ccTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAc                                              <  1:1362642/100‑1 (MQ=255)
                                                        ctgctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACg                                             <  2:1346579/100‑1 (MQ=255)
                                                         tgctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa                                            >  1:1528608/1‑100 (MQ=255)
                                                         tgctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa                                            <  1:122878/100‑1 (MQ=255)
                                                         tgctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa                                            <  2:497817/100‑1 (MQ=255)
                                                          gctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa                                           >  2:97525/1‑100 (MQ=255)
                                                          gctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa                                           <  2:1105961/100‑1 (MQ=255)
                                                           ctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                          >  2:2220053/1‑100 (MQ=255)
                                                           ctgGAAGAACTGAAAAAAGAAAAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                          <  1:413633/100‑1 (MQ=255)
                                                            tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGc                                         <  2:623619/100‑1 (MQ=255)
                                                               aagaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt                                      <  2:197474/100‑1 (MQ=255)
                                                                gcaggcgggACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                                     <  1:1335933/94‑1 (MQ=255)
                                                                agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                                     <  1:22019/100‑1 (MQ=255)
                                                                agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                                     <  1:2240015/100‑1 (MQ=255)
                                                                agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                                     <  1:1907975/100‑1 (MQ=255)
                                                                      gggACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTg                               >  1:120916/1‑100 (MQ=255)
                                                                       ggACCGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTCGCGCTGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTg                               >  1:745960/1‑99 (MQ=255)
                                                                        gACATGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGc                             <  2:1278461/100‑1 (MQ=255)
                                                                        gACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGc                             <  2:2230740/100‑1 (MQ=255)
                                                                          cAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGa                           >  1:2587643/1‑100 (MQ=255)
                                                                            ggTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                         >  2:1380811/1‑100 (MQ=255)
                                                                            ggTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                         >  2:2499162/1‑100 (MQ=255)
                                                                             gTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAAt                        >  1:599373/1‑100 (MQ=255)
                                                                               aGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCg                      >  2:2554506/1‑100 (MQ=255)
                                                                                gAAAAAGCAGTACTGCGCATAGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGt                     >  1:1576494/1‑100 (MQ=255)
                                                                                  aaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc                   <  2:926029/100‑1 (MQ=255)
                                                                                  aaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc                   >  2:58945/1‑100 (MQ=255)
                                                                                   cagttcgTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCg                  >  1:976013/7‑100 (MQ=38)
                                                                                    aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCgg                 <  1:1244329/100‑1 (MQ=255)
                                                                                      ttcgTACAGCGCAATGCGCAGTACGAGCTGTCTAAACGTAGCTATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgc                >  1:1157822/4‑99 (MQ=255)
                                                                                        aGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCTAAATCGTTCGGCGCa             >  2:1025771/1‑100 (MQ=255)
                                                                                         gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag            >  2:664555/1‑100 (MQ=255)
                                                                                         gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag            >  2:1600502/1‑100 (MQ=255)
                                                                                         gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag            >  2:593160/1‑100 (MQ=255)
                                                                                         gTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag            >  2:1557643/1‑100 (MQ=255)
                                                                                         gTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag            >  2:1558381/1‑100 (MQ=255)
                                                                                         gTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag            >  2:1661935/1‑100 (MQ=255)
                                                                                          tACTGCGCATTGCGCTGTACGCAATGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACAAAGCGATAGAACTGGCGAAATCGTTCGGCGCaga           >  1:1144073/1‑100 (MQ=255)
                                                                                           aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa          >  2:2019993/1‑100 (MQ=255)
                                                                                           aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa          >  2:1184798/1‑100 (MQ=255)
                                                                                                 cATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGc    <  1:2459533/100‑1 (MQ=255)
                                                                                                  aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGcc   >  1:2399094/1‑100 (MQ=255)
                                                                                                  aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGcc   <  1:1161911/100‑1 (MQ=255)
                                                                                                   ttGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCa  >  2:91775/1‑100 (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: