Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 435,390 T→A 16.7% 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,3900TA16.7% 288.8 / 46.2 114L85Q (CTG→CAG) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (52/43);  new base A (7/12);  total (59/55)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 115 total aligned reads displayed.

TGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAA  >  NC_000913/435291‑435487
                                                                                                   |                                                                                                 
   ccGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTgcg                                                                                                <  1:1098829/100‑1 (MQ=255)
    cGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTgcg                                                                                                <  1:2556150/99‑1 (MQ=255)
      gggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCAt                                                                                             <  1:140277/100‑1 (MQ=255)
          ggCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcg                                                                                         <  2:251333/100‑1 (MQ=255)
          ggCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcg                                                                                         >  1:1785399/1‑100 (MQ=255)
           gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                        <  2:326329/100‑1 (MQ=255)
           gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                        <  2:1034323/100‑1 (MQ=255)
           gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                        >  2:128674/1‑100 (MQ=255)
            cGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCt                                                                                       >  2:2431783/1‑100 (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)
                  aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActg                                                                                 <  1:626739/100‑3 (MQ=255)
                  aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActg                                                                                 <  2:542128/100‑3 (MQ=255)
                  aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActg                                                                                 <  1:284930/100‑3 (MQ=255)
                   tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                <  2:682455/100‑4 (MQ=255)
                   tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                <  2:353792/100‑4 (MQ=255)
                   tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                <  2:2132366/100‑4 (MQ=255)
                     ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACt                                                                              >  2:1401781/1‑100 (MQ=255)
                     ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                              <  1:192757/100‑6 (MQ=38)
                     ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                              <  2:2434037/100‑6 (MQ=38)
                     ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                              <  1:1481152/100‑6 (MQ=38)
                     ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgctt                                                                              <  2:1666523/100‑6 (MQ=38)
                      cGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTg                                                                             >  1:1940156/1‑100 (MQ=255)
                      cGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTg                                                                             <  2:1599689/100‑1 (MQ=255)
                      cGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgcttt                                                                             <  2:414955/100‑7 (MQ=38)
                            ccTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                       >  2:1468398/1‑100 (MQ=255)
                            ccTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                       >  1:1012315/1‑100 (MQ=255)
                            ccTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                       >  1:197472/1‑100 (MQ=255)
                            ccTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                       >  1:1370998/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)
                                   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:777582/1‑100 (MQ=255)
                                          tgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGcc                                                         >  2:196620/1‑100 (MQ=255)
                                           gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCa                                                        <  2:1783216/100‑1 (MQ=255)
                                            aaGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCAt                                                       >  1:431834/1‑100 (MQ=255)
                                            aaGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCAt                                                       >  1:2271369/1‑100 (MQ=255)
                                             aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATa                                                      <  1:162373/100‑1 (MQ=255)
                                             aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATa                                                      <  1:2507535/100‑1 (MQ=255)
                                                 aTACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACaaa                                                  >  2:1815803/1‑100 (MQ=255)
                                                    ccTGTCCCGCCTGCTGTAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGt                                                <  1:2217841/99‑1 (MQ=255)
                                                     cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTgg                                              <  1:1143438/100‑1 (MQ=255)
                                                      tGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGc                                             >  1:678968/1‑100 (MQ=255)
                                                       gTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                            <  2:2014214/100‑1 (MQ=255)
                                                       gTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                            <  2:1386957/100‑1 (MQ=255)
                                                        tcacGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACTAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCa                                           <  1:2009824/97‑1 (MQ=255)
                                                         cccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCAt                                          >  2:2394532/1‑100 (MQ=255)
                                                            gCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTa                                        <  1:1955527/99‑1 (MQ=255)
                                                             ccTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAc                                      <  1:1362642/100‑1 (MQ=255)
                                                               tgctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa                                    <  2:497817/100‑1 (MQ=255)
                                                               tgctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa                                    >  1:1528608/1‑100 (MQ=255)
                                                               tgctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa                                    <  1:122878/100‑1 (MQ=255)
                                                                gctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa                                   >  2:97525/1‑100 (MQ=255)
                                                                gctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa                                   <  2:1105961/100‑1 (MQ=255)
                                                                 ctgGAAGAACTGAAAAAAGAAAAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                  <  1:413633/100‑1 (MQ=255)
                                                                  tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGc                                 <  2:623619/100‑1 (MQ=255)
                                                                  tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGc                                 <  2:623049/100‑1 (MQ=255)
                                                                     aagaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt                              <  2:197474/100‑1 (MQ=255)
                                                                      gcaggcgggACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                             <  1:1335933/94‑1 (MQ=255)
                                                                      agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                             <  1:1907975/100‑1 (MQ=255)
                                                                      agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                             <  1:22019/100‑1 (MQ=255)
                                                                      agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                             <  1:2240015/100‑1 (MQ=255)
                                                                            gggACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTg                       >  1:120916/1‑100 (MQ=255)
                                                                            gggACAGGTAGAAAAAGCAGTACTGCGCATTCCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTg                       >  1:787921/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)
                                                                                  ggTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTTGCGaaa                 >  1:15488/1‑100 (MQ=255)
                                                                                  ggTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                 >  2:1380811/1‑100 (MQ=255)
                                                                                     aGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCg              >  2:2554506/1‑100 (MQ=255)
                                                                                      gAAAAAGCAGTACTGCGCATAGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGt             >  1:1576494/1‑100 (MQ=255)
                                                                                        aaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc           <  2:808927/100‑1 (MQ=255)
                                                                                        aaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc           >  2:58945/1‑100 (MQ=255)
                                                                                        aaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc           <  2:926029/100‑1 (MQ=255)
                                                                                         cagttcgTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCg          >  1:976013/7‑100 (MQ=38)
                                                                                          aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCgg         <  1:1244329/100‑1 (MQ=255)
                                                                                            ttcgTACAGCGCAATGCGCAGTACGAGCTGTCTAAACGTAGCTATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgc        >  1:1157822/4‑99 (MQ=255)
                                                                                              aGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCa     >  2:114516/1‑100 (MQ=255)
                                                                                               gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag    >  2:1600502/1‑100 (MQ=255)
                                                                                               gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag    >  2:664555/1‑100 (MQ=255)
                                                                                               gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag    >  2:593160/1‑100 (MQ=255)
                                                                                               gTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag    >  2:1557643/1‑100 (MQ=255)
                                                                                               gTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCag    >  2:1661935/1‑100 (MQ=255)
                                                                                                tACTGCGCATTGCGCTGTACGCAATGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACAAAGCGATAGAACTGGCGAAATCGTTCGGCGCaga   >  1:1144073/1‑100 (MQ=255)
                                                                                                 aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa  >  2:1184798/1‑100 (MQ=255)
                                                                                                 aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa  >  2:2019993/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
TGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAA  >  NC_000913/435291‑435487

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: