Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,635,658 T→C D46D (GAT→GAC creA → PF05981 family protein CreA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,635,6580TC94.2% 153.0 / 0.3 52D46D (GAT→GACcreAPF05981 family protein CreA
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (32/17);  total (34/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.66e-01

TGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGC  >  NC_000913/4635425‑4635903
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
tGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGGTGTGGAAGCCTTTGACGATCCCGATGTGaaa                                                                                                                                                                                                                                         <  1:281082/249‑1 (MQ=255)
                      aTTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCgg                                                                                                                                                                                                                 >  1:322436/1‑251 (MQ=255)
                              cTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTAGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAAc                                                                                                                                                                                                         >  1:533866/1‑251 (MQ=255)
                                 aaTCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTGATAGATAAAAATGGAAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGACTTTGACGATCCCGCTGTGAAGAATGTAACCTGTTATGTGGGCCGGGCGAAAAcagg                                                                                                                                                                                                      >  1:245547/1‑248 (MQ=255)
                                              tttACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAggg                                                                                                                                                                                         >  1:55306/1‑251 (MQ=255)
                                               ttACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAgggg                                                                                                                                                                                        >  1:91211/1‑251 (MQ=255)
                                                          aaCTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTg                                                                                                                                                                             >  1:365037/1‑251 (MQ=255)
                                                           actcaTTGCTCGTGTAATAGATAAAAATGGAAACACCATGAAATACACGCTTTTGCTCCTGTCTTTAAGCCTGATAATGCTGGGGCCTATGGCTCATACAGAAGCGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCAGGATCACAACATTGTTGTGGAAGCCTTTGACGCTCACGATGGGAAAAATGTCACCTGATATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTgg                                                                                                                                                                            >  2:65349/5‑251 (MQ=255)
                                                            cTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGAtg                                                                                                                                                                                                                                              <  2:315162/184‑1 (MQ=255)
                                                            cTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGAtg                                                                                                                                                                                                                                              >  1:315162/1‑184 (MQ=255)
                                                                 tGCTCGAGTAATCGACAAAAATGATAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAACCTTTGCCGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGGAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGacg                                                                                                                                                                      >  2:280008/1‑249 (MQ=255)
                                                                   cTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGGGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGAGGAAggt                                                                                                                                                                    >  2:314933/1‑249 (MQ=255)
                                                                               aTAAAAATGGTAACAATTAGCAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGTGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCTACACCGTAGTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATACCGATGTGAAAAATGTCACCTGTGATGTGAGTCGTGCGAAAACCGGTGGTATTAAAGGAGGATTGGGTCTGGCGGAAGATTCCTCCGATGCg                                                                                                                                                        <  2:32488/251‑1 (MQ=255)
                                                                                  aaaaTGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTATGGCGGAAGATACCTCCGATGCGGct                                                                                                                                                     >  1:107079/1‑250 (MQ=255)
                                                                                        gTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGCTAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAACTGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGGACCCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCGGGCGGAAGATACCTCCGGTGCGGCTATTTCt                                                                                                                                               >  2:338298/1‑251 (MQ=255)
                                                                                            cAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGCAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTc                                                                                                                                           >  2:183103/1‑251 (MQ=255)
                                                                                             aaTATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTAACCTGTTATGTGAGCAGGGCGACAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTc                                                                                                                                           >  2:139473/1‑250 (MQ=255)
                                                                                                        cAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTTTTTCTTGTCAGCAAGTCGGGc                                                                                                                               >  2:267117/1‑251 (MQ=255)
                                                                                                          aGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAACATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCCAGTCGGGCCg                                                                                                                             >  1:137845/1‑251 (MQ=255)
                                                                                                             aTTTGATCCTGGCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGAACCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGAtt                                                                                                                          <  2:365037/251‑1 (MQ=255)
                                                                                                             aTTTGATCCTGGCTTTAAGCCTGATAAAGCTGGGGCCATTTGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGAtt                                                                                                                          <  2:137845/251‑1 (MQ=255)
                                                                                                               ttGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTg                                                                                                                         <  1:271767/250‑1 (MQ=255)
                                                                                                                       gTCTTTAAGCCTGCTAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCAGATTGAACTGTCgg                                                                                                                >  2:238911/1‑251 (MQ=255)
                                                                                                                       gTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGACGATTGAACTGTCgg                                                                                                                >  1:428365/1‑251 (MQ=255)
                                                                                                                        tCTTAAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAACTGTCCCCTGTAATGTGAGCCGGGCGAAACCAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGCAAACTCCGATGCGGCTATTTATTGGCAGCAAGTAGGGACGAGTGAACTATCgg                                                                                                                >  2:542619/1‑250 (MQ=255)
                                                                                                                           ttAAGCCTGATAAGGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTGCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGGGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATc                                                                                                             <  2:55306/250‑1 (MQ=255)
                                                                                                                              aGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTAt                                                                                                         <  1:197910/251‑1 (MQ=255)
                                                                                                                                ccTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGctat                                                                                                                                                   <  1:114787/207‑4 (MQ=255)
                                                                                                                                     tAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGGATTAAAAAc                                                                                                  >  2:313885/1‑251 (MQ=255)
                                                                                                                                         gCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGc                                                                                               <  1:478778/250‑1 (MQ=255)
                                                                                                                                          cTGGGGCCATTGGCTCATGCAGAAGAGTTTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCt                                                                                                                                               >  2:484767/1‑201 (MQ=255)
                                                                                                                                          cTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCt                                                                                                                                               <  1:484767/201‑1 (MQ=255)
                                                                                                                                              ggCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCt                                                                                         >  1:5528/1‑251 (MQ=255)
                                                                                                                                                 cATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCAGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGAACCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTc                                                                                                                  <  1:303966/223‑1 (MQ=255)
                                                                                                                                                 cATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCAGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGAACCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTc                                                                                                                  >  2:303966/1‑223 (MQ=255)
                                                                                                                                                     ggCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGAGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCg                                                                                  >  1:499496/1‑251 (MQ=255)
                                                                                                                                                              agaagaGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTAt                                                                         <  2:105634/251‑1 (MQ=255)
                                                                                                                                                               gaagaGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTAtt                                                                        <  2:5528/251‑1 (MQ=255)
                                                                                                                                                                    gaTTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCt                                                                                         >  2:304862/1‑229 (MQ=255)
                                                                                                                                                                    gaTTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCt                                                                                         <  1:304862/229‑1 (MQ=255)
                                                                                                                                                                        ggTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAACTTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGGAGTATTCAAAAAACg                                                              >  1:574908/1‑251 (MQ=255)
                                                                                                                                                                             ggTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCAACTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGt                                                         >  1:443697/1‑251 (MQ=255)
                                                                                                                                                                                 gACACCGAATTTCAAATGATCGGCCCGAATCACCAAATTGTTGTGGAAGCCTTTGATGATACCGATGTGAACAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGCTGCGGCCCTTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGAGCGTATTAAAAACGGCAAAGCGCAGGGCGAGGTAGTCTTCCAAAAAACGCACTTccc                                                      >  2:338737/1‑251 (MQ=255)
                                                                                                                                                                                           ttAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTa                                             >  1:344595/1‑249 (MQ=255)
                                                                                                                                                                                             aaaaTGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCCAAGCTCAGGGCGAGGGAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTc                                         >  1:636349/1‑251 (MQ=255)
                                                                                                                                                                                             aaaaTGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTc                                         >  1:619387/1‑251 (MQ=255)
                                                                                                                                                                                               aaggATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTc                                                        <  2:235637/231‑1 (MQ=255)
                                                                                                                                                                                               aaTGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTc                                                        >  1:235637/1‑234 (MQ=255)
                                                                                                                                                                                                   aTCGGCCCGGATCACACAATTGTTGTGGAAGTCTTTGACGATCCCGATGTGAAAAAGGTCACCGGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACa                                   <  2:2438/251‑1 (MQ=255)
                                                                                                                                                                                                          cGGATCACAAAATTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTg                            <  1:141845/251‑1 (MQ=255)
                                                                                                                                                                                                                     aTTGTTGTGGAAGCCTTTGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTatg                  >  1:326840/1‑250 (MQ=255)
                                                                                                                                                                                                                                    tttGACGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACAGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCTATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGACGCCAaacgcaacgc  >  1:432344/1‑251 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
TGACGGAAATAAAAGTATTGAGATTTTGTTCTTAATCAATATGTTATTTACCGTGACGAACTAATTGCTCGTGTAATAGATAAAAATGGTAACAATATGAAATACAAGCATTTGATCCTGTCTTTAAGCCTGATAATGCTGGGGCCATTGGCTCATGCAGAAGAGATTGGTTCGGTCGACACCGTATTTAAAATGATCGGCCCGGATCACAAAATTGTTGTGGAAGCCTTTGATGATCCCGATGTGAAAAATGTCACCTGTTATGTGAGCCGGGCGAAAACCGGTGGTATTAAAGGGGGATTGGGTCTGGCGGAAGATACCTCCGATGCGGCCATTTCTTGTCAGCAAGTCGGGCCGATTGAACTGTCGGATCGTATTAAAAACGGCAAAGCTCAGGGCGAGGTAGTATTCAAAAAACGCACGTCCCTGGTCTTTAAGTCGTTACAGGTCGTGCGCTTTTATGATGCCAAACGCAACGC  >  NC_000913/4635425‑4635903

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