Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,095,328 C→A R37L (CGT→CTT)  insE4 ← IS3 element protein InsE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,095,3280CA98.0% 155.4 / NA 49R37L (CGT→CTT) insE4IS3 element protein InsE
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (29/19);  total (30/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCTGGCGTTTGAGACGTGCAATCTCGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGG  >  NC_000913/1095207‑1095622
                                                                                                                         |                                                                                                                                                                                                                                                                                                      
cAGCTTGCGTTTTAGACGTGCAATCTCGGTAGACATCTCCAGTGCACGTTCAGAAGACGTCTGCTGATTTGGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGCCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGc                                                                                                                      <  2:39637/300‑1 (MQ=9)
     ggCGTTTGAGACGTGCAATCTCGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGc                                                                                                                 >  2:150893/1‑300 (MQ=17)
     agcattttagaCTTGCAATCTCGGTAGACATCTCCAGTTCACGTTCAGAAGACGNCTGCTGATTTTGCTGTTTACTGCACCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTTATGTGAGTCACCTCTGACTNAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGc                                                                                                                 <  2:868355/292‑1 (MQ=255)
       cGTTTGAGACGTGCAATCTCGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGcac                                                                                                               >  1:496825/1‑300 (MQ=18)
       cGTTTGAGACGTGCAATCTCGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTGTGTCATGTGAGTCACCTATGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGcac                                                                                                               >  2:793147/1‑300 (MQ=18)
       cGTTTGAGACGTGCAATCTCGGTAGACATCTACAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGcac                                                                                                               >  1:714974/1‑300 (MQ=18)
          ttGAGACGTGCAATCTCGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACAcg                                                                                                            <  1:336927/300‑1 (MQ=21)
                 gTGCAATCTCGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGa                                                                                                     >  2:581504/1‑300 (MQ=37)
                      aTCTCGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGaa                                                                                                >  1:15043/1‑300 (MQ=39)
                       tctcGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAg                                                                                               >  1:713666/1‑300 (MQ=255)
                         tcGGGAGACAGCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCATTGGTAGAGTTTTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCCCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGAc                                                                                             <  2:584830/300‑1 (MQ=255)
                          cGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACg                                                                                            >  1:453072/1‑300 (MQ=255)
                                    ctcCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGc                                                                                  <  2:880098/300‑1 (MQ=255)
                                      ccAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCaa                                                                                >  2:319234/1‑300 (MQ=255)
                                      ccAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCATTCAGGCaa                                                                                >  2:318849/1‑300 (MQ=255)
                                                 tCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGAttct                                                                     >  2:85466/1‑298 (MQ=255)
                                                  cAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                                    <  1:354774/300‑1 (MQ=255)
                                                  cAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGACTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATc                                                                    >  1:411794/1‑300 (MQ=255)
                                                            ctgctgATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACt                                                          >  2:347640/1‑300 (MQ=255)
                                                                  aTTTTGCTGTCTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGatat                                                    <  2:74226/300‑1 (MQ=255)
                                                                     ttGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGc                                                 <  1:744068/300‑1 (MQ=255)
                                                                       gCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTNCGGGGTTTTTTACTGGTTGATACTGTNTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAATACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCtt                                               >  2:615881/1‑300 (MQ=255)
                                                                                      aGTTGTAGAGTTGTGATTCATACAGGTTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATGCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTTACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAgg                                <  2:370773/300‑1 (MQ=255)
                                                                                       gTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGt                               >  2:838058/1‑300 (MQ=255)
                                                                                       gTTGTAGAGTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGACAGAGCATCAGGt                               >  2:361723/1‑300 (MQ=255)
                                                                                       gTTGTAGAGTTGTGATTCATACAGGCTGAGTNCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGNTTGATACTGTTTTTGTCATGTGAGTCGCCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCGTCAGGt                               >  2:644110/1‑300 (MQ=255)
                                                                                       gTTGTAGAGTTGTACTTCATACATGCTGATTTTAAGTGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACNGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGt                               <  2:226171/300‑1 (MQ=255)
                                                                                            agagTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACNTCTGACTGAGAGTTTACTCACTTAGCNGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtga                          <  2:622562/300‑1 (MQ=255)
                                                                                              agTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatg                        <  1:439485/300‑1 (MQ=255)
                                                                                              agTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGAGGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAGTTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatg                        <  2:103143/300‑1 (MQ=255)
                                                                                               gTTGTGATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCtgatga                       >  1:333118/1‑300 (MQ=255)
                                                                                                 tgtgATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAAc                     <  1:108178/300‑1 (MQ=255)
                                                                                                    gATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGACGCGTTCAGCAAGCTTCAGGGCTTCACTGCGACATTCAGGCGAATGCTGTTTATGGGGTCTTTTACTGGCCGATACCGTTTTTGTCATGTGAGTCAACTCTGACTGCGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGCTCAGTGCCTGGCACACGCGTTTGATGGAAGAAGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCAAGAGCCTCAGGTTCTGACGACCGTa                  >  2:297213/1‑300 (MQ=18)
                                                                                                    gATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTa                  >  1:428285/1‑300 (MQ=255)
                                                                                                    gATTCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTa                  >  2:124842/1‑300 (MQ=255)
                                                                                                      ttCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAAc                <  1:854199/300‑1 (MQ=255)
                                                                                                      ttCATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTGACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAAc                <  1:861637/300‑1 (MQ=255)
                                                                                                        cATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCTATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCGTGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAg              <  1:845470/300‑1 (MQ=255)
                                                                                                        cATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAg              <  1:364391/300‑1 (MQ=255)
                                                                                                        cATACAGGCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAg              <  2:829795/300‑1 (MQ=255)
                                                                                                              ggCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCaaaaa        >  2:724811/1‑300 (MQ=255)
                                                                                                              ggCTGAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTAGCCAGAGCATCAGGTTCTGATGAACGTAACAGCaaaaa        >  2:53220/1‑300 (MQ=255)
                                                                                                                  gAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATTAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGa    >  2:6927/1‑300 (MQ=255)
                                                                                                                  gAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATTGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGa    >  2:263057/1‑300 (MQ=255)
                                                                                                                  gAGTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGa    >  1:424533/1‑300 (MQ=255)
                                                                                                                    gTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACNTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTNGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAgg  >  2:745416/1‑300 (MQ=255)
                                                                                                                    gTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAgg  >  1:352668/1‑300 (MQ=255)
                                                                                                                    gTTCAAGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAgg  >  2:576569/1‑300 (MQ=255)
                                                                                                                    gTTCAAGGGCTGCGGCAGTAACACCGACGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAgg  >  1:502212/1‑300 (MQ=255)
                                                                                                                         |                                                                                                                                                                                                                                                                                                      
CAGCTGGCGTTTGAGACGTGCAATCTCGGTAGACATCTCCAGTTCACGTTCAGAAGACGTCTGCTGATTTTGCTGTTTACTGCGCCAGTTGTAGAGTTGTGATTCATACAGGCTGAGTTCACGGGCTGCGGCAGTAACACCGATGCGTTCAGCAAGCTTCAGGGCTTCACTGCGAAATTCAGGCGAATGCTGTTTACGGGGTTTTTTACTGGTTGATACTGTTTTTGTCATGTGAGTCACCTCTGACTGAGAGTTTACTCACTTAGCCGCGTGTCCACTATTGCTGGGTAAGATCAGTGCCTGGCACACGCGTTTGATGGAAGACGGCGTTCAGGCAATGCGCGATTATCTTGCCGGACTGGATATCGCTTCGCCAGAGCATCAGGTTCTGATGAACGTAACAGCAAAAAGCGAGG  >  NC_000913/1095207‑1095622

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