Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,018,321 T→C 6.1% intergenic (‑22/‑164) ycbZ ← / → matP putative peptidase/Ter macrodomain organizer matS‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,018,3210TC 6.1% 205.2 / 3.5 66intergenic (‑22/‑164)ycbZ/matPputative peptidase/Ter macrodomain organizer matS‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (35/27);  new base C (2/2);  total (37/29)
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.55e-01

GCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTAC  >  NC_000913/1018173‑1018464
                                                                                                                                                    |                                                                                                                                               
gCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCtgtt                                                                                                                                                 <  3:140730/149‑1 (MQ=255)
        ggTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGc                                                                                                                                         <  3:415724/149‑1 (MQ=255)
              gTATCACTGAATATAGGATCGTTTTCGTCAATCAAATGTGGCTGCGCAAATATTTCGTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGAATCGTAATGGTCAAAGTTGATGTTTTTTAGTTTGTTGTCAAAGCCGCGAt                                                                                                                                   >  2:69432/1‑149 (MQ=255)
               tATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGAtt                                                                                                                                  >  4:61052/1‑149 (MQ=255)
                 tCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtat                                                                                                                                >  1:155843/1‑149 (MQ=255)
                  cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata                                                                                                                               <  4:247235/149‑1 (MQ=255)
                  cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata                                                                                                                               <  2:67522/149‑1 (MQ=255)
                  cACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATtata                                                                                                                               <  1:586343/149‑1 (MQ=255)
                    cTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAcc                                                                                                                             <  1:450241/149‑1 (MQ=255)
                     tGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCAAGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCg                                                                                                                            <  2:700989/149‑1 (MQ=255)
                       taTAAAGGATCGTTTTGGTAAATCAAATGTGTCTGACCAAATATTTCCAGATAGCTATCGGTATCAGGAAACAGGTCACGTCTTGCTAGTGTCGTAATGGTCTAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTa                                                                                                                          >  2:176829/2‑149 (MQ=255)
                       aaTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTa                                                                                                                          >  3:223905/1‑149 (MQ=255)
                       aaTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTa                                                                                                                          >  1:128259/1‑149 (MQ=255)
                       aaTAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTa                                                                                                                          >  4:189283/1‑149 (MQ=255)
                         tAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCtgttgt                                                                                                                                               >  1:302663/1‑126 (MQ=255)
                         tAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCtgttgt                                                                                                                                               <  2:302663/126‑1 (MQ=255)
                           aaGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCg                                                                                                                      >  3:28321/1‑149 (MQ=255)
                            aGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCgg                                                                                                                     <  1:361956/149‑1 (MQ=255)
                              gATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCa                                                                                                                   >  1:318862/1‑149 (MQ=255)
                              gATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCAGCGATTATACCGTAACCGGCa                                                                                                                   >  4:492682/1‑149 (MQ=255)
                               aTCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCAc                                                                                                                  >  4:335411/1‑149 (MQ=255)
                                tCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACt                                                                                                                 >  4:43170/1‑149 (MQ=255)
                                tCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACt                                                                                                                 >  2:405933/1‑149 (MQ=255)
                                  gTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTAc                                                                                                               <  3:97770/149‑1 (MQ=255)
                                  gTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTAc                                                                                                               >  4:87176/1‑149 (MQ=255)
                                   ttttCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACa                                                                                                              >  2:460262/1‑149 (MQ=255)
                                   ttttCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACa                                                                                                              <  4:28321/149‑1 (MQ=255)
                                   ttttCGGCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAGGTGGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACAGTAACCGGCACTACa                                                                                                              >  4:358106/1‑149 (MQ=255)
                                    tttCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAg                                                                                                             <  3:492682/149‑1 (MQ=255)
                                       cGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGcac                                                                                                          <  1:487640/149‑1 (MQ=255)
                                        gTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGcaca                                                                                                         >  2:103876/1‑149 (MQ=255)
                                         tcaatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATTAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTCGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCacac                                                                                                        >  4:84191/1‑149 (MQ=255)
                                         tcaatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGCATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTCGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCacac                                                                                                        >  2:242470/1‑149 (MQ=255)
                                           aatcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGt                                                                                                      >  4:269957/1‑149 (MQ=255)
                                             tcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAg                                                                                                    >  4:296893/1‑149 (MQ=255)
                                             tcaaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAg                                                                                                    >  3:387759/1‑149 (MQ=255)
                                              caaATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGa                                                                                                   >  3:590162/1‑149 (MQ=255)
                                                aaTGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGaaa                                                                                                 <  1:612630/149‑1 (MQ=255)
                                                     ggCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCAcc                                                                                            <  2:605908/149‑1 (MQ=255)
                                                          aGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACaa                                                                                       <  2:128259/149‑1 (MQ=255)
                                                          aGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACaa                                                                                       >  4:206980/1‑149 (MQ=255)
                                                               atatTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGGACACGTAGAAAGCACCGACAATACTc                                                                                  <  4:583594/149‑1 (MQ=255)
                                                                             cTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTATCGTAACGGTCAAAGTTGATGTTTTTCAGTCTGTTGTCAATGCCGCGATTATACCGCAACCGGCACTACAGCACACGTAGAACGCAGAGATAATACTCCTGGCATGGGCCtt                                                                    <  4:424254/149‑1 (MQ=255)
                                                                              tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATgg                                                                         <  2:562690/143‑1 (MQ=255)
                                                                              tATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATgg                                                                         >  1:562690/1‑143 (MQ=255)
                                                                                tCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAGTACTCCTGGCATGGGCGTTaaa                                                                 <  3:269957/149‑1 (MQ=255)
                                                                                tCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTaaa                                                                 <  4:387759/149‑1 (MQ=255)
                                                                                 cGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAg                                                                <  4:487346/149‑1 (MQ=255)
                                                                                      tCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTcac                                                           >  3:116991/1‑149 (MQ=255)
                                                                                          gAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGa                                                       >  4:486905/1‑149 (MQ=255)
                                                                                           aaCCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAcc                                                                                                                             <  3:437995/78‑1 (MQ=255)
                                                                                           aaCCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATAcc                                                                                                                             >  4:437995/1‑78 (MQ=255)
                                                                                              cAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGCCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGga                                                   >  4:633153/1‑149 (MQ=255)
                                                                                                 gTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGAt                                                >  4:478402/1‑149 (MQ=255)
                                                                                                 gTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGAt                                                >  1:667737/1‑149 (MQ=255)
                                                                                                          tGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGAttcttttctt                                       <  4:590162/149‑1 (MQ=255)
                                                                                                             aaGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCAc                                    >  2:348998/1‑149 (MQ=255)
                                                                                                             aaGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCAc                                    <  4:223905/149‑1 (MQ=255)
                                                                                                                tttCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTgg                                 <  2:20001/149‑1 (MQ=255)
                                                                                                                     tAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTaa                            >  1:166139/1‑149 (MQ=255)
                                                                                                                     tAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTaa                            >  3:121148/1‑149 (MQ=255)
                                                                                                                          gTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGc                       >  3:487211/1‑149 (MQ=255)
                                                                                                                          gTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGc                       >  1:568607/1‑149 (MQ=255)
                                                                                                                             aaaGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGa                    <  1:460262/149‑1 (MQ=255)
                                                                                                                                         ttttAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAag        <  4:428902/149‑1 (MQ=255)
                                                                                                                                          tttAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTCATATTCTGTAAaga       <  3:49198/149‑1 (MQ=255)
                                                                                                                                          tttAGTCTGTCGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAaga       <  3:84191/149‑1 (MQ=255)
                                                                                                                                          tttAGTCTGTCGTCAAAGCCGCGATTATACCGTAAACGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAaga       <  1:242470/149‑1 (MQ=255)
                                                                                                                                               tCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTac  >  4:657053/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GCAGCCGCGGTTGAGTATCACTGAATAAAGGATCGTTTTCGTCAATCAAATGTGGCTGAGCAAATATTTCCTGATAGCTATCGGTATCAGGAACCAGGTCACGCCATGCAAGTTTCGTAATGGTCAAAGTTGATGTTTTTTAGTCTGTTGTCAAAGCCGCGATTATACCGTAACCGGCACTACAGCACACGTAGAAAGCACCGACAATACTCCTGGCATGGGCGTTAAAGCTCACAGGATGGAGATTCTTTTCTTCACTGGCCTAAAAAGCTGATATTCTGTAAAGAGTTAC  >  NC_000913/1018173‑1018464

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