Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 986,050 A→T 7.1% V311V (GTT→GTA 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,0500AT 7.1% 264.2 / 6.8 93V311V (GTT→GTAompFouter membrane porin 1a (Ia;b;F)
Reads supporting (aligned to +/- strand):  ref base A (48/38);  new base T (3/4);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

TGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATTTATTAGTGA  >  NC_000913/985901‑986198
                                                                                                                                                     |                                                                                                                                                    
tGGTCAACGATACCCACAGCACCGGTGTCGGCTGAACCAACGCCCAGTTTGTTGTCAGAAGCGATCTGGTTGATGTTGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCa                                                                                                                                                      <  2:295873/150‑1 (MQ=38)
 ggTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGCTCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTGGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCaa                                                                                                                                                     >  2:222025/1‑150 (MQ=255)
 ggTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCaa                                                                                                                                                     >  1:448373/1‑150 (MQ=255)
 ggTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGATGTCAGAATCGATCTGGTTGATGATGTAGTCAAAATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCaa                                                                                                                                                     >  1:1083787/1‑150 (MQ=255)
  gTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTGTAGAAGTAGTAGGTTGCGCCAACTTCAAAGTAGTTCACCAGATCAACa                                                                                                                                                   <  1:917116/150‑1 (MQ=255)
    aaaCGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCTACAt                                                                                                                                                  >  2:1374103/1‑150 (MQ=255)
    aaaCGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCTACAt                                                                                                                                                  >  2:416254/1‑150 (MQ=255)
    aaaCGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCTACAt                                                                                                                                                  >  2:1016810/1‑150 (MQ=255)
    aaaCGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACAt                                                                                                                                                  >  2:399619/1‑150 (MQ=255)
         aTACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACAAGATCAACATCACCg                                                                                                                                             >  1:379578/1‑150 (MQ=255)
           aCCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGAt                                                                                                                                           <  1:461008/150‑1 (MQ=255)
            cccACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATa                                                                                                                                          >  1:1294746/1‑150 (MQ=255)
            cccACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATa                                                                                                                                          >  2:395596/1‑150 (MQ=255)
              cacaGCAACGGTGTCGTCTGAACCTACGCCCAGGTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTGTTTGTTGAAGTAGTAGGTTGCGCCCCCTTCAAAGTAGTTCACCAGATCAACATCCCCGAGAcc                                                                                                                                        >  1:1100572/1‑150 (MQ=255)
               acaGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCt                                                                                                                                       <  2:1515104/150‑1 (MQ=255)
                caGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCtt                                                                                                                                      <  1:1626055/150‑1 (MQ=255)
                 aGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAATATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTc                                                                                                                                     <  2:1572590/150‑1 (MQ=255)
                  tcAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCt                                                                                                                                    <  1:909107/149‑1 (MQ=255)
                  gCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTAAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCt                                                                                                                                    <  2:1497230/150‑1 (MQ=255)
                    aaCGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTAc                                                                                                                                  <  1:524081/150‑1 (MQ=255)
                      cGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGt                                                                                                                                <  2:788080/150‑1 (MQ=255)
                           tcgtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGAGCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTc                                                                                                                           <  1:1198650/150‑1 (MQ=255)
                           tagtcTGAACCTACGCCCAGTTTGTTGTCTGAATCGAGATGGTGGATGTTGTAGTAAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTGGCGCCCACATCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTc                                                                                                                           <  1:281262/148‑1 (MQ=37)
                            cgtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCg                                                                                                                          <  1:1464714/150‑1 (MQ=255)
                             ttcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGc                                                                                                                         <  1:584333/149‑1 (MQ=255)
                             gtcTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGc                                                                                                                         <  2:391864/150‑1 (MQ=255)
                                tGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCttt                                                                                                                      <  2:4689/150‑1 (MQ=255)
                                  aaCCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAg                                                                                                                    >  1:233118/1‑150 (MQ=255)
                                  aaCCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAg                                                                                                                    >  1:1156496/1‑150 (MQ=255)
                                    ccaaagCCCAGGTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACCTAGGTGGACATGTTTTTGTAGAATTAGTAGGTTTCGCCCACTTCAAAGTAGTTCACCAGATAAACATCACCTATAAATTCTACGTCTTTAGCTTTAGAt                                                                                                                  >  2:240583/6‑150 (MQ=255)
                                    ccTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTATACGTCTTTCGCTTTAGAt                                                                                                                  >  1:734275/1‑150 (MQ=255)
                                       aCGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTg                                                                                                               >  2:1129263/1‑150 (MQ=255)
                                         gCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGgt                                                                                                             <  2:387732/150‑1 (MQ=255)
                                         gCCCAGTCTGTTGTCAGAATCGATCTGGTTGATGATCTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGAGCCCACTTCAAAGTAGTTCCCCAGATCAACATCACATATACCTTCTACGTCTTTCGCTTTAGATTTGgt                                                                                                             <  1:149159/150‑1 (MQ=39)
                                          cccAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGgtg                                                                                                            <  2:857871/150‑1 (MQ=255)
                                                  gttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGAt                                                                                                    <  2:1437454/150‑1 (MQ=255)
                                                  gttgtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGAt                                                                                                    >  1:1076855/1‑150 (MQ=255)
                                                     gtCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATgga                                                                                                 >  1:926638/1‑150 (MQ=255)
                                                      tCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATggac                                                                                                >  2:487916/1‑150 (MQ=255)
                                                       cAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATggacg                                                                                               <  1:520454/150‑1 (MQ=255)
                                                        aGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATggacgg                                                                                              >  2:1429776/1‑150 (MQ=255)
                                                          aaTCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATggacggac                                                                                            >  2:821741/1‑150 (MQ=255)
                                                              gATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCTACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAg                                                                                        <  1:416254/150‑1 (MQ=255)
                                                              gATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCTACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAg                                                                                        <  1:1374103/150‑1 (MQ=255)
                                                              gATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCTACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAg                                                                                        <  1:1016810/150‑1 (MQ=255)
                                                              gATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAg                                                                                        >  2:992800/1‑150 (MQ=255)
                                                                 cTGGTTGATGATGTAGTCAACATAGGTCGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTGCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGAcc                                                                                     <  2:158700/150‑1 (MQ=255)
                                                                  tGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCg                                                                                    <  2:1294746/150‑1 (MQ=255)
                                                                      tgatgaTGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAAt                                                                                <  2:511250/150‑1 (MQ=255)
                                                                      tgatgaTGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAAt                                                                                <  2:448373/150‑1 (MQ=255)
                                                                       tatgatgTAGTCAACATAGGTGGCCCTGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCCCCGATTCCTTCTCCGTCTTTCGCTTAAGATTTGGTGTAAGCGATCGACGGACGCCGACCGAAATc                                                                               >  1:190084/2‑150 (MQ=38)
                                                                       gatgatGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATc                                                                               <  2:1662398/150‑1 (MQ=255)
                                                                              tAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTTCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATCCCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTgg                                                                        <  2:25557/150‑1 (MQ=255)
                                                                              tAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTgg                                                                        >  2:1346914/1‑150 (MQ=255)
                                                                               aGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCttt                                                                                                                      >  2:118504/1‑103 (MQ=255)
                                                                               aGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCttt                                                                                                                      <  1:118504/103‑1 (MQ=255)
                                                                               aGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTggg                                                                       >  2:667487/1‑149 (MQ=255)
                                                                               aGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGt                                                                       >  1:381282/1‑150 (MQ=255)
                                                                                  cAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTAtt                                                                    >  1:1369294/1‑150 (MQ=255)
                                                                                    aCATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTgc                                                                  >  2:1260285/1‑150 (MQ=255)
                                                                                    aCATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTgc                                                                  >  2:1707712/1‑150 (MQ=255)
                                                                                             gACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACAAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAAc                                                         >  1:734191/1‑150 (MQ=255)
                                                                                              aCATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCATCATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACa                                                        >  1:32638/1‑150 (MQ=255)
                                                                                              aCATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACa                                                        >  2:949938/1‑150 (MQ=255)
                                                                                                  gTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAAc                                                    >  2:33327/1‑150 (MQ=255)
                                                                                                  gTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAAc                                                    >  1:1004212/1‑150 (MQ=255)
                                                                                                   tttttGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACg                                                   <  1:1322219/150‑1 (MQ=255)
                                                                                                       tgtttAAGTAGTACGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCTCCGATACCTTCTACGACTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGCCCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCtt                                               <  2:267104/150‑1 (MQ=255)
                                                                                                       tgttgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGTCCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCtt                                               >  2:1302333/1‑150 (MQ=255)
                                                                                                        gttgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTg                                              >  1:1350648/1‑150 (MQ=255)
                                                                                                         ttgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTATCAGAAAGTCTTgg                                             >  2:1097683/1‑149 (MQ=255)
                                                                                                          tgAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCg                                            <  1:1506658/150‑1 (MQ=255)
                                                                                                           gAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGt                                           >  1:1458558/1‑150 (MQ=255)
                                                                                                                  tagGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTttgttg                                    <  1:1293873/150‑1 (MQ=255)
                                                                                                                   agGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAAAATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTgg                                   <  2:1163253/150‑1 (MQ=255)
                                                                                                                     gTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAATTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCg                                 >  2:1451489/1‑150 (MQ=255)
                                                                                                                         cgcCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAATCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGc                             <  2:734275/150‑1 (MQ=255)
                                                                                                                            ccACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGc                          >  2:1319341/1‑150 (MQ=255)
                                                                                                                               cTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTgg                       >  1:4812/1‑150 (MQ=255)
                                                                                                                               cTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTgg                       >  2:1316663/1‑150 (MQ=255)
                                                                                                                                  cAAAGTAGTTCACCAGATCAACATCCCCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCATGCGTTTTGTTGGCGAAGCCGCTGgtgt                    >  2:1067924/1‑150 (MQ=255)
                                                                                                                                  cAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGATATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGgtgt                    <  2:1076855/150‑1 (MQ=255)
                                                                                                                                  cAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGgtgt                    >  2:641433/1‑150 (MQ=255)
                                                                                                                                   aaaGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGtt                   >  2:1376814/1‑150 (MQ=255)
                                                                                                                                     agtagtTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTg                 <  1:23384/150‑1 (MQ=255)
                                                                                                                                       tagtTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTa               >  2:177243/1‑150 (MQ=255)
                                                                                                                                        agtTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTaa              >  2:1152111/1‑150 (MQ=255)
                                                                                                                                         gtTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTaaa             <  1:1316663/150‑1 (MQ=255)
                                                                                                                                          ttCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAAt            <  1:399619/150‑1 (MQ=255)
                                                                                                                                           tCACCAGATCTACATCACCCATACCTTCTACGTCTTTCGCTTTAGATTTGTTCTACGCGATGGACGGATGCAGAGCGAAATCGATCTGGTATTGTTTATCTAACAGAACGTCTTGGGTTTTGTTGGTGAAGCCGCTGGTGTTTGTAAttt           <  1:234663/150‑1 (MQ=255)
                                                                                                                                           tCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGTGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAAtt           >  1:1431323/1‑150 (MQ=255)
                                                                                                                                           tCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAAtt           >  2:539145/1‑150 (MQ=255)
                                                                                                                                            cACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAAttt          >  2:653549/1‑150 (MQ=255)
                                                                                                                                               cAGATCAACCTCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATttatt       <  2:623227/150‑1 (MQ=255)
                                                                                                                                                aGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATttatta      >  2:1620891/1‑150 (MQ=255)
                                                                                                                                                 gATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATttattat     >  2:760291/1‑149 (MQ=255)
                                                                                                                                                    cAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATTTATTAGTGa  >  2:1204409/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGTAAACGATACCCACAGCAACGGTGTCGTCTGAACCTACGCCCAGTTTGTTGTCAGAATCGATCTGGTTGATGATGTAGTCAACATAGGTGGACATGTTTTTGTTGAAGTAGTAGGTTGCGCCCACTTCAAAGTAGTTCACCAGATCAACATCACCGATACCTTCTACGTCTTTCGCTTTAGATTTGGTGTAAGCGATGGACGGACGCAGACCGAAATCGAACTGGTATTGCGCAACTAACAGAACGTCTTGCGTTTTGTTGGCGAAGCCGCTGGTGTTTGTAAATTTATTAGTGA  >  NC_000913/985901‑986198

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: