Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,597,272 T→A 5.7% pseudogene (716/3861 nt) yneO ← pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,597,2720TA 5.7% 430.8 / 12.7 123pseudogene (716/3861 nt)yneOpseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base T (59/57);  new base A (4/3);  total (63/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

Only 100 of 127 total aligned reads displayed.

AGACATTAAGTTCACCGTCGCCAGCATCACCGATATAGATATAGCGGAAAGCTTCGCCCGTTCCGAGGAAATTCCAATGCCCTTTATCGGTCACGTTGACGACACCGTGGGAACCGTCCTGAACGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCATAGCTACCTATTTCGAATAATTCGGTCGTCAGAAC  >  NC_000913/1597003‑1597540
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                            
agaCATTAAGTTCACCGTCGCCAGCATCACCGATATAGATATAGCGGAAAGCTTCGCCCGTTCCGAGGAAATTCCAATGCCCTTTATCGGTCACGTTGACGACACCGTGGGAACCGTCCTGAACGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGCTGGCATTGCCa                                                                                                                                                                                                                                                                         >  2:242148/1‑275 (MQ=255)
                                 tatagatatagCGGAAAGCTTCGCCCGTTCCGAGGAAATTCCAATGCCCTTTATCGGTCACGTTGACGACACCGTGGGAACCGTCCTGAACGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAcccc                                                                                                                                                                                                                                        <  2:292546/275‑1 (MQ=255)
                                                                                           aaCGTTGACCAAACCCTGGGAACCGTCCTGAACGCCCACCAAGGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAgg                                                                                                                                                                              <  4:189240/274‑1 (MQ=255)
                                                                                                     acacCGTGGGAACCGTCCTGAACGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattatt                                                                                                                                                                    <  2:386014/275‑1 (MQ=255)
                                                                                                                         aaCGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCAt                                                                                                                                                                                                                                                                              <  6:676631/149‑1 (MQ=255)
                                                                                                                            gCCCACAAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTg                                                                                                                                                                                                                                                                            <  7:306526/148‑1 (MQ=255)
                                                                                                                              ccaccaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCa                                                                                                                                                                                                                                                                         >  5:539800/1‑149 (MQ=255)
                                                                                                                              ccaccaATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGcc                                                                                                                                                                                                                      >  4:315215/1‑200 (MQ=255)
                                                                                                                                  caATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGa                                                                                                                                                                                                                                                                     >  6:455319/1‑149 (MQ=255)
                                                                                                                                  caATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGa                                                                                                                                                                                                                                                                     >  5:667687/1‑149 (MQ=255)
                                                                                                                                  caATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTccc                                                                                                                                                                                                                  <  2:315215/200‑1 (MQ=255)
                                                                                                                                   aaTCAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGAt                                                                                                                                                                                                                                                                    >  5:36604/1‑149 (MQ=255)
                                                                                                                                       aaTATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAt                                                                                                                                                                                                                                                                 <  8:184191/148‑1 (MQ=255)
                                                                                                                                       aaTATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATc                                                                                                                                                                                                                                                                >  5:536452/1‑149 (MQ=255)
                                                                                                                                            tCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAtt                                                                                                                                                                                                                                                           <  6:666351/149‑1 (MQ=255)
                                                                                                                                            tCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAt                                                                                                                                                                                                                                                            <  7:301099/148‑1 (MQ=255)
                                                                                                                                             cTTTGGTATTAAGCAGTCCATTTTTGGAAATATTGACTGTGCCATTACCGAAATATCGATTATGAGGTCGGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCAt                                                                                                                        <  1:14142/279‑1 (MQ=255)
                                                                                                                                              tttGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTAttt                                                                                                                                                                                                                                                          <  7:493692/148‑1 (MQ=255)
                                                                                                                                               ttGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTg                                                                                                                             >  4:265938/1‑272 (MQ=255)
                                                                                                                                                tGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATTCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCa                                                                                                                                                                                                                                                       <  6:923928/149‑1 (MQ=255)
                                                                                                                                                tGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCa                                                                                                                                                                                                                                                       <  6:536452/149‑1 (MQ=255)
                                                                                                                                                 gTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCa                                                                                                                                                                                                                                                       <  7:966124/148‑1 (MQ=255)
                                                                                                                                                        aTCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTaac                                                                                                                                                                                                                                               >  5:820881/1‑149 (MQ=255)
                                                                                                                                                                 ttattGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAAccccc                                                                                                                                                                                                                                       >  8:344005/1‑148 (MQ=255)
                                                                                                                                                                   attGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGc                                                                                                      >  4:291990/1‑275 (MQ=255)
                                                                                                                                                                       gAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACCACAACCTGGCCATTACTGCCCGCCTGATCGCCTAAAAtt                                                                                              >  1:221563/1‑278 (MQ=255)
                                                                                                                                                                            atTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAg                                                                                                                                                                                                                           >  6:239404/1‑149 (MQ=255)
                                                                                                                                                                             ttGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                                                                                                                                                                                          <  5:257237/149‑1 (MQ=255)
                                                                                                                                                                              tGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGt                                                                                                                                                                                                                          <  8:9262/148‑1 (MQ=255)
                                                                                                                                                                              tGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTc                                                                                                                                                                                                                         <  5:239404/149‑1 (MQ=255)
                                                                                                                                                                              tGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGCCAACCTGCCCATTACTGCCCGCCTGATAGCCTACAATGGcgacga                                                                                       >  1:177920/1‑279 (MQ=255)
                                                                                                                                                                                tctttgCCAGTACTGATACAGCCATTATATACTCGCCGTACCGAGATAGCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCg                                                                                                                                                                                                                        <  7:783553/144‑1 (MQ=255)
                                                                                                                                                                                  tgtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCt                                                                                                                                                                                                                     >  6:344510/1‑149 (MQ=255)
                                                                                                                                                                                  tgtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGcgacga                                                                                       >  2:554739/1‑275 (MQ=255)
                                                                                                                                                                                  tgtgCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCACCCCAGTTCCTTGATTACCAAATTGAATTTCACTTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCCCATGCTAGCCTACAATGCccacga                                                                                       >  4:45491/1‑275 (MQ=255)
                                                                                                                                                                                      ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTcccc                                                                                                                                                                                                                 <  6:1031326/149‑1 (MQ=255)
                                                                                                                                                                                      ccATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCCCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGAAAGCCTAAAAGGGCGCAGAttt                                                                                    >  4:581222/1‑273 (MQ=255)
                                                                                                                                                                                        aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATAACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAg                                                                                                                                                                                                                                                      >  6:640491/1‑110 (MQ=255)
                                                                                                                                                                                        aTTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATAACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAg                                                                                                                                                                                                                                                      <  5:640491/110‑1 (MQ=255)
                                                                                                                                                                                         ttACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAg                                                                                                                                                                                                               <  7:885316/148‑1 (MQ=255)
                                                                                                                                                                                            ccGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                                                                                                                                                                           <  6:820881/149‑1 (MQ=255)
                                                                                                                                                                                             cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                                                                                                                                                                                                     >  5:866423/1‑106 (MQ=255)
                                                                                                                                                                                             cGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGc                                                                                                                                                                                                                                                     <  6:866423/106‑1 (MQ=255)
                                                                                                                                                                                              gAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCtt                                                                                                                                                                                                         <  5:932173/149‑1 (MQ=255)
                                                                                                                                                                                               aaaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                                                                                                                                                                        <  5:43500/149‑1 (MQ=255)
                                                                                                                                                                                                aaTATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTg                                                                                                                                                                                                        <  8:590466/148‑1 (MQ=255)
                                                                                                                                                                                                 atagCCATTATAGAGTCTGCGGACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGt                                                                        <  2:581222/275‑1 (MQ=255)
                                                                                                                                                                                                        ttATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                                                                                                                                                                                               >  6:355834/1‑149 (MQ=255)
                                                                                                                                                                                                        ttATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAAt                                                                                                                                                                                               >  6:207940/1‑149 (MQ=255)
                                                                                                                                                                                                            agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAAt                                                                                                                                                                                                                                                                 >  7:851787/1‑79 (MQ=255)
                                                                                                                                                                                                            agagTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGa                                                                                                                                                                                           >  5:276398/1‑149 (MQ=255)
                                                                                                                                                                                                                tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATATCCTAAAATGGCGACGATTTAACTCGTGACTTAACCCTAATCCGTaa                                                         >  4:112309/1‑275 (MQ=255)
                                                                                                                                                                                                                tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTaa                                                         >  2:296690/1‑275 (MQ=255)
                                                                                                                                                                                                                tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTACCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTaa                                                         >  2:408808/1‑275 (MQ=255)
                                                                                                                                                                                                                tCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACAAACCCCCCCTCGCGAATAGTCGCCTACCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGATTCATTATTTTTTTTTAGCCACTCTCCACCCTTTTCAACGACACCCTGGCCTTTACTCCCCGCCTGTTAGCCTAATATGCCCACTATTGACCCCGTACCGTAACCCTTTTCCGTaa                                                         >  4:626889/1‑275 (MQ=255)
                                                                                                                                                                                                                  tacGTACCGTGATGGCCGGGCTGTTATCCTGGACATTCAGAGTTCCCATACCGGTAGCCTTGCCACCGATAATCGTAATTTCAGCCATAACCAAACCCCCCTCGCGAATAGTCGCCTCCCCAATTCCTTAATTACCAATGTGAAATTCAATTGAGGAATAATTATTTTTTATTAGCCACTCACCACCCTTTTCAGCGACAGCCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGGCGTAACCCTTATCCGTAata                                                       <  4:527407/273‑1 (MQ=255)
                                                                                                                                                                                                                  tGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGa                                                                                                                                    >  2:676388/1‑198 (MQ=255)
                                                                                                                                                                                                                    cGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAt                                                                                                                                                                                   <  5:125657/149‑1 (MQ=255)
                                                                                                                                                                                                                      tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAAtt                                                                                                                                                                                  <  8:308827/148‑1 (MQ=255)
                                                                                                                                                                                                                      tACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc                                                                                                                                <  4:676388/198‑1 (MQ=255)
                                                                                                                                                                                                                        ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                                                                                                                                                               <  5:436686/149‑1 (MQ=255)
                                                                                                                                                                                                                        ccGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAg                                                                                                                                                                               <  6:541946/149‑1 (MQ=255)
                                                                                                                                                                                                                         cgtgtgtcagtggTCTTCATGCTGCACATTCAGNNNCCCGCTACCGGTGGCATTGCCACCGCTCATCGTTTTTTACGCCGTAAACAACCCCCCCCCGCGAATAGTCGCCTCCCCAGTTACTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTCTATTAGCCACTCGCCACCCTTTTCAACGACAACCTTGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAAcc                                            <  3:221563/267‑1 (MQ=255)
                                                                                                                                                                                                                           tgatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATAGAGGaa                                                                                                                                                                            >  6:785042/1‑149 (MQ=255)
                                                                                                                                                                                                                            gatgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAggga                                                                                                                                                                            >  8:875767/1‑146 (MQ=255)
                                                                                                                                                                                                                             atgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATc                                                                                                                                                                          >  5:271541/1‑149 (MQ=255)
                                                                                                                                                                                                                              tgaCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGACATTCAATTGAGGAATCCTTATTTTTTATTAGCCACTCGCTACCCTTTTCAACGACAACCTGGCCATTACTTCCCGCCTGATAGCCTAAAATGGCTACGATTGCACTCGTACCGCAACCCTTAACCGTAATATTTAATCAAcc                                            >  4:221492/1‑274 (MQ=255)
                                                                                                                                                                                                                               gatagagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                                                                                                                                                                        >  6:424182/4‑149 (MQ=255)
                                                                                                                                                                                                                                aCAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCat                                                                                                                                                                        <  8:722020/148‑1 (MQ=255)
                                                                                                                                                                                                                                 cAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatt                                                                                                                                                                       >  8:940598/1‑148 (MQ=255)
                                                                                                                                                                                                                                  agagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCatta                                                                                                                                                                      <  7:344005/148‑1 (MQ=255)
                                                                                                                                                                                                                                   gagTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCattat                                                                                                                                                                     >  7:151655/1‑148 (MQ=255)
                                                                                                                                                                                                                                     gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGCATCATTAtttt                                                                                                                                                                  <  6:271541/149‑1 (MQ=255)
                                                                                                                                                                                                                                     gTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTAtttt                                                                                                                                                                  >  5:626911/1‑149 (MQ=255)
                                                                                                                                                                                                                                        ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                                                                                                                                                               <  6:556981/149‑1 (MQ=255)
                                                                                                                                                                                                                                        ttGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTtta                                                                                                                                                               >  5:917017/1‑149 (MQ=255)
                                                                                                                                                                                                                                          gATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttat                                                                                                                                                              <  8:68476/148‑1 (MQ=255)
                                                                                                                                                                                                                                           aTCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTttatta                                                                                                                                                            <  5:993880/149‑1 (MQ=255)
                                                                                                                                                                                                                                                gCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCAc                                                                                                                                                       >  6:1045897/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      tCAGGGTTTCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTTCCCAGTTCATTGATTACCAATTTGAAATCCAATTGAGGAATCATTATTTTTTATTAGCCACTCGcc                                                                                                                                                  >  8:561639/1‑148 (MQ=255)
                                                                                                                                                                                                                                                          ggTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCt                                                                                                                                             >  5:484542/1‑149 (MQ=255)
                                                                                                                                                                                                                                                            ttCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCtt                                                                                                                                            <  8:1022817/148‑1 (MQ=255)
                                                                                                                                                                                                                                                              ccGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACCCGTGACGTTACCCTTCATCGTAATCTTTAATTAACCTTTGCACTATCTACCTATTTCAAATAATTCGgt           >  4:644127/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                gATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCaa                                                                                                                                       >  5:190242/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 ttatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACt                                                                                                                                                      <  7:181606/129‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ttatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACt                                                                                                                                                      >  8:181606/5‑133 (MQ=255)
                                                                                                                                                                                                                                                                 ttatcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAAc                                                                                                                                      >  5:5026/5‑149 (MQ=255)
                                                                                                                                                                                                                                                                 aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCaa                                                                                                                                       <  7:309137/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aTACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAAc                                                                                                                                      <  5:207940/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tatcGGTGGCAATGCCACCGGTAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGCTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCTACTCCCCACCCTTTTCAAc                                                                                                                                      >  7:830006/4‑148 (MQ=255)
                                                                                                                                                                                                                                                                    tcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                                                                                                                                                                      >  5:201149/2‑98 (MQ=255)
                                                                                                                                                                                                                                                                    tcGGTGGCAATGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTc                                                                                                                                                                                      <  6:201149/97‑1 (MQ=255)
                                                                                                                                                                                                                                                                        tGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAAc                                                                                                                                >  7:129252/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                        tGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTTCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCATAGCTACCTATTTCGAATAATTCGGTCGTCAGAAc  >  2:3757/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                         ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTcc                                                                                                                                                                                                           >  4:303072/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                         ggCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCaa                                                                                                                                                                                    >  2:495787/1‑95 (MQ=255)
                                                                                                                                                                                                                                                                           cATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTgg                                                                                                                            <  6:667687/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aTTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGc                                                                                                                           >  6:1022267/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                             ttGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGa                                                                                                                                                                                                       <  2:303072/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                             ttGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGa                                                                                                                                                                                <  4:495787/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                            
AGACATTAAGTTCACCGTCGCCAGCATCACCGATATAGATATAGCGGAAAGCTTCGCCCGTTCCGAGGAAATTCCAATGCCCTTTATCGGTCACGTTGACGACACCGTGGGAACCGTCCTGAACGCCCACCAATGAATATTCTTTGTTATTAATCAGTCCATTATTGGAAATATTGACTGTGCCATTACCGAAATATCCATTATAGAGTCTGCGTACCGTGATGACAGAGTCTTGATCCTGCACATTCAGGGTTCCGATACCGGTGGCATTGCCACCGATAATCGTATTTTCAGCCGTAACCAACCCCCCCTCGCGAATAGTCGCCTCCCCAGTTCCTTGATTACCAATTTGAAATTCAATTGAGGAATCATTATTTTTTATTAGCCACTCGCCACCCTTTTCAACGACAACCTGGCCATTACTGCCCGCCTGATAGCCTAAAATGGCGACGATTGAACTCGTGACGTAACCCTTATCCGTAATATTTAATGAACCTGTGCCATAGCTACCTATTTCGAATAATTCGGTCGTCAGAAC  >  NC_000913/1597003‑1597540

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