Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 946,381 G→A 6.5% G171R (GGG→AGG)  ycaD → putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913946,3810GA 6.5% 192.1 / 2.5 62G171R (GGG→AGG) ycaDputative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (29/28);  new base A (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-01

CAGTGGGACGTCACGTAACCGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTCTGGGATCTCTATA  >  NC_000913/946233‑946529
                                                                                                                                                    |                                                                                                                                                    
cAGTGGGACGTCACGTAACCGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAg                                                                                                                                                      <  7:251179/149‑1 (MQ=255)
   tGGGACGTCACGTAACCGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGa                                                                                                                                                   >  8:302907/1‑149 (MQ=255)
     ggACGTCACGTAACCGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATc                                                                                                                                                 >  5:43459/1‑149 (MQ=255)
     ggACGTCACGTAACCGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGGTACAGGTTTGACGTTGGCAGGGATc                                                                                                                                                 >  6:144546/1‑149 (MQ=255)
     ggACGTCACGTAACCGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACCTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGAt                                                                                                                                                  <  3:208231/148‑1 (MQ=255)
      gACGTCACGTAACCGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGa                                                                                                                                                   <  7:133990/146‑1 (MQ=255)
      gACGTCACGTAACCGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGa                                                                                                                                                   >  8:133990/1‑146 (MQ=255)
                   cGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGATACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTgttgtt                                                                                                                                   >  2:169827/1‑149 (MQ=255)
                       ggCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTAc                                                                                                                                <  4:134982/148‑1 (MQ=255)
                        gCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGTAAAGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTAcgc                                                                                                                              <  4:221881/149‑1 (MQ=255)
                              gcttgctGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGtgtg                                                                                                                         <  6:43459/148‑1 (MQ=255)
                                   ctgcGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGaa                                                                                                                    >  6:108805/1‑148 (MQ=255)
                                    tgcGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGGTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAAt                                                                                                                   >  5:216767/1‑148 (MQ=255)
                                      cGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATc                                                                                                                  <  3:65489/147‑1 (MQ=255)
                                        tataTGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATcagca                                                                                                              <  8:33954/149‑1 (MQ=255)
                                        gatatGATGCTGTCTTACGTGCGCACCTTTTTAGGCCAGTGACGGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTTCCGTGGGTTACAGGTTTGACTTTGGCAGGGATCTTACCGAGGTTGTTTACGCGTGTGCTGAATcagca                                                                                                              <  2:107029/148‑1 (MQ=255)
                                         atatGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATcagcag                                                                                                             >  2:329436/1‑149 (MQ=255)
                                          tatGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAgg                                                                                                            <  5:108805/149‑1 (MQ=255)
                                              atgGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCTTGTGCTGAATCAGCAGGCTg                                                                                                         <  7:8293/148‑1 (MQ=255)
                                                        aCGTGGGAACGTTTTTATGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTATTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGa                                                                                               <  6:149219/148‑1 (MQ=255)
                                                          gTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTc                                                                                            >  1:85197/1‑149 (MQ=255)
                                                                          gCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTc                                                                             >  7:76590/1‑148 (MQ=255)
                                                                                 aCTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTa                                                                      >  5:304989/1‑148 (MQ=255)
                                                                                  cTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTaa                                                                     >  6:318613/1‑148 (MQ=255)
                                                                                   tGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTaaa                                                                    >  2:295995/1‑148 (MQ=255)
                                                                                     gTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGTCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACt                                                                 >  6:79968/1‑149 (MQ=255)
                                                                                     gTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTaaaa                                                                   <  1:329436/147‑1 (MQ=255)
                                                                                     gTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTaaaa                                                                   <  6:214801/147‑1 (MQ=255)
                                                                                     gTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACt                                                                 >  3:298943/1‑149 (MQ=255)
                                                                                        aGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTcc                                                               <  1:295995/148‑1 (MQ=255)
                                                                                          cAAAGTTTCAACCGAGCTTATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTc                                                            <  6:226045/149‑1 (MQ=255)
                                                                                          cAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGt                                                             >  6:91073/1‑148 (MQ=255)
                                                                                          cAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAAGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGt                                                             >  2:325981/1‑148 (MQ=255)
                                                                                          cAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAAGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGt                                                             >  2:321091/1‑148 (MQ=255)
                                                                                                 tCATCCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGTCGTTGGCAGGGATCTTACCGCTGTTGTTTACGCTTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGAAGTCTATTACGTCAATGCTAAAACTCCGTCAGgcgc                                                      <  8:270724/148‑1 (MQ=255)
                                                                                                 tCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGgcgc                                                      <  5:318613/148‑1 (MQ=255)
                                                                                                   aaCCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGc                                                   >  6:233117/1‑149 (MQ=255)
                                                                                                          cTGATGTCCGTATTGCCGTGGGATACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGt                                            >  6:223031/1‑149 (MQ=255)
                                                                                                             aTGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAAGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGaa                                         <  1:325981/149‑1 (MQ=255)
                                                                                                             aTGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAAGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGaa                                         <  1:321091/149‑1 (MQ=255)
                                                                                                              tGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGaa                                         <  8:76590/148‑1 (MQ=255)
                                                                                                              tGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAAt                                        <  5:79968/149‑1 (MQ=255)
                                                                                                               gTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATg                                       >  6:88213/1‑149 (MQ=255)
                                                                                                                 ccGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGc                                     >  6:285952/1‑149 (MQ=255)
                                                                                                                  cGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCt                                    <  6:304989/149‑1 (MQ=255)
                                                                                                                         ccGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTAt                             >  7:303757/1‑149 (MQ=255)
                                                                                                                          cGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTAtc                            >  2:326947/1‑149 (MQ=255)
                                                                                                                                ttACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGa                      <  7:302907/149‑1 (MQ=255)
                                                                                                                                 tACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGaa                     <  6:216767/149‑1 (MQ=255)
                                                                                                                                  aCAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCGGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAAt                    >  4:241051/1‑149 (MQ=255)
                                                                                                                                  aCAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGaa                     <  5:91073/148‑1 (MQ=255)
                                                                                                                                      gTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGtt                >  6:48139/1‑149 (MQ=255)
                                                                                                                                       tttGACGTTGGCATGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTc               <  1:230726/149‑1 (MQ=255)
                                                                                                                                       tttGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTc               <  1:169827/149‑1 (MQ=255)
                                                                                                                                         tGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTCTg             >  2:127456/1‑149 (MQ=255)
                                                                                                                                             gTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTCTGGGAt         >  2:252624/1‑149 (MQ=255)
                                                                                                                                              ttGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGt                                            <  4:49298/113‑1 (MQ=255)
                                                                                                                                              ttGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGt                                            >  3:49298/1‑113 (MQ=255)
                                                                                                                                              ttGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTCTGGGAtc        <  1:252624/149‑1 (MQ=255)
                                                                                                                                              ttGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTCTGGGAt         >  1:143848/1‑148 (MQ=255)
                                                                                                                                                 gCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATcagc                                                                                                               <  7:129966/43‑1 (MQ=255)
                                                                                                                                                 gCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATcagc                                                                                                               >  8:129966/1‑43 (MQ=255)
                                                                                                                                                 gCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGTAATCGTTCTGGGATctct     <  5:88213/149‑1 (MQ=255)
                                                                                                                                                    gggATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTCTGGGATCTCtata  <  8:303757/149‑1 (MQ=255)
                                                                                                                                                    gggAGCTTACCGCTGCTTTCTACGCGTGTGCTGAATCAGCAGGCGGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTCTGGGATCTCtata  <  6:92767/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAGTGGGACGTCACGTAACCGTGGGCGTTTGCTTGCTGCGTATATGATGGTTTATTACGTGGGAACGTTTTTAGGCCAGTTACTGGTCAGCAAAGTTTCAACCGAGCTGATGTCCGTATTGCCGTGGGTTACAGGTTTGACGTTGGCAGGGATCTTACCGCTGTTGTTTACGCGTGTGCTGAATCAGCAGGCTGAAAACCATGATTCGACGTCAATTACGTCAATGCTAAAACTCCGTCAGGCGCGGCTTGGCGTGAATGGCTGCATTATCTCAGGAATCGTTCTGGGATCTCTATA  >  NC_000913/946233‑946529

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