Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 11,156 G→A A67A (GCC→GCT
P109P (CCG→CCA
yaaW ←
mbiA →
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091311,1560GA86.1% 91.3 / 6.8 43A67A (GCC→GCT
P109P (CCG→CCA
yaaW
mbiA
putative enzyme‑specific chaperone YaaW
uncharacterized protein MbiA
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base A (14/23);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGT  >  NC_000913/10917‑11382
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
tCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTg                                                                                                                                                                                                                           >  1:163471/1‑251 (MQ=255)
 cATTCACCCTGGCATCGACCGCGGGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTg                                                                                                                                                                                                                           <  1:286077/250‑1 (MQ=255)
      aCCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGtt                                                                                                                                                                                                                     >  1:256519/1‑251 (MQ=255)
        ccTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTgc                                                                                                                                                                                                                   <  1:265133/251‑1 (MQ=255)
          tGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTgcgg                                                                                                                                                                                                                 >  1:291824/1‑251 (MQ=255)
           ggCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTgcggc                                                                                                                                                                                                                <  1:323050/251‑1 (MQ=255)
            gCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTgcggcg                                                                                                                                                                                                               <  1:321619/251‑1 (MQ=255)
                                   cgtcttctgattcTCGTCCATACTCTTCCATGTGTTTAGCATAAATGGTTCCAGTACCTGTTTATTAATCTCAAACTTAGACAGCTCTTTGTCGTCGTTTAGCGTCAATCGTTTAGAAACATCGAGCAAAATTGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGTCACCAAAATGCTGTAATGCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTGCCa                                                                                                                                                                                        <  2:202864/239‑1 (MQ=255)
                                     gCTTATGTTCCTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTTCGGCGATGTTGCTCGGGATGCCCTTCCAt                                                                                                                                                                                       <  1:122317/250‑1 (MQ=255)
                                     gCTTATGTTCCTCGTCCATATTCTTCCAGGGGTTACTCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCTATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCAt                                                                                                                                                                                       <  1:213282/250‑1 (MQ=255)
                                        tttGTTCCTCGTCCCTATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTTAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCATCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATTCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGa                                                                                                                                                                                    <  1:149992/248‑1 (MQ=255)
                                              ccTCGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAAc                                                                                                                                                                             <  1:378606/251‑1 (MQ=255)
                                                 cGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACa                                                                                                                                                                            <  1:337061/249‑1 (MQ=255)
                                                 cGTCCATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAAGTTACCGTGTTCGCGCAGTTTGTTGGCGATACTATCGCCTCCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAg                                                                                                                                                                           <  1:353339/250‑1 (MQ=255)
                                                    ccATATTCTTCCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTc                                                                                                                                                                        <  1:275591/250‑1 (MQ=255)
                                                     catattcttcCAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCACTCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGt                                                                                                                                                                      >  2:270683/5‑251 (MQ=255)
                                                              ccAGGTGTTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGGCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATc                                                                                                                                                              <  1:248831/250‑1 (MQ=255)
                                                                  gTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATAAGTAc                                                                                                                                                         >  2:75828/1‑251 (MQ=255)
                                                                         gCAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCTATACTATCGCCACCAAAATGCTGTAATTCTTCAGCAATCAGCTGCCAGTTGCGGCGATGGTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGgag                                                                                                                                                  <  1:163722/251‑1 (MQ=255)
                                                                          cAGAAATTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGGTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGATTGTTGGCGATACTCTCGCCACCAAAATGCTGTAATTCTCAAGCAATCAGATGACAGATGCGGCGAATCTGCACGGGATGCCCTACCATCGATTTAAACAGTTCGTTGCGTATCAGTAAGCTGgagg                                                                                                                                                 >  2:96702/1‑250 (MQ=255)
                                                                            gAAAGTGTTCCACTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAAGTTCCCGTGTCCGCGCAGTTTGTTGGCGATACTTTCGCCACCAAAATGCTGTAAGTCTCCAGCAATCATCTGCCAGTTGCGGCGTTGTTGCTCGGGATGCCCTTCCATCGATTTAAACATTTCGTTGCGTATCAGTTCGCTGGAGAgg                                                                                                                                              <  2:73749/251‑1 (MQ=255)
                                                                              aaTTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGGTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCg                                                                                                                                            <  1:359698/251‑1 (MQ=255)
                                                                               aTTGTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCg                                                                                                                                            >  1:26147/1‑250 (MQ=255)
                                                                                 tGTTCCAGTAACTGTTGCTCACTTTCAAAAGTAGACATCTCTTTGTCGGCTTTCAACTTCAATCGCTTTGAAACATCGAGCAACATGGCCCGATAAAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGACACCAAAATGCTGTAATTCTCCGGCAATCATCTGCCAGTTGCGGCGATGTTGCTAGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCCGTACGATGGAGCGGCGAgg                                                                                                                                         >  2:44055/1‑250 (MQ=255)
                                                                                  gTTCCAGTAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGGCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCTAGtt                                                                                                                                        <  1:254103/251‑1 (MQ=255)
                                                                                         tAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGCTACAATTTACCGTGTCCGCGCAGTTTGTTGGCGCTACTATCGCCACCAAAATGCTGTACTTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGAGGCCCTTCCATCGATTTAAACAGTTCGTTGTGCATCATTACGCTGGAGAGGCGAGTTTTACCtt                                                                                                                                 >  2:110411/1‑251 (MQ=255)
                                                                                         tAACTGTTGCTCAATCTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCt                                                                                                                                                      >  1:107447/1‑231 (MQ=255)
                                                                                                           aaaCGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAACCAGTTCGTTGCGTATCAGTACGCTGGGACGGCGAGTTTTGCCTTTTTCATTTTGGGGGAGc                                                                                                               >  2:270080/1‑250 (MQ=255)
                                                                                                             aCGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAAt                                                                                                            <  1:270080/251‑1 (MQ=255)
                                                                                                                       ctctTTGTCGGCTTTCAGCCTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACCGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTTAGCAATCggggggaa                                                                                                   >  2:181540/1‑245 (MQ=255)
                                                                                                                               cGGCTTTCAGCTTCAAGCGGTTTGAAACATCGAGCAGAATGGCCCGATACAAGTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCaa                                                                                           <  1:184675/250‑1 (MQ=255)
                                                                                                                                             aaTCGTTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACa                                                                            <  1:64789/251‑1 (MQ=255)
                                                                                                                                               tcgttttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGATTGTTGGCGATAATATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGAGAGCAATCGGGCGAAATTTGCCAACTGTGCCTCataacaa                                                                           >  2:99405/5‑244 (MQ=255)
                                                                                                                                               ccGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAAt                                                                          <  1:181540/250‑1 (MQ=255)
                                                                                                                                                cgttttGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATg                                                                         >  1:96716/4‑251 (MQ=255)
                                                                                                                                                           aTCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAAt                                                              >  1:63835/1‑251 (MQ=255)
                                                                                                                                                                    aaTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTg                                                     <  1:270683/251‑1 (MQ=255)
                                                                                                                                                                       ggCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGaa                                                   <  1:316429/250‑1 (MQ=255)
                                                                                                                                                                                     gaCCGTGTTCGCGCAGTATGTTGTCGATACTATCGCTACCAAAATGCTGTAATGCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGACTTATGCATTATGGGGGAGCAATCGGGCGAAATTTGACAACTGGTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCATGTAAtt                                     <  2:70130/249‑1 (MQ=255)
                                                                                                                                                                                      aCCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTa                                    <  1:44055/250‑1 (MQ=255)
                                                                                                                                                                                        cGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAAc                                  >  1:203656/1‑250 (MQ=255)
                                                                                                                                                                                                       tgttgGCGATACTATCGCCACGAACATGCTGTAATTCGCAAGCAATCAGCTGCCAGATGCGGCGATGATGCTCGGAATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCACTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGATGTGGCAATCGGGCGAAGATTGCAACCTGTTCCTCACTACAATGCTGAAGCAATGACAGCTCTGAATCTTTCAGGAAATTAACATTAATTTTTTTTggt                  >  2:64219/1‑250 (MQ=255)
                                                                                                                                                                                                                       gCCACCAAAATGCTGTAATTCTCCAGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGTATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCACTCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGACGAAAATCCAGATCTGAATCATTCAGGACATTAACATTACTTTTTTGGGGTTTCTATCTTCTGGCGt  >  1:236353/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                    
TCATTCACCCTGGCATCGACCGCGTGCAGAAACTCCTGCTTATGTTCCTCGTCCATTTTCTTCCAGGTATTACGCAGAAATTGTTCCAGTAACTGTTGCTCAATTTCAAACGTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCATCAGTACGCTGGAGAGGCGAGTTTTGCCTTTTTCATTATGGGTGAGCAATCGGGCGAAATTTGCCAACTGTTCCTCACTACAATGCTGAAGAAAATCCAGATCTGAATCATTCAGGTAATTAACATTCATTTTTTGTGGCTTCTATATTCTGGCGT  >  NC_000913/10917‑11382

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