Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,015,333 A→T 7.1% intergenic (+19/+21) metE → / ← ysgA cobalamin‑independent homocysteine transmethylase/putative dienelactone hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,015,3330AT 7.1% 352.1 / 11.6 113intergenic (+19/+21)metE/ysgAcobalamin‑independent homocysteine transmethylase/putative dienelactone hydrolase
Reads supporting (aligned to +/- strand):  ref base A (50/55);  new base T (3/5);  total (53/60)
Fisher's exact test for biased strand distribution p-value = 7.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 116 total aligned reads displayed.

CGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACA  >  NC_000913/4015184‑4015481
                                                                                                                                                     |                                                                                                                                                    
cGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTaa                                                                                                                                                      >  2:118177/1‑150 (MQ=255)
cGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTaa                                                                                                                                                      >  1:1293637/1‑150 (MQ=255)
 ggCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAAt                                                                                                                                                     <  2:1266175/150‑1 (MQ=255)
 ggCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAAt                                                                                                                                                     <  2:1391674/150‑1 (MQ=255)
   cAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATac                                                                                                                                                   >  1:51693/1‑150 (MQ=255)
   cAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATac                                                                                                                                                   <  1:177620/150‑1 (MQ=255)
    aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacc                                                                                                                                                  >  1:1626548/1‑150 (MQ=255)
    aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacc                                                                                                                                                  >  1:142261/1‑150 (MQ=255)
    aaaaCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATacc                                                                                                                                                  >  1:215425/1‑150 (MQ=255)
       aCGCATTCCGGCAGAGCGCCTGTGGGTCAAGCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACACGCGCGGCACTGGTGAACATGGTGCAGTCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATaccacc                                                                                                                                               <  2:775594/150‑1 (MQ=255)
       aCGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATaccacc                                                                                                                                               <  2:664024/150‑1 (MQ=255)
        cGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAATACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGATCATGGTGCAGTCGGCGCATAACTTGCGTCGGTGGTAAAATCCAAACCGGGTGGTAATACCAccc                                                                                                                                              <  2:1612368/150‑1 (MQ=255)
          cATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                                                                            <  1:1211918/150‑1 (MQ=255)
          cATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCgg                                                                                                                                            <  1:1478920/150‑1 (MQ=255)
          cATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCgg                                                                                                                                            <  1:664067/150‑1 (MQ=255)
           aTTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGt                                                                                                                                           <  1:1541287/150‑1 (MQ=255)
            ttCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTc                                                                                                                                          <  2:150894/150‑1 (MQ=255)
            ttCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTc                                                                                                                                          <  2:243751/150‑1 (MQ=255)
              ccGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttg                                                                                                                                        <  2:1303085/150‑4 (MQ=255)
              ccGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTATTACCACCCGGtttg                                                                                                                                        <  1:1305116/150‑4 (MQ=255)
              ccGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCtt                                                                                                                                        >  1:737075/1‑150 (MQ=255)
                 gCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTtc                                                                                                                                     <  1:1235833/150‑1 (MQ=255)
                    gagCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCa                                                                                                                                  <  1:701979/150‑1 (MQ=255)
                      gcgcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCAtt                                                                                                                                <  2:1152945/150‑1 (MQ=255)
                       cgcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTa                                                                                                                               >  1:1182551/1‑150 (MQ=255)
                        gcCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTAc                                                                                                                              <  2:1341726/150‑1 (MQ=255)
                         ccTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACa                                                                                                                             >  1:116037/1‑150 (MQ=255)
                             tgGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGa                                                                                                                         >  1:546878/1‑150 (MQ=255)
                               ggTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGAct                                                                                                                       >  2:1447098/1‑150 (MQ=255)
                                gTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActt                                                                                                                      >  2:379805/1‑150 (MQ=255)
                                 tCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActtc                                                                                                                     <  1:258726/150‑1 (MQ=255)
                                   aaCCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGActtctt                                                                                                                   >  1:1051614/1‑150 (MQ=255)
                                       cGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCcca                                                                                                               >  1:1155951/1‑150 (MQ=255)
                                         gACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCccacc                                                                                                             <  1:150989/150‑1 (MQ=255)
                                               ggCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                                                                       >  1:206045/1‑150 (MQ=255)
                                               ggCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTg                                                                                                       >  1:1013829/1‑150 (MQ=255)
                                                gCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGc                                                                                                      >  2:1398274/1‑150 (MQ=255)
                                                gCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGc                                                                                                      >  2:1582380/1‑150 (MQ=255)
                                                 ccTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCt                                                                                                     <  2:1626548/150‑1 (MQ=255)
                                                   tGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTa                                                                                                   <  1:108977/150‑1 (MQ=255)
                                                    gAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaa                                                                                                  >  1:1531988/1‑150 (MQ=255)
                                                     aaaaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTaaa                                                                                                 <  1:1610132/150‑1 (MQ=255)
                                                       aaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAcc                                                                                               >  2:61284/1‑150 (MQ=255)
                                                       aaCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGGCGGGGGTAAAATCCAAACCGGGGGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAAcc                                                                                               <  2:1104721/150‑1 (MQ=255)
                                                        aCGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGTGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCa                                                                                              >  2:539850/1‑150 (MQ=255)
                                                           cgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTc                                                                                           <  1:220088/150‑1 (MQ=255)
                                                           cgcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTc                                                                                           <  1:1398274/150‑1 (MQ=255)
                                                            gcgGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTcc                                                                                          <  1:1096283/150‑1 (MQ=255)
                                                              ggCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAg                                                                                        <  2:487818/150‑1 (MQ=255)
                                                               gCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGc                                                                                       <  2:51693/150‑1 (MQ=255)
                                                               gCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCAGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGc                                                                                       <  2:969410/150‑1 (MQ=255)
                                                                   gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATAc                                                                                   >  2:1042598/1‑150 (MQ=255)
                                                                   gCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATAc                                                                                   >  2:674928/1‑150 (MQ=255)
                                                                       gAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTg                                                                               >  2:96757/1‑150 (MQ=255)
                                                                        aaaCCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgc                                                                              <  2:1477098/150‑1 (MQ=255)
                                                                         aaCCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgcc                                                                             <  2:1531988/150‑1 (MQ=255)
                                                                           cccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccag                                                                           <  2:1013829/150‑1 (MQ=255)
                                                                            ccGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagc                                                                          <  2:1500894/150‑1 (MQ=255)
                                                                             cgcgcgGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTgccagcc                                                                         <  2:1474797/150‑1 (MQ=255)
                                                                                          cGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGAt                                                            >  2:1002590/1‑150 (MQ=255)
                                                                                            aaCATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGc                                                          <  1:1451302/150‑1 (MQ=255)
                                                                                             aCATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCg                                                         <  1:1518249/150‑1 (MQ=255)
                                                                                                tGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCa                                                      >  2:761273/1‑150 (MQ=255)
                                                                                                tGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCa                                                      >  2:1020218/1‑150 (MQ=255)
                                                                                                 cgTGCAGGCGGCGCAGAACTTGCGGCGGGGGTAAAATCCAAACCGGGGGGTAATACCACCCGGTCTTTTCTCATTACAGGGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCAt                                                     <  2:848240/149‑1 (MQ=255)
                                                                                                     cAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATa                                                 <  2:664334/150‑1 (MQ=255)
                                                                                                       ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGc                                                                                 <  2:137507/116‑1 (MQ=255)
                                                                                                       ggcggcGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGc                                                                                 >  1:137507/1‑116 (MQ=255)
                                                                                                        gcggcgCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCt                                              >  1:1234938/1‑150 (MQ=255)
                                                                                                           gcgcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATATCTCgg                                           >  1:471653/1‑150 (MQ=255)
                                                                                                            ggcAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCggg                                          <  2:922353/149‑1 (MQ=255)
                                                                                                              cAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCg                                        >  2:421999/1‑150 (MQ=255)
                                                                                                              cAGAACTTGCGGCGGGTGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATCCTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCg                                        <  1:623530/150‑1 (MQ=255)
                                                                                                                     gtcGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCa                                 >  2:1199387/3‑150 (MQ=255)
                                                                                                                        gTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCg                              >  1:698833/1‑150 (MQ=255)
                                                                                                                           gggggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTg                           <  2:546878/150‑1 (MQ=255)
                                                                                                                            ggggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGa                          >  2:431266/1‑150 (MQ=255)
                                                                                                                              ggTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAat                        <  2:1523257/150‑1 (MQ=255)
                                                                                                                               gTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatg                       <  2:24284/150‑1 (MQ=255)
                                                                                                                               gTAAAATCCAAACCGGGTCGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCACACCATACTGCTTAAACCATTCCAGCAAACGCTGCCAGCCATCATCTGCAGATGAGGCATGATAGCTCGGGCGATAATCAGCGTTCAatg                       >  2:1426573/1‑150 (MQ=255)
                                                                                                                                  aaaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgcat                    <  2:206045/150‑1 (MQ=255)
                                                                                                                                   agaaaagaCCGGGTGGTATTACCACCCGGTCTTTTGTGATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgcatg                   >  2:631385/8‑150 (MQ=255)
                                                                                                                                   agaaaagaCCGGGTGGTATTACCACCCCGTCTTTTCTCATTACAGCGACATCTTCCCACCAAACTCCTTAAACCATTACACCATACGCTGCCAGCCAACATCTGCAGAAGCGGCATGAAAGCTCGGGCGATAATCAGCGTTGAaagcatg                   >  2:636722/8‑150 (MQ=255)
                                                                                                                                   aaTCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAatgcatg                   >  2:1320566/1‑150 (MQ=255)
                                                                                                                                     tCCAAACCGGGTGGTAATACCCCCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGcc                 >  2:170567/1‑150 (MQ=255)
                                                                                                                                     tCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGcc                 >  2:120370/1‑150 (MQ=255)
                                                                                                                                      ccAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGccc                >  2:117906/1‑150 (MQ=255)
                                                                                                                                      ccAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGccc                <  1:1089357/150‑1 (MQ=255)
                                                                                                                                       cAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCg               >  1:369222/1‑150 (MQ=255)
                                                                                                                                        aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCgg              <  1:985626/150‑1 (MQ=255)
                                                                                                                                        aaaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCgg              <  2:1246932/150‑1 (MQ=255)
                                                                                                                                         gaCCGGGTGGTATTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGc             >  1:1199178/2‑150 (MQ=255)
                                                                                                                                         aaCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGc             >  2:329121/1‑150 (MQ=255)
                                                                                                                                          aCCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCg            <  1:1320566/150‑1 (MQ=255)
                                                                                                                                            cGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTc          >  1:618578/1‑150 (MQ=255)
                                                                                                                                              ggtggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCg        >  2:504770/1‑150 (MQ=255)
                                                                                                                                                 ggtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGt     >  1:79482/1‑150 (MQ=255)
                                                                                                                                                  gtAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTa    <  1:539850/150‑1 (MQ=255)
                                                                                                                                                   tAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTac   <  1:1199387/150‑1 (MQ=255)
                                                                                                                                                    aaTACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTaca  >  2:241666/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGCAAAACGCATTCCGGCAGAGCGCCTGTGGGTCAACCCGGACTGTGGCCTGAAAACGCGCGGCTGGCCAGAAACCCGCGCGGCACTGGCGAACATGGTGCAGGCGGCGCAGAACTTGCGTCGGGGGTAAAATCCAAACCGGGTGGTAATACCACCCGGTCTTTTCTCATTACAGCGACTTCTTCCCACCATACTGCTTAAACCATTCCAGCATACGCTGCCAGCCATCTTCTGCAGATGCGGCATGATAGCTCGGGCGATAATCAGCGTTGAATGCATGCCCGGCGTCCGGGTACA  >  NC_000913/4015184‑4015481

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