Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,555,047 2 bp→TG coding (1963‑1964/2706 nt) malT → DNA‑binding transcriptional activator MalT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,555,0470GT94.9% 86.0 / ‑2.9 40V655C (GTC→TGC) malTDNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (17/20);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01
*NC_0009133,555,0480TG94.9% 105.9 / ‑3.5 40V655C (GTC→TGC) malTDNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (17/20);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CATTCGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAAGT‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGC  >  NC_000913/3554806‑3555284
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                             
cATTCGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTgg                                                                                                                                                                                                                                      >  2:111104/1‑250 (MQ=255)
   tCGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCTGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTTGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGAtt                                                                                                                                                                                                                                    <  1:194374/250‑1 (MQ=255)
     gttcgcAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTAc                                                                                                                                                                                                                                 <  1:398509/248‑1 (MQ=255)
     gTGCGCAGCTGTTATGGGCCTGGGCGCAGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGGAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCAGGCAATGATGATTCAATGCTCGCAGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACAGTCTGGAAAACCTGCTGGGGAATGGCAAATCTCAAAGCGACTGGCTCTCGAACGCCACACAATG‑‑CCGGGTGGTTTAc                                                                                                                                                                                                                                 >  2:325248/1‑251 (MQ=255)
        cgcAGCTGCTATGGGCCTAGGCGCGGCTGGATGAAGCCGAAGCGCCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCCAAGCAACAGCTTCAGTGCCTGGCCATGTTGTTTCAAAGCTCGCTGGCCCGTGGTGATTTAGGTAACGCCCGTAGCCAGCTGAACCGTCTGGAAACCCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTATAACCCTA‑ACAAATG‑CCGGGGGATTTACTgg                                                                                                                                                                                                                              >  2:47257/1‑251 (MQ=255)
         acAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCCACAGCTTCAGTGCCTGGCAATGTTGATTAAATGCTCGCTGGACCGTGGAGATTTAGATAACGCCCGTAGCCAGCTGAAACGACTGGAAAACCTGCTGAGGAATGGCAAATAAAACAGCAACTGGATACATAACGTCA‑ACAAATG‑CCGGGTGGTTTACTGGc                                                                                                                                                                                                                             >  2:391255/2‑251 (MQ=255)
                 tATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGAcc                                                                                                                                                                                                                     >  2:394846/1‑251 (MQ=255)
                      gCCTGGGCTCGGCTGGATGAACCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGAGATTTAGGTAACGCCAGTAGCCAGCTGAACCGTCTGGCAAACCTGCTGGGGAATGGCAAATATCACAGCGCCTGGTTCTCTAACGCCC‑A‑AAACTGCCGGGTGCTTTCCTGGTAACTGACCGGCg                                                                                                                                                                                                                 >  2:231521/1‑250 (MQ=255)
                       ccTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGaagatcggg                                                                                                                                                                                                               >  2:163297/1‑243 (MQ=255)
                       ccTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGa                                                                                                                                                                                                                       <  1:163297/243‑1 (MQ=255)
                          gggCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGAGCTCTAACGCCA‑ACACATG‑CCGGGGGATTTACTGGCAAATGACCAGCGATaa                                                                                                                                                                                                             >  2:474780/1‑250 (MQ=255)
                                           gCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTg                                                                                                                                                                                           <  1:248169/251‑1 (MQ=255)
                                             cGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGt                                                                                                                                                                                         <  1:437728/251‑1 (MQ=255)
                                              gAAGCGTCGGCGCGTAGCGGGATTGAAGCATTGTCGTCTTACCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGtt                                                                                                                                                                                        <  1:179557/251‑1 (MQ=255)
                                                aGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGc                                                                                                                                                                                      >  2:147988/1‑251 (MQ=255)
                                                             aGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATAc                                                                                                                                                                               >  1:354283/1‑245 (MQ=255)
                                                             aGCGGGATTGAAGTCTTGCCTTCTTACCAGCCCCAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATAc                                                                                                                                                                               <  2:354283/245‑1 (MQ=255)
                                                                   aTTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAAcc                                                                                                                                                                      >  1:358875/1‑248 (MQ=255)
                                                                       aaGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGttt                                                                                                                                                               <  1:991/251‑1 (MQ=255)
                                                                                 gtcTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAAGT‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATCCGGCTAAACCAGATGTTGCGAACAAcc                                                                                                                                                     >  2:157626/1‑251 (MQ=255)
                                                                                     taccAGCCACAGCAACAGCTTCAGGGCCTGGAAATGTTGATTCAATGCGCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCAA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCa                                                                                                                                                    <  1:134854/245‑1 (MQ=255)
                                                                                      aTCAGCCACAGCAACAGCTTAAGTGCCTGGCAATGTTGATTCAAGGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAAGGGCAAATATCACAGCGAATGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTc                                                                                                                                                <  2:425081/251‑1 (MQ=255)
                                                                                                     aGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATgg                                                                                                                                 <  2:100831/251‑1 (MQ=255)
                                                                                                     aGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATgg                                                                                                                                 >  1:100831/1‑251 (MQ=255)
                                                                                                           aGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACt                                                                                                                                                                                            <  2:151205/186‑1 (MQ=255)
                                                                                                           aGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACt                                                                                                                                                                                            >  1:151205/1‑186 (MQ=255)
                                                                                                           aGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAAc                                                                                                                           <  1:111104/251‑1 (MQ=255)
                                                                                                            gTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCTTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACa                                                                                                                          <  1:474780/251‑1 (MQ=255)
                                                                                                                                     gcGCTGGCCCTGGGTAATTTAGAAAACGCCCGTAGCCAGCTGAACCTTCGGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTggg                                                                                                 <  2:120549/250‑1 (MQ=255)
                                                                                                                                                           aTAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACAGGCAGAAAtt                                                                           >  1:76164/1‑251 (MQ=255)
                                                                                                                                                                           gCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCATAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGTCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCa                                                           <  1:61083/251‑1 (MQ=255)
                                                                                                                                                                                        gAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCg                                              >  1:74007/1‑251 (MQ=255)
                                                                                                                                                                                        gAAAACCTGCGGGGGAATGGCAAATAACACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGGGATTTACGGGCAAAGGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCg                                              <  2:37464/251‑1 (MQ=255)
                                                                                                                                                                                                       aaTGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTtgat                               <  1:394846/251‑1 (MQ=255)
                                                                                                                                                                                                        aTGGCAAATATCACAGCGAATGGATCTCTAACGCCA‑CCAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCACCTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACCTTGCCCGTGCACAAATCATGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGACCATGCCCGGAGTCTGCGGTtgatg                              >  2:297088/1‑251 (MQ=255)
                                                                                                                                                                                                          ggCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTtgatg                              >  2:136902/1‑249 (MQ=255)
                                                                                                                                                                                                              aataccacagacactggacttctAACGCCA‑ACAAATG‑CGGGGTGATTTCCTGGCAAATGACCGGCGATAAATCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCATAGTTTGCGAACAACCCCTGCCTGCAAGGGCAATGGCGCATCATGGCCCTTGCACAAATCTTGCTTGTCGAGTTTGAACCGGCAGATATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGTTTGATGAGCGAt                        <  2:24082/231‑1 (MQ=255)
                                                                                                                                                                                                                                  tctAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTtgct    >  1:283134/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ccaaCGCCA‑ACACATG‑CCGGTTGTTTTACGGGCATATGCCCGGCGTTAAAGCCGTCGCTGCCAACTTGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACCTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTCGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCATCCGTATCctgctgctg   <  2:169026/249‑1 (MQ=255)
                                                                                                                                                                                                                                    tAACGCCA‑ACAAATG‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTtgctgc  >  1:190612/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                             
CATTCGTGCGCAGCTGTTATGGGCCTGGGCGCGGCTGGATGAAGCCGAAGCGTCGGCGCGTAGCGGGATTGAAGTCTTGTCGTCTTATCAGCCACAGCAACAGCTTCAGTGCCTGGCAATGTTGATTCAATGCTCGCTGGCCCGTGGTGATTTAGATAACGCCCGTAGCCAGCTGAACCGTCTGGAAAACCTGCTGGGGAATGGCAAATATCACAGCGACTGGATCTCTAACGCCA‑ACAAAGT‑CCGGGTGATTTACTGGCAAATGACCGGCGATAAAGCCGCCGCTGCCAACTGGTTGCGTCATACGGCTAAACCAGAGTTTGCGAACAACCACTTCCTGCAAGGTCAATGGCGCAACATTGCCCGTGCACAAATCTTGCTGGGCGAGTTTGAACCGGCAGAAATTGTTCTCGAAGAACTCAATGAAAATGCCCGGAGTCTGCGGTTGATGAGCGATCTCAACCGTAACCTGTTGCTGC  >  NC_000913/3554806‑3555284

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