Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 434,687 A→G E109E (GAA→GAG nusB → transcription termination

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913434,6870AG84.6% 700.1 / 114.6 306E109E (GAA→GAGnusBtranscription termination
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base G (145/114);  total (167/139)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 308 total aligned reads displayed.

AGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACAAAAAGT  >  NC_000913/434596‑434778
                                                                                           |                                                                                           
aGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAc                                                                                            <  2:2980664/93‑1 (MQ=255)
 ggcaGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGTTCGAGCt                                                                                           >  2:2390433/4‑93 (MQ=255)
  gcaGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTg                                                                                          <  1:3506439/91‑1 (MQ=255)
   caGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTgg                                                                                         <  1:6283504/92‑1 (MQ=255)
    aGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGc                                                                                        >  1:6478612/1‑93 (MQ=255)
    aGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGc                                                                                        >  1:5561369/1‑93 (MQ=255)
    aGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGc                                                                                        <  1:2862905/93‑1 (MQ=255)
      aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGa                                                                                      >  1:1229308/1‑93 (MQ=255)
      aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGa                                                                                      >  2:5000306/1‑93 (MQ=255)
        aaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGaaa                                                                                    <  2:943876/93‑1 (MQ=255)
        aaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGaaa                                                                                    <  2:4787150/93‑1 (MQ=255)
        aaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGaaa                                                                                    <  2:1697968/93‑1 (MQ=255)
         aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAAt                                                                                   <  2:395930/93‑1 (MQ=255)
           gCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCg                                                                                 >  1:455144/1‑93 (MQ=255)
             aGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGtt                                                                               >  2:5818682/1‑93 (MQ=255)
             aGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGtt                                                                               <  2:2507464/93‑1 (MQ=255)
              gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTc                                                                              >  2:4235913/1‑93 (MQ=255)
              gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTc                                                                              <  1:3168441/93‑1 (MQ=255)
                aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCgg                                                                            <  1:2514924/93‑1 (MQ=255)
                aCTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCgg                                                                            >  1:1541863/1‑93 (MQ=255)
                 cTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgc                                                                           <  2:5362599/93‑1 (MQ=255)
                  tGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGgcg                                                                          >  2:3574644/1‑93 (MQ=255)
                  tGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGgcg                                                                          >  2:605129/1‑93 (MQ=255)
                  tGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgcg                                                                          >  2:3122746/1‑93 (MQ=255)
                   gcgcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGgcgc                                                                         >  1:3925419/1‑93 (MQ=255)
                     gcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCag                                                                       >  2:4148233/1‑93 (MQ=255)
                      cATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCaga                                                                      >  2:2840610/1‑93 (MQ=255)
                       aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCagaa                                                                     <  1:5138277/93‑1 (MQ=255)
                       aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCagaa                                                                     >  1:6343548/1‑93 (MQ=255)
                       aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCagaa                                                                     <  2:766498/93‑1 (MQ=255)
                       aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATCAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa                                                                     <  1:6225392/93‑1 (MQ=255)
                         tGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCagaaga                                                                   <  1:275891/93‑1 (MQ=255)
                          gcgcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGAc                                                                  <  2:1063379/93‑1 (MQ=255)
                          gcgcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGAc                                                                  >  1:5732386/1‑93 (MQ=255)
                           cgcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACa                                                                 <  1:3293425/93‑1 (MQ=255)
                           cgcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACa                                                                 <  1:4393861/93‑1 (MQ=255)
                           cgcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACa                                                                 >  2:4652550/1‑93 (MQ=255)
                           cgcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACa                                                                 <  2:3216717/93‑1 (MQ=255)
                              tGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGcc                                                              >  1:4730468/1‑93 (MQ=255)
                                tACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCAt                                                            >  1:6449683/1‑93 (MQ=255)
                                  cGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATaa                                                          <  1:1576102/93‑1 (MQ=255)
                                  cGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATaa                                                          >  1:7101509/1‑93 (MQ=255)
                                   gAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAg                                                         >  2:5985890/1‑93 (MQ=255)
                                   gAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAg                                                         <  2:6714690/93‑1 (MQ=255)
                                    aaCTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGt                                                        >  2:7131899/1‑93 (MQ=255)
                                    aaCTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGt                                                        >  1:38459/1‑93 (MQ=255)
                                     aCTGTCTAAACGTAGCGATGTGCCATACAAACTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGtt                                                       <  2:3298108/93‑1 (MQ=255)
                                      cTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTc                                                      <  2:3368044/93‑1 (MQ=255)
                                         tCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTc                                                   <  2:5887262/93‑1 (MQ=255)
                                         tCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTc                                                   >  1:2196622/1‑93 (MQ=255)
                                         tCTAAACGTAGCGATATGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTc                                                   >  1:1330804/1‑93 (MQ=255)
                                          cTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCa                                                  <  2:1807674/93‑1 (MQ=255)
                                           tAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCaa                                                 >  2:3772536/1‑93 (MQ=255)
                                           tAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCaa                                                 >  1:4305132/1‑93 (MQ=255)
                                           tAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCaa                                                 >  1:659281/1‑93 (MQ=255)
                                            aaaCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAAc                                                <  1:4860501/93‑1 (MQ=255)
                                            aaaCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAAc                                                <  1:6745298/93‑1 (MQ=255)
                                             aaCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                               >  2:563598/1‑93 (MQ=255)
                                              aCGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACgg                                              <  2:7404401/93‑1 (MQ=255)
                                                gTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCg                                            <  2:2950625/93‑1 (MQ=255)
                                                    cGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACt                                        <  1:2061035/93‑1 (MQ=255)
                                                     gATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTc                                       >  2:1293674/1‑93 (MQ=255)
                                                       tgtgCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATTGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGa                                     >  1:5460800/1‑93 (MQ=255)
                                                       tgtgCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGa                                     <  1:6649104/93‑1 (MQ=255)
                                                         tgCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATa                                   >  1:2343995/1‑93 (MQ=255)
                                                          gCCATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATaa                                  <  1:6103651/93‑1 (MQ=255)
                                                            cATACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAag                                >  2:2262068/1‑93 (MQ=255)
                                                             aTACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAaga                               >  1:5931364/1‑92 (MQ=255)
                                                              tACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                              >  2:7293229/1‑93 (MQ=255)
                                                              tACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                              <  1:1017195/93‑1 (MQ=255)
                                                              tACAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                              >  2:1404232/1‑93 (MQ=255)
                                                               aCAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcag                             >  2:7038520/1‑93 (MQ=255)
                                                               aCAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcag                             <  1:2296164/93‑1 (MQ=255)
                                                                cAAAGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcagc                            <  2:44745/93‑1 (MQ=255)
                                                                cAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcagc                            >  1:1696410/1‑93 (MQ=255)
                                                                  aaGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCAc                          <  1:4125148/93‑1 (MQ=255)
                                                                   aGTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCAcc                         >  1:2131732/1‑93 (MQ=255)
                                                                    gTGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCt                        <  2:6267832/93‑1 (MQ=255)
                                                                     tGGCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCtg                       <  1:2558572/93‑1 (MQ=255)
                                                                      ggCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCtgt                      <  2:7343053/93‑1 (MQ=255)
                                                                      ggCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCtgt                      <  2:520887/93‑1 (MQ=255)
                                                                       gCCATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCtgtg                     <  1:506416/93‑1 (MQ=255)
                                                                        ccATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGa                    <  2:4809518/93‑1 (MQ=255)
                                                                        ccATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGa                    >  1:3464908/1‑93 (MQ=255)
                                                                        ccATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGa                    >  2:2045483/1‑93 (MQ=255)
                                                                        ccATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGa                    <  1:1887017/93‑1 (MQ=255)
                                                                         cATTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAt                   >  2:5038712/1‑93 (MQ=255)
                                                                          aTTAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAtt                  >  2:1137282/1‑93 (MQ=255)
                                                                           ttAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTc                 >  2:6349671/1‑93 (MQ=255)
                                                                           ttAACGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTc                 >  2:570397/1‑93 (MQ=255)
                                                                               cGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc             <  1:784734/93‑1 (MQ=255)
                                                                               cGAAGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc             <  1:2384068/93‑1 (MQ=255)
                                                                               cGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc             >  2:1959700/1‑93 (MQ=255)
                                                                                 aaGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTa           <  1:6956937/93‑1 (MQ=255)
                                                                                  aGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTaa          <  1:4456077/93‑1 (MQ=255)
                                                                                  aGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTaa          >  2:1331464/1‑93 (MQ=255)
                                                                                  aGCGATCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTaa          <  1:2978106/93‑1 (MQ=255)
                                                                                   gCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAAc         >  1:4906632/1‑93 (MQ=255)
                                                                                      aTCGAGCTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACaaa      >  2:4587983/1‑93 (MQ=255)
                                                                                        cGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACaaaaa    <  1:7139767/93‑1 (MQ=255)
                                                                                           |                                                                                           
AGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACAAAAAGT  >  NC_000913/434596‑434778

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