Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 986,411 G→T 6.0% S191Y (TCT→TAT) ‡ 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,4110GT 6.0% 260.5 / 3.1 82S191Y (TCT→TAT) ‡ompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base T (3/2);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

CAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCA  >  NC_000913/986263‑986560
                                                                                                                                                    |                                                                                                                                                     
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:124910/150‑1 (MQ=255)
                   tCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGc                                                                                                                                   >  2:4590/1‑150 (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                                                                                              <  2:213969/150‑1 (MQ=255)
                                                        gCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGc                                                                                              <  1:499500/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                                                                       >  1:389401/1‑150 (MQ=255)
                                                                               ccATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTg                                                                       <  2:70978/150‑1 (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               <  2:717273/150‑1 (MQ=255)
                                                                                                                                       acCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTa               <  1:119067/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)
                                                                                                                                                    gAACGGCGTGCAGTGTCATGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATcacca  >  2:407484/1‑150 (MQ=255)
                                                                                                                                                    gAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGCCCAACGAAGAAGTCATCGCTGTATGCATTATcacca  >  1:1008900/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGACCAGTAGCCCACTGTTCAGCTTTTTTACCGTTGCCAAGAGGTTGAGCTTCTTGCAGGTTGGTACGGTCAGCTGCACCATAAGCACCAACGATACCAAAGCCTTCGTATTCGTAGCTGATAGAACCGCCAACACCGTCGCCGTTAGAACGGCGTGCAGTGTCACGCTCGTTTTTACCCAGGTACTGAACAGCGAAGTTCAGGCCATCAACCAGACCAAAGAAGTTGGAGTTACGATAGGTAGCAACGCCGCCAACACGACCAACGAAGAAGTCATCGCTGTATGCAGTATCACCA  >  NC_000913/986263‑986560

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: