Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,085,434 C→T 6.4% R31H (CGC→CAC)  yeeE ← inner membrane protein YeeE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,085,4340CT 6.4% 292.3 / 5.8 94R31H (CGC→CAC) yeeEinner membrane protein YeeE
Reads supporting (aligned to +/- strand):  ref base C (49/38);  new base T (4/2);  total (53/41)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACTGGTGGAATGC  >  NC_000913/2085285‑2085581
                                                                                                                                                     |                                                                                                                                                   
cACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATcgc                                                                                                                                                     >  2:1460490/1‑150 (MQ=255)
cACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATcgc                                                                                                                                                     >  1:1085276/1‑150 (MQ=255)
     ataACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCAGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGaaaa                                                                                                                                                <  1:1035434/150‑1 (MQ=255)
      taACAGTACCAAGCCACGGGAACGCCCAGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAACCCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAac                                                                                                                                               <  1:215285/150‑1 (MQ=255)
       aaCAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAacc                                                                                                                                              <  2:886237/150‑1 (MQ=255)
       aaCAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAacc                                                                                                                                              >  2:1384075/1‑150 (MQ=255)
          aGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAaccacc                                                                                                                                           >  2:862028/1‑150 (MQ=255)
           ttACCAATCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGTCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTGCCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCGCGATATGTATATCGTGaaaaataatc                                                                                                                                          <  1:1178091/149‑7 (MQ=255)
           gTACCAAGCCACTGGAACGCCCCGGCTTCGTAAGTCAGTAGTCCCGCGTGAATTAACGCAAAAACCCGCACGCTTTGTACCGAAATTGCAATCAGCAATGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCt                                                                                                                                          <  2:1279229/150‑1 (MQ=255)
             aacAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATAGACATATCGCGAAACCCACCTGt                                                                                                                                        >  2:1280697/3‑150 (MQ=255)
              ccAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAACCCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                                                                                                       >  1:645116/1‑150 (MQ=255)
              ccAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                                                                                                       >  1:297754/1‑150 (MQ=255)
              ccAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                                                                                                       >  1:1353441/1‑150 (MQ=255)
              ccAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                                                                                                       >  2:68187/1‑150 (MQ=255)
               cAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCa                                                                                                                                      >  2:1096429/1‑150 (MQ=255)
                  gCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGc                                                                                                                                   >  2:1316143/1‑150 (MQ=255)
                   acACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCTCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCa                                                                                                                                  <  1:1426402/149‑1 (MQ=255)
                        ggAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAAc                                                                                                                             >  2:1301769/1‑150 (MQ=255)
                               cccGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACg                                                                                                                      <  2:1450677/150‑1 (MQ=255)
                                     ttCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCAt                                                                                                                <  1:203074/150‑1 (MQ=255)
                                      tCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATg                                                                                                               <  1:1316143/150‑1 (MQ=255)
                                       cGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGa                                                                                                              >  2:1411021/1‑150 (MQ=255)
                                          aaGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACaa                                                                                                           >  1:790327/1‑150 (MQ=255)
                                           aGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACaaa                                                                                                          >  1:1181138/1‑150 (MQ=255)
                                           aGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACaaa                                                                                                          >  2:1228481/1‑150 (MQ=255)
                                           aGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACaaa                                                                                                          >  1:539708/1‑150 (MQ=255)
                                             tCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATc                                                                                                        <  2:58332/150‑1 (MQ=255)
                                              cAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATcc                                                                                                       >  1:1270035/1‑150 (MQ=255)
                                                 tAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcac                                                                                                    <  1:991966/150‑2 (MQ=255)
                                                   ggCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTCCCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcagca                                                                                                  >  2:1490555/1‑150 (MQ=255)
                                                   ggCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATg                                                                                                               >  1:364317/1‑137 (MQ=255)
                                                   ggCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATg                                                                                                               <  2:364317/137‑1 (MQ=255)
                                                   ggCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcagca                                                                                                  >  1:834184/1‑150 (MQ=255)
                                                   ggCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCcagca                                                                                                  >  1:1287485/1‑150 (MQ=255)
                                                         ccTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGcac                                                                                            >  2:1053062/1‑150 (MQ=255)
                                                                 aaCGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGCCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAAtt                                                                                    >  2:198526/1‑150 (MQ=255)
                                                                   cGCAAAAACCCCCACGCGTTGTACCGAAATTCCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAg                                                                                  <  2:1232372/150‑1 (MQ=255)
                                                                    gCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGc                                                                                 <  1:68187/150‑1 (MQ=255)
                                                                    gCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGc                                                                                 <  1:1384075/150‑1 (MQ=255)
                                                                     cAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGcc                                                                                >  1:697180/1‑150 (MQ=255)
                                                                     cAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGcc                                                                                >  1:370708/1‑150 (MQ=255)
                                                                     cAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGcc                                                                                >  1:297716/1‑150 (MQ=255)
                                                                     cAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGcc                                                                                >  1:1481611/1‑150 (MQ=255)
                                                                           cccccACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTa                                                                          >  1:699690/1‑150 (MQ=255)
                                                                           cccccACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATAAATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTa                                                                          >  2:1185040/1‑150 (MQ=255)
                                                                                aCGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATc                                                                     <  1:1034104/150‑1 (MQ=255)
                                                                                 cGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGCGCACCACAAATTAGCCCGCTTAATATCa                                                                    >  2:1485370/1‑150 (MQ=255)
                                                                                   cTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCAtt                                                                  <  1:903750/150‑1 (MQ=255)
                                                                                     ttGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGa                                                                <  1:1365237/150‑1 (MQ=255)
                                                                                       gTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGaaa                                                              <  1:1367461/150‑1 (MQ=255)
                                                                                          ccGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACa                                                           >  2:82220/1‑150 (MQ=255)
                                                                                           cGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGGGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACAt                                                          <  2:688232/150‑1 (MQ=255)
                                                                                            gAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATa                                                         <  1:862028/150‑1 (MQ=255)
                                                                                             aaaTTGCATTCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATaa                                                        <  1:15830/150‑1 (MQ=255)
                                                                                                  gCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAaattaa                                                   >  1:1451502/1‑150 (MQ=255)
                                                                                                   cAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGTCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAaattaat                                                  <  1:1236413/150‑1 (MQ=255)
                                                                                                   cAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAaattaat                                                  >  2:1466116/1‑150 (MQ=255)
                                                                                                   cAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAaattaat                                                  >  2:1462044/1‑150 (MQ=255)
                                                                                                   cAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAaattaat                                                  >  2:1362888/1‑150 (MQ=255)
                                                                                                    aaTCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAaattaatt                                                 >  1:462208/1‑150 (MQ=255)
                                                                                                     aTCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAattaatta                                                >  2:849225/1‑150 (MQ=255)
                                                                                                      tCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAAttaattaa                                               <  1:887561/150‑1 (MQ=255)
                                                                                                        aGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAAcc                                             >  2:187747/1‑150 (MQ=255)
                                                                                                         gCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACACCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCa                                            >  2:652992/1‑150 (MQ=255)
                                                                                                          cAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAg                                           <  2:1085276/150‑1 (MQ=255)
                                                                                                          cAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAg                                           >  1:1418171/1‑150 (MQ=255)
                                                                                                           aaGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGa                                          <  1:1096429/150‑1 (MQ=255)
                                                                                                                 tAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATg                                    >  1:50923/1‑150 (MQ=255)
                                                                                                                    aaCATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTa                                 >  2:336598/1‑150 (MQ=255)
                                                                                                                     aCATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGTCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTaa                                <  2:540308/150‑1 (MQ=255)
                                                                                                                     aCATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTaa                                <  2:106979/150‑1 (MQ=255)
                                                                                                                       aTGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACTCTTCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAAGATAAATTAATTAACCAGATGAATGTTTATg                              <  1:612794/150‑1 (MQ=255)
                                                                                                                       aTGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTAAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATg                              <  1:1053062/150‑1 (MQ=255)
                                                                                                                          cGATTATTTTTCACGATATACATATCGCGATAACCACCTGTCAGGCAAAAACGCCCCCGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAGCCAGATGAATGTTAATGAgg                           >  1:755388/1‑150 (MQ=255)
                                                                                                                           aggtggttttCGCGATATGTATATCTTGAAAACCACATGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCGAGCTTAAAATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGa                          >  2:910816/7‑150 (MQ=255)
                                                                                                                            ggtggttttCGCGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAg                            >  2:1073436/6‑147 (MQ=255)
                                                                                                                            ggtggttttCGCGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAg                            <  1:1073436/142‑1 (MQ=255)
                                                                                                                              tattTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCtt                                                                           >  2:508688/1‑98 (MQ=255)
                                                                                                                              tattTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCtt                                                                           <  1:508688/98‑1 (MQ=255)
                                                                                                                               attTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATg                              >  1:1310728/1‑142 (MQ=255)
                                                                                                                               attTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATg                              <  2:1310728/142‑1 (MQ=255)
                                                                                                                                gttttCGCGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAatta                     >  2:1001069/2‑150 (MQ=255)
                                                                                                                                  tttCGCGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAattatt                   >  2:217292/1‑150 (MQ=255)
                                                                                                                                  tttCACGATATACAAATCGCGAAAACCACCTGTAAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATAACCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAattatt                   >  2:326069/1‑150 (MQ=255)
                                                                                                                                    tCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCa                 <  2:790327/150‑1 (MQ=255)
                                                                                                                                     cACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCAt                <  2:1418171/150‑1 (MQ=255)
                                                                                                                                     cACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCAt                <  1:1460490/150‑1 (MQ=255)
                                                                                                                                      aCGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATg               <  2:219715/150‑1 (MQ=255)
                                                                                                                                      aCGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATg               >  2:668995/1‑150 (MQ=255)
                                                                                                                                      aCGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATg               <  2:539708/150‑1 (MQ=255)
                                                                                                                                        gATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGAc             <  2:370708/150‑1 (MQ=255)
                                                                                                                                        gATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGAc             <  2:1481611/150‑1 (MQ=255)
                                                                                                                                             aCATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACtggtg        >  2:58105/1‑150 (MQ=255)
                                                                                                                                              cATATCGCGAAAACCACCTGTCAGGCAAAAACGCACACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACtggtgg       >  1:1179069/1‑150 (MQ=255)
                                                                                                                                                  tCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACTGGTGGAATg   <  2:50923/150‑1 (MQ=255)
                                                                                                                                                  tCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACTGGTGGAATg   >  1:1338131/1‑150 (MQ=255)
                                                                                                                                                   cgcgAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACTGGTGGAATGc  >  1:748684/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACTGGTGGAATGC  >  NC_000913/2085285‑2085581

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