Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,677,405 A→T 6.1% I319I (ATT→ATA)  yphE ← putative sugar ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,677,4050AT 6.1% 214.9 / 3.6 65I319I (ATT→ATA) yphEputative sugar ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base T (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CGATTGGTCAGCACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGATATCCTCCAGCTTGGGCTTATGGCGTAAC  >  NC_000913/2677265‑2677549
                                                                                                                                            |                                                                                                                                                
cGATTGGTCAGCACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGcc                                                                                                                                          <  3:282239/149‑1 (MQ=255)
 gATTGGTCAGCACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGccc                                                                                                                                         >  3:54896/1‑149 (MQ=255)
    tGGTCAGCACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGt                                                                                                                                      >  8:407822/1‑149 (MQ=255)
      gTCAGCACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGtt                                                                                                                                     >  6:326798/1‑148 (MQ=255)
       tCAGCACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCa                                                                                                                                   >  4:35232/1‑149 (MQ=255)
         aGCACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATa                                                                                                                                 <  5:333432/149‑1 (MQ=255)
           cACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATAct                                                                                                                               >  3:52912/1‑149 (MQ=255)
           cACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATAct                                                                                                                               >  1:252852/1‑149 (MQ=255)
              tGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACTATTTCGCCCTGTTCATActcct                                                                                                                            >  6:284752/1‑149 (MQ=255)
              tGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGAGCGTGATTTTCTCGCCGTTGATAACAATTTCGACCTGTTCATActcct                                                                                                                            >  7:118617/1‑149 (MQ=255)
              tGTATTTTCGTCAACGCCCAACCAGGGAATGATACCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACTATTTCGCCCTGTTCATActcct                                                                                                                            >  2:353523/1‑149 (MQ=255)
                tATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGc                                                                                                                                           <  1:114563/132‑1 (MQ=255)
                tATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGc                                                                                                                                           >  2:114563/1‑132 (MQ=255)
                tATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATActcctc                                                                                                                           <  4:52912/148‑1 (MQ=255)
                tATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATActcctc                                                                                                                           <  7:407822/148‑1 (MQ=255)
                tATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATActcctc                                                                                                                           <  8:616/148‑1 (MQ=255)
                     tCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGcc                                                                                                                      <  4:273000/148‑1 (MQ=255)
                           aCGCCCAACCAGCGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCACTGCCGCGTTTCAGCATGTCGACGTACTCGGGGCGCGTGATTTTCTCGCCGGTGATAACAATTTCGCCCGGTTCAGACTACTCCAGCGCACCAAt                                                                                                               >  6:21585/1‑149 (MQ=255)
                             gCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTc                                                                                                                                    <  8:363240/126‑1 (MQ=255)
                             gCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTc                                                                                                                                    >  7:363240/1‑126 (MQ=255)
                               ccAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCAcc                                                                                                           <  3:35232/149‑1 (MQ=255)
                                 aaCCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATActcct                                                                                                                            >  8:401303/1‑130 (MQ=255)
                                 aaCCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATActcct                                                                                                                            <  7:401303/130‑1 (MQ=255)
                                  aCCATGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTc                                                                                                        <  6:137757/149‑1 (MQ=255)
                                  aCCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTc                                                                                                        <  2:252852/149‑1 (MQ=255)
                                   ccAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCa                                                                                                       <  5:150677/149‑1 (MQ=255)
                                     aGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAg                                                                                                      >  7:26564/1‑148 (MQ=255)
                                       ggAATGATCCCATCTTCTTTTCGGTTTTCTGGCGTATATCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCATGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCaa                                                                                                   <  8:118617/149‑1 (MQ=255)
                                           tGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATCTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCa                                                                                               <  3:173191/149‑1 (MQ=255)
                                                cccGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTgc                                                                                           <  5:326798/148‑1 (MQ=255)
                                                 ccGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTgcgc                                                                                         >  3:96916/1‑149 (MQ=255)
                                                 ccGCGTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTgcgc                                                                                         >  1:87842/1‑149 (MQ=255)
                                                  cGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGcc                                                                                        >  8:124946/1‑149 (MQ=255)
                                                    cttcttTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACTATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCt                                                                                      <  1:353523/149‑1 (MQ=255)
                                                    cttcttTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACTATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCt                                                                                      <  5:284752/149‑1 (MQ=255)
                                                          tGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAAttt                                                                                                                                              <  7:212886/87‑1 (MQ=255)
                                                          tGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAAttt                                                                                                                                              >  8:212886/1‑87 (MQ=255)
                                                          tGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGcccc                                                                                 >  6:139299/1‑148 (MQ=255)
                                                               ttttCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCcag                                                                              >  1:376191/1‑146 (MQ=255)
                                                                      gCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGc                                                                     <  3:327245/148‑1 (MQ=255)
                                                                       cGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGAGCAGCaa                                                                   <  5:21585/149‑1 (MQ=255)
                                                                           taGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGcc                                                               <  1:216269/149‑1 (MQ=255)
                                                                              ccAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTgc                                                                                           <  1:74391/118‑1 (MQ=255)
                                                                              ccAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTgc                                                                                           >  2:74391/1‑118 (MQ=255)
                                                                              ccAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAg                                                            <  4:96916/149‑1 (MQ=255)
                                                                                aaTGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCa                                                                                                       <  5:47645/104‑1 (MQ=255)
                                                                                aaTGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCa                                                                                                       >  6:47645/1‑104 (MQ=255)
                                                                                 aTGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCAc                                                         <  4:54896/149‑1 (MQ=255)
                                                                                      gcgTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCGTACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGc                                                    >  5:223999/1‑149 (MQ=255)
                                                                                           tCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCAacga                                               >  2:197136/1‑149 (MQ=255)
                                                                                                   tCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGt                                       <  8:26564/149‑1 (MQ=255)
                                                                                                    cGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTAACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGt                                       >  7:263673/1‑148 (MQ=255)
                                                                                                      ccGTAATCGGGGCGCGTGATTTTCTCCCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTaaa                                    <  2:87842/149‑1 (MQ=255)
                                                                                                           aTCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGa                               <  5:139299/149‑1 (MQ=255)
                                                                                                           aTCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGa                               <  7:124946/149‑1 (MQ=255)
                                                                                                            tCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGa                               >  3:299323/1‑148 (MQ=255)
                                                                                                                    cgTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAAAGCCGAGCACTTCGCCACGACGAAGCGTAAAACTGATAtcatcc                       >  1:286148/1‑148 (MQ=255)
                                                                                                                                 ccGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGATATCCTCCAGCTTGGGCTTATg         <  4:299323/149‑1 (MQ=255)
                                                                                                                                  cGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGa                               >  2:45601/1‑126 (MQ=255)
                                                                                                                                  cGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGa                               <  1:45601/126‑1 (MQ=255)
                                                                                                                                       aTAACAATTTCGCCCTGTTCGTACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGATATCCTCCAGCTTGGGCTTATGGCGTa    <  6:223999/148‑1 (MQ=255)
                                                                                                                                       aTAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGATATCCTCCAGCTTGGGCTTATGGCGTaa   >  2:121767/1‑149 (MQ=255)
                                                                                                                                        tAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGATATCCTCCAGCTTGGGCTTATGGCGTAAc  >  2:399110/1‑149 (MQ=255)
                                                                                                                                         aaCAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGccc                                                                                       <  3:48901/63‑1 (MQ=255)
                                                                                                                                         aaCAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGccc                                                                                       >  4:48901/1‑63 (MQ=255)
                                                                                                                                            |                                                                                                                                                
CGATTGGTCAGCACTGTATTTTCGTCAACGCCCAACCAGGGAATGATCCCCGCTTCTTTGCGGTTTTCTGGCGTATAGCCAATGCCGCGTTTCAGCATGTCGCCGTAATCGGGGCGCGTGATTTTCTCGCCGTTGATAACAATTTCGCCCTGTTCATACTCCTCCAGCCCAACAATCGCCTTCAGCAATTCACTGCGCCCTGCCCCCAGCAGACCAGCAATGCCGAGCACTTCGCCACGACGTAGCGTAAAACTGATATCCTCCAGCTTGGGCTTATGGCGTAAC  >  NC_000913/2677265‑2677549

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑