Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,411,691 A→G E13E (GAA→GAG yhdJ → DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,411,6910AG100.0% 245.7 / NA 70E13E (GAA→GAGyhdJDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/37);  total (33/37)

GCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAG  >  NC_000913/3411557‑3411827
                                                                                                                                      |                                                                                                                                        
gCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATTGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTaa                                                                                                                                     <  1:859806/140‑1 (MQ=255)
gCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTaa                                                                                                                                     <  1:762233/140‑1 (MQ=255)
gCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTaa                                                                                                                                     <  2:191358/140‑1 (MQ=255)
    gAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTCTGGTAATGAGGCTAAGAcc                                                                                                                                 <  2:823953/140‑1 (MQ=255)
     aaCTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCa                                                                                                                                <  2:476944/140‑1 (MQ=255)
          aTTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCattatt                                                                                                                           <  1:871801/140‑1 (MQ=255)
               ggTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACgg                                                                                                                      >  2:162860/1‑140 (MQ=255)
               ggTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACgg                                                                                                                      >  1:733588/1‑140 (MQ=255)
               ggTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACgg                                                                                                                      >  2:360365/1‑140 (MQ=255)
                gTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGt                                                                                                                     >  2:659403/1‑140 (MQ=255)
                    gCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATg                                                                                                                 >  1:648154/1‑140 (MQ=255)
                          tGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTg                                                                                                           >  1:847706/1‑140 (MQ=255)
                             ttGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCg                                                                                                        <  2:638532/140‑1 (MQ=255)
                               gATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGaa                                                                                                      >  1:88934/1‑140 (MQ=255)
                               gATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGaa                                                                                                      >  1:319382/1‑140 (MQ=255)
                                          cgcTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGAt                                                                                           >  1:601049/1‑140 (MQ=255)
                                           gcTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATc                                                                                          >  1:130994/1‑140 (MQ=255)
                                              aCTCGGCATGGGGAAGCGCTTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCTGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATcccc                                                                                       <  2:672057/140‑1 (MQ=255)
                                               cTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGa                                                                                                       <  1:708773/123‑1 (MQ=255)
                                               cTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGa                                                                                                       >  2:708773/1‑123 (MQ=255)
                                               cTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCccg                                                                                      >  2:437527/1‑140 (MQ=255)
                                                  ggCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCccgccg                                                                                   <  1:162860/140‑1 (MQ=255)
                                                   gCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGa                                                                                  <  1:380570/140‑1 (MQ=255)
                                                   gCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGa                                                                                  <  1:901670/140‑1 (MQ=255)
                                                       ggggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAgt                                                                              >  2:426590/1‑140 (MQ=255)
                                                        gggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAgtg                                                                             <  1:290108/140‑1 (MQ=255)
                                                         ggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAgtgt                                                                            >  2:693698/1‑140 (MQ=255)
                                                         ggAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAgtgt                                                                            >  1:914878/1‑140 (MQ=255)
                                                          gAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTc                                                                           <  1:369568/140‑1 (MQ=255)
                                                                  ttttttATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCgatctgat                                                                   >  1:967681/1‑140 (MQ=255)
                                                                   tttttATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCgatctgatc                                                                  >  1:53259/1‑140 (MQ=255)
                                                                     tttATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCtt                                                                <  2:319382/140‑1 (MQ=255)
                                                                        ataGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGc                                                             <  2:771325/140‑1 (MQ=255)
                                                                           ggTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGa                                                          <  2:914878/140‑1 (MQ=255)
                                                                           ggTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGa                                                          <  2:827132/140‑1 (MQ=255)
                                                                             tgtCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                        <  1:437527/140‑1 (MQ=255)
                                                                             tgtCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                        <  2:130994/140‑1 (MQ=255)
                                                                             tgtCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                        <  1:197690/140‑1 (MQ=255)
                                                                             tgtCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                        <  2:648154/140‑1 (MQ=255)
                                                                               tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      >  2:455889/1‑140 (MQ=255)
                                                                               tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      <  2:53259/140‑1 (MQ=255)
                                                                               tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      >  2:440238/1‑140 (MQ=255)
                                                                               tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      <  2:381603/140‑1 (MQ=255)
                                                                               tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      >  2:661471/1‑140 (MQ=255)
                                                                               tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      <  2:847706/140‑1 (MQ=255)
                                                                               tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      >  2:707232/1‑140 (MQ=255)
                                                                               tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      <  1:360365/140‑1 (MQ=255)
                                                                               tCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      <  1:440238/140‑1 (MQ=255)
                                                                               ccacaAAGGGAGTGACCGTGAGGACAGGATGGGATCCAGCCCGGTTTGGTAATGGGGCTAAGACCATTATTCAGGGTGATGCGCTTGCGGACCTAAAAAAGCTCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACcca                                                      <  2:614172/139‑1 (MQ=255)
                                                                                       ggAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTTCCGACCCACCATATaa                                              >  2:315633/1‑140 (MQ=255)
                                                                                         aGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACa                                            <  1:426590/140‑1 (MQ=255)
                                                                                         aGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACa                                            <  1:470324/140‑1 (MQ=255)
                                                                                         aGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACa                                            >  2:814296/1‑140 (MQ=255)
                                                                                           tGACCTTGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATc                                          >  1:469820/1‑140 (MQ=255)
                                                                                             aCCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                        >  1:917484/1‑124 (MQ=255)
                                                                                             aCCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGAcc                                                        <  2:917484/124‑1 (MQ=255)
                                                                                             aCCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCgg                                        >  2:216444/1‑140 (MQ=255)
                                                                                                 tGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTaaa                                    >  1:973363/1‑140 (MQ=255)
                                                                                                  gagaACAGGATGTGAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTaaaa                                   >  1:359251/1‑140 (MQ=255)
                                                                                                            tgtgAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATgg                         >  2:236726/1‑140 (MQ=255)
                                                                                                              tgAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTc                       <  1:661471/140‑1 (MQ=255)
                                                                                                               gAACCGACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCt                      >  2:980421/1‑140 (MQ=255)
                                                                                                                    gACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCg                 >  2:198819/1‑140 (MQ=255)
                                                                                                                    gACCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCg                 >  1:803335/1‑140 (MQ=255)
                                                                                                                     aCCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGa                <  1:549263/140‑1 (MQ=255)
                                                                                                                     aCCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGa                <  1:455889/140‑1 (MQ=255)
                                                                                                                     aCCCGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGa                <  2:967681/140‑1 (MQ=255)
                                                                                                                      cccGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGaa               <  1:90199/140‑1 (MQ=255)
                                                                                                                        cGGTTTGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGc             <  1:236726/140‑1 (MQ=255)
                                                                                                                             tGGTAATGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGa        <  1:216444/140‑1 (MQ=255)
                                                                                                                                   tGAGGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAaagaag  >  2:480865/1‑140 (MQ=255)
                                                                                                                                      |                                                                                                                                        
GCATGAACTAATTCAGGTTAGCTAAATGCTTGATTAAAAAGGCGCTACTCGGCATGGGGAAGCGCCTTTTTTATAGGTGTCACAAAGGGAGTGACCATGAGAACAGGATGTGAACCGACCCGGTTTGGTAATGAAGCTAAGACCATTATTCACGGTGATGCGCTTGCCGAACTTAAAAAGATCCCCGCCGAAAGTGTCGATCTGATCTTTGCCGACCCACCATATAACATCGGTAAAAATTTTGATGGTCTGATCGAAGCCTGGAAAGAAG  >  NC_000913/3411557‑3411827

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