Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,085,426 A→G 6.4% Y34H (TAT→CAT)  yeeE ← inner membrane protein YeeE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,085,4260AG 6.4% 283.9 / 7.4 94Y34H (TAT→CAT) yeeEinner membrane protein YeeE
Reads supporting (aligned to +/- strand):  ref base A (49/39);  new base G (4/2);  total (53/41)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

TATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACTGGTG  >  NC_000913/2085279‑2085575
                                                                                                                                                   |                                                                                                                                                     
tataCCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAATCCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACa                                                                                                                                                     <  2:865671/150‑1 (MQ=255)
    cccACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATc                                                                                                                                                 >  1:688232/1‑150 (MQ=255)
      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                                                                                                                                 >  2:68187/1‑150 (MQ=255)
                    ccAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                                                                                                 >  1:297754/1‑150 (MQ=255)
                    ccAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTc                                                                                                                                 >  1:1353441/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:1481611/1‑150 (MQ=255)
                                                                           cAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGcc                                                                          >  1:297716/1‑150 (MQ=255)
                                                                           cAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGcc                                                                          >  1:370708/1‑150 (MQ=255)
                                                                           cAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGcc                                                                          >  1:697180/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                                     >  1:1418171/1‑150 (MQ=255)
                                                                                                                cAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAg                                     <  2:1085276/150‑1 (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                      <  1:1073436/142‑1 (MQ=255)
                                                                                                                                  ggtggttttCGCGATATGTATATCGTGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAg                      >  2:1073436/6‑147 (MQ=255)
                                                                                                                                    tattTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCtt                                                                     >  2:508688/1‑98 (MQ=255)
                                                                                                                                    tattTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCtt                                                                     <  1:508688/98‑1 (MQ=255)
                                                                                                                                     attTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATg                        <  2:1310728/142‑1 (MQ=255)
                                                                                                                                     attTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATg                        >  1:1310728/1‑142 (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          <  1:1460490/150‑1 (MQ=255)
                                                                                                                                           cACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCAt          <  2:1418171/150‑1 (MQ=255)
                                                                                                                                            aCGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATg         >  2:668995/1‑150 (MQ=255)
                                                                                                                                            aCGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATg         <  2:219715/150‑1 (MQ=255)
                                                                                                                                            aCGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATg         <  2:539708/150‑1 (MQ=255)
                                                                                                                                              gATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGAc       <  2:1481611/150‑1 (MQ=255)
                                                                                                                                              gATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGAc       <  2:370708/150‑1 (MQ=255)
                                                                                                                                                   aCATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACtggtg  >  2:58105/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TATACCCACCTATAACAGTACCAAGCCACGGGAACGCCCCGGCTTCGTAAGTCAGTAGGCCCGCCTGAATTAACGCAAAAACCCCCACGCTTTGTACCGAAATTGCAATCAGCAAGGCGTAAAACATGCGATTATTTTTCACGATATACATATCGCGAAAACCACCTGTCAGGCAAAAACGCCCACGCTGCATGACAAATCCCAGCAGAGCACCACAAATTAGCCCGCTTAATATCATTGAAAACATAAATTAATTAACCAGATGAATGTTAATGAGGAAATTATTCATGACTGGTG  >  NC_000913/2085279‑2085575

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: