Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 986,410 A→C 6.3% S191S (TCT→TCG) ‡ ompF ← outer membrane porin 1a (Ia;b;F)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913986,4100AC 6.3% 221.4 / 6.6 80S191S (TCT→TCG) ‡ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (39/36);  new base C (3/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.97e-01

TTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACC  >  NC_000913/986261‑986559
                                                                                                                                                     |                                                                                                                                                     
ttCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTa                                                                                                                                                       <  2:735720/150‑1 (MQ=255)
  cAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGa                                                                                                                                                     >  1:990907/1‑150 (MQ=255)
   aGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATAACAAAGCCTTCGTATTCGTAGCTTATAGAACCGCCAACACCGTCGACGTTACaa                                                                                                                                                    >  2:18386/1‑150 (MQ=255)
        aGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTCTAACGGCg                                                                                                                                               <  1:669930/150‑1 (MQ=255)
         gTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGACGCCGTTAGAACGGCGt                                                                                                                                              >  2:157232/1‑150 (MQ=255)
            gCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTGCGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCa                                                                                                                                           >  2:809915/1‑150 (MQ=255)
             cccACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAg                                                                                                                                          <  1:1356059/150‑1 (MQ=255)
               cACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAgtg                                                                                                                                        <  1:1070268/150‑1 (MQ=255)
                   gTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTccc                                                                                                                                    >  2:716148/1‑148 (MQ=255)
                     tCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGc                                                                                                                                  >  2:4590/1‑150 (MQ=255)
                     tCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGc                                                                                                                                  <  2:124910/150‑1 (MQ=255)
                       aGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTc                                                                                                                                >  2:673372/1‑150 (MQ=255)
                          ttttttACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGtt                                                                                                                             <  1:374089/150‑1 (MQ=255)
                          ttttttACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATAACAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGtt                                                                                                                             >  1:535784/1‑150 (MQ=255)
                           tttttACCGTTGCCAAGAGGTTGAGCTTCTGGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGttt                                                                                                                            <  2:1310261/150‑1 (MQ=255)
                            ttttACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTTCGGTCAGCTGCCCCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGtttt                                                                                                                           >  1:84606/1‑150 (MQ=255)
                                aCCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTAcc                                                                                                                       <  2:1391336/150‑1 (MQ=255)
                                       ccAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTAc                                                                                                                >  1:1303698/1‑150 (MQ=255)
                                          agagGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGATCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGa                                                                                                             >  1:837510/1‑150 (MQ=255)
                                             ggTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACa                                                                                                          <  1:1328125/150‑1 (MQ=255)
                                                       cttGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCa                                                                                                <  2:688876/150‑1 (MQ=255)
                                                         tGCAGGTTGGTACGGTCAGCTGCAGCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATTGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACACAGGTACTGAACAGCGTAGTTCAgg                                                                                              >  2:261550/1‑150 (MQ=255)
                                                         tGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCTACGATACCAAAGCGGTCGTATTCGGATCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCAAGCACGTTTTTACACAGGTACTGAACAGCGAAGTTCAgg                                                                                              >  1:159321/1‑150 (MQ=255)
                                                          gCAGGTTGGTACGGTCAGCTGCACCATAAGCCCCAACGATGCCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGACCGGCGTGCAGTGTCCCGCTCGTTTTTACCCACGTACTGAACAGCGAAGTTCAGGc                                                                                             >  2:244341/1‑150 (MQ=255)
                                                          gCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                                                                             <  1:499500/150‑1 (MQ=255)
                                                          gCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                                                                             <  2:213969/150‑1 (MQ=255)
                                                          gCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGACAACACCGTCGCCGTAAGAACGGCGTGCAGTTACACGCTCGTTTTTACCCAAGTACTGAACAGCGAAGTTGAGGc                                                                                             >  2:890441/1‑150 (MQ=255)
                                                           cAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGcc                                                                                            >  1:1331353/1‑150 (MQ=255)
                                                            aGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACTTCGCTGTTCAgtacc                                                                                           <  2:1296292/150‑5 (MQ=255)
                                                            aGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCa                                                                                           >  1:1652227/1‑150 (MQ=255)
                                                            aGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCa                                                                                           >  1:1590705/1‑150 (MQ=255)
                                                              gTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACTTCGCTGTTCAgtacctg                                                                                         <  2:1303698/150‑7 (MQ=255)
                                                               ttGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACTTCGCTGTTCAgtacctgg                                                                                        <  2:135985/150‑8 (MQ=255)
                                                               ttGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCa                                                                                        <  1:1534511/150‑1 (MQ=255)
                                                                  gTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTTGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAcc                                                                                     <  1:4748/150‑1 (MQ=255)
                                                                   aaCGGGCAGCTGCACCCTAAGCACCACAGATACCAACGACTTCGTATTCGTAGCTGCTAGACCCGCCAACACCGTCGCCGTTAGACCGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCa                                                                                    <  2:330137/149‑1 (MQ=255)
                                                                     cGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGa                                                                                  >  1:70978/1‑150 (MQ=255)
                                                                      ggTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGAc                                                                                 <  1:590098/150‑1 (MQ=255)
                                                                      ggTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCACCAACCAGAc                                                                                 <  2:1327758/150‑1 (MQ=255)
                                                                         ccgCTGCACCATAAGCAGCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACAGCGTCGCCGTTAGTACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCaa                                                                              <  2:1149242/148‑1 (MQ=255)
                                                                                aCCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCATCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAAACAGACCAAAGAAGtt                                                                       >  1:663310/1‑150 (MQ=255)
                                                                                 ccATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                                                      <  2:70978/150‑1 (MQ=255)
                                                                                 ccATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                                                      >  1:389401/1‑150 (MQ=255)
                                                                                 ccATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGCCCATCACACAGACCAAAGAAGATg                                                                      >  1:356207/1‑150 (MQ=255)
                                                                                  cATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCAGTTAGAACTGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTgg                                                                     <  1:157232/150‑1 (MQ=255)
                                                                                   aTAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGa                                                                    <  2:1192476/150‑1 (MQ=255)
                                                                                     aaGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGt                                                                  <  1:1108670/150‑1 (MQ=255)
                                                                                      aGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGtt                                                                 <  1:419975/150‑1 (MQ=255)
                                                                                       gCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTa                                                                >  1:1610402/1‑150 (MQ=255)
                                                                                       gCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTa                                                                <  1:492662/150‑1 (MQ=255)
                                                                                         aCCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACg                                                              <  1:673372/150‑1 (MQ=255)
                                                                                              cGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAgg                                                         <  1:1402129/150‑1 (MQ=255)
                                                                                                   ccAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCATCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCa                                                    >  2:975473/1‑150 (MQ=255)
                                                                                                     aaaGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGATTTACGATAGGTAGCAAc                                                  >  2:1604072/1‑150 (MQ=255)
                                                                                                      aaGCCTTCGTATTCGTAGCTGATATAACCTCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcg                                                 <  2:535784/150‑1 (MQ=255)
                                                                                                       aGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgc                                                <  1:809915/150‑1 (MQ=255)
                                                                                                          cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgccgc                                             >  1:1557126/1‑150 (MQ=255)
                                                                                                          cTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAAcgccgc                                             >  2:119067/1‑150 (MQ=255)
                                                                                                           ttcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACgccgcc                                            <  1:1515218/150‑1 (MQ=255)
                                                                                                            tcgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCa                                           <  2:84606/150‑1 (MQ=255)
                                                                                                             cgtattcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCaa                                          >  2:1142300/1‑150 (MQ=255)
                                                                                                                tttcgtaGCTGATAGAACCGCCACCACCGGCGCCGGTAGAACGGCGTGCAGTGGCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGACGTTGGAGTTACGATAGGTAGCAACGCCGCCAacac                                       <  2:305868/149‑1 (MQ=38)
                                                                                                                attcgtaGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAacac                                       >  1:717273/1‑150 (MQ=255)
                                                                                                                           aTAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACgaag                            <  1:18386/150‑1 (MQ=255)
                                                                                                                             ataaCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTAATGAACAGCGAAGTTAATTCCATCAAACAGACCAAAGAAGTTTGAGAGACGATAGGTAGCAACTCCGCCAACACAACCAACgaagca                          >  1:1492980/3‑148 (MQ=255)
                                                                                                                             aGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACgaagaa                          >  2:541514/1‑150 (MQ=255)
                                                                                                                                aCCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGGCATCAACCAGCCCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCAGCCAACAGGACCAACGAAGAAGTc                       >  2:1012752/1‑150 (MQ=255)
                                                                                                                                    ccAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCg                   >  1:206771/1‑150 (MQ=255)
                                                                                                                                    ccAACACCGTCGCCGTTAGAACGGCGTGCAATGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCg                   >  1:1311059/1‑150 (MQ=255)
                                                                                                                                     cAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGt                        >  2:149640/1‑144 (MQ=255)
                                                                                                                                     cAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGt                        <  1:149640/144‑1 (MQ=255)
                                                                                                                                       cactgcacGCCGTTCTAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTg                >  1:83253/8‑150 (MQ=255)
                                                                                                                                        cacCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCCGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACCCCGCCCACACGACCAACGAAGAAGTCATCGCTGt               >  1:1055889/1‑150 (MQ=255)
                                                                                                                                         acCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTa              <  1:119067/150‑1 (MQ=255)
                                                                                                                                         acCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTa              <  2:717273/150‑1 (MQ=255)
                                                                                                                                           gcacGCCGTTCTAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATg            >  2:702059/4‑150 (MQ=255)
                                                                                                                                           gcacGCCGTTCTAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCAGTATg            >  1:267796/4‑150 (MQ=255)
                                                                                                                                              cGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAg         >  1:1562906/1‑150 (MQ=255)
                                                                                                                                               gCCGTTCTAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGt        <  2:83253/150‑1 (MQ=255)
                                                                                                                                                 cGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGATCGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGAATGCAGTAt      >  1:190446/1‑150 (MQ=255)
                                                                                                                                                  gTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATc     >  1:233646/1‑150 (MQ=255)
                                                                                                                                                     aGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACTCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATcacc  >  1:1018440/1‑150 (MQ=255)
                                                                                                                                                     aGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATcacc  <  2:389401/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACC  >  NC_000913/986261‑986559

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