Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,315,438 A→T 6.2% F34Y (TTT→TAT) ‡ phnO ← aminoalkylphosphonate N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,315,4380AT 6.2% 300.3 / 4.8 95F34Y (TTT→TAT) ‡phnOaminoalkylphosphonate N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (48/41);  new base T (2/4);  total (50/45)
Fisher's exact test for biased strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTG  >  NC_000913/4315289‑4315587
                                                                                                                                                     |                                                                                                                                                     
ccaccaACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCTAATGCAGGCCGATCATGCCGACAACTTCGCCATCTAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGaa                                                                                                                                                       <  1:1242122/150‑1 (MQ=255)
  accaACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAc                                                                                                                                                     <  1:1229871/150‑1 (MQ=255)
  accaACTCCTGAATTTCGCCGATCCAGTTGACATGACGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAc                                                                                                                                                     <  1:1467135/150‑1 (MQ=255)
     aaCTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAcgcg                                                                                                                                                  <  2:744541/150‑1 (MQ=255)
      aCTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCgt                                                                                                                                                 <  2:558814/150‑1 (MQ=255)
          cTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggt                                                                                                                                             <  1:1379428/150‑3 (MQ=255)
          cTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggt                                                                                                                                             <  2:879522/150‑3 (MQ=255)
           tGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggtt                                                                                                                                            <  1:263896/150‑4 (MQ=255)
            gAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAACCCCACGCGAAACGCGAgggggt                                                                                                                                           >  2:909203/1‑150 (MQ=255)
              aTTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGTTTCGCgtgggtttta                                                                                                                                         <  1:1270157/150‑7 (MQ=255)
                 tCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGGGGTCAAAc                                                                                                                                      >  1:255321/1‑150 (MQ=255)
                  cGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAAca                                                                                                                                     >  2:512257/1‑149 (MQ=255)
                  cGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACt                                                                                                                                     >  2:1549505/1‑150 (MQ=255)
                   gCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTc                                                                                                                                    >  1:1264038/1‑150 (MQ=255)
                   gCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTc                                                                                                                                    <  1:103771/150‑1 (MQ=255)
                    ccGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTcc                                                                                                                                   <  1:53567/150‑1 (MQ=255)
                      gATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGc                                                                                                                                 >  1:1255617/1‑150 (MQ=255)
                      gATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGc                                                                                                                                 >  1:523474/1‑150 (MQ=255)
                           aGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGtt                                                                                                                            >  2:1609890/1‑150 (MQ=255)
                           aGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGtt                                                                                                                            >  2:260989/1‑150 (MQ=255)
                           aGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGtt                                                                                                                            >  2:1404821/1‑150 (MQ=255)
                              tGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTa                                                                                                                         <  1:815596/150‑1 (MQ=255)
                               gACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTg                                                                                                                              >  1:159165/1‑144 (MQ=255)
                               gACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTg                                                                                                                              <  2:159165/144‑1 (MQ=255)
                               gACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAg                                                                                                                        >  1:1038995/1‑150 (MQ=255)
                               gACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAg                                                                                                                        >  1:1004710/1‑150 (MQ=255)
                                    gatgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTcac                                                                                                                   >  1:290341/1‑150 (MQ=255)
                                     atgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTcaca                                                                                                                  <  2:390313/150‑1 (MQ=255)
                                      tgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCGTCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACaa                                                                                                                 >  1:520069/1‑150 (MQ=255)
                                      tgCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACaa                                                                                                                 >  1:374892/1‑150 (MQ=255)
                                       gCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACaaa                                                                                                                >  2:1331693/1‑150 (MQ=255)
                                       gCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACaaa                                                                                                                >  2:1452023/1‑150 (MQ=255)
                                        cAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGATAAAACCCACGCGAAACGCGAGGTGGTCAAACTCCGACTGTTTTAGCTCACaaaa                                                                                                               >  2:705346/1‑149 (MQ=255)
                                          gATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCa                                                                                                             >  2:515365/1‑150 (MQ=255)
                                            tGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCATGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgc                                                                                                           >  2:1566704/1‑150 (MQ=255)
                                            tGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgc                                                                                                           >  1:35078/1‑150 (MQ=255)
                                            tGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgc                                                                                                           >  1:655031/1‑150 (MQ=255)
                                              aaaCTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAgcgc                                                                                                         >  1:1118358/1‑150 (MQ=255)
                                                  tGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTaa                                                                                                     >  1:1361689/1‑150 (MQ=255)
                                                    cAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCAAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAAc                                                                                                   >  1:1248087/1‑150 (MQ=255)
                                                        tGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGGGTTAAAACCCACGCGAAACGCGTCGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACcgcg                                                                                               >  2:852730/1‑150 (MQ=255)
                                                        tGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACcgcg                                                                                               >  1:585970/1‑150 (MQ=255)
                                                          cAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTc                                                                                             <  2:403592/150‑1 (MQ=255)
                                                            ggCCGATCATGCCGACAACTTCGCCGTCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCgg                                                                                           <  2:520069/150‑1 (MQ=255)
                                                            ggCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGGGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCgg                                                                                           <  2:255321/150‑1 (MQ=255)
                                                                gATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTc                                                                                       >  2:853424/1‑150 (MQ=255)
                                                                   cATGCCGACAACTTCGCCATCAAGCAGCGACAGATGGTAGCGAATCTTTGGGTCCCGCAGATTGGAGTTAAAACCCACGATAAACGCGTGGTGGTCAAACTCCGACTGTTTTAGCTCACAAATCAGCGCGTAAAACGCGTCGGATTCGTa                                                                                    >  2:1116814/1‑150 (MQ=255)
                                                                              cTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCg                                                                         <  2:1361689/150‑1 (MQ=255)
                                                                              cTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCg                                                                         <  2:1004710/150‑1 (MQ=255)
                                                                               ttCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCgg                                                                        <  2:1118358/150‑1 (MQ=255)
                                                                                  gCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCg                                                                     <  2:925966/150‑1 (MQ=255)
                                                                                   ccATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGa                                                                    >  2:374067/1‑150 (MQ=255)
                                                                                         agcagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTca                                                              >  2:920767/1‑150 (MQ=255)
                                                                                           cagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTcaca                                                            >  1:843552/1‑150 (MQ=255)
                                                                                           cagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTcaca                                                            <  2:655031/150‑1 (MQ=255)
                                                                                           cagcGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTcaca                                                            <  2:1038995/150‑1 (MQ=255)
                                                                                                   gATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCa                                                    <  1:1331693/150‑1 (MQ=255)
                                                                                                      ggTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATgg                                                 >  1:1167954/1‑150 (MQ=255)
                                                                                                      ggTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATgg                                                 <  1:998513/150‑1 (MQ=255)
                                                                                                        tAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAAAGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgtg                                               >  2:798742/1‑150 (MQ=255)
                                                                                                         aGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGgtgt                                              <  1:515365/150‑1 (MQ=255)
                                                                                                          gcgcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCGCAAATCAGCGCGTAAAACGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGtt                                             >  2:231279/1‑150 (MQ=255)
                                                                                                           cgcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGttt                                            <  1:512257/150‑1 (MQ=255)
                                                                                                            gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtc                                           <  1:260989/150‑1 (MQ=255)
                                                                                                            gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtc                                           >  2:1003277/1‑150 (MQ=255)
                                                                                                            gcATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACCCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtc                                           <  1:463218/150‑1 (MQ=255)
                                                                                                              aTGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTtctc                                         >  1:595894/1‑150 (MQ=255)
                                                                                                               tGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTcc                                        <  2:838821/150‑1 (MQ=255)
                                                                                                               tGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTcc                                        <  1:1549505/150‑1 (MQ=255)
                                                                                                                gTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCt                                       <  1:920767/150‑1 (MQ=255)
                                                                                                                gTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCt                                       <  1:808442/150‑1 (MQ=255)
                                                                                                                  ttGGGGCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTGTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTc                                     <  1:1416500/150‑1 (MQ=255)
                                                                                                                   tGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCt                                    >  1:451771/1‑150 (MQ=255)
                                                                                                                   tGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCt                                    <  2:451771/150‑1 (MQ=255)
                                                                                                                    gggTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTg                                   <  2:1597376/150‑1 (MQ=255)
                                                                                                                      gTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGAt                                 >  1:57878/1‑150 (MQ=255)
                                                                                                                       tCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATg                                <  1:421453/150‑1 (MQ=255)
                                                                                                                        cgcgCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATgg                               >  1:1384051/1‑150 (MQ=255)
                                                                                                                        cgcgCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATgg                               <  2:1333942/150‑1 (MQ=255)
                                                                                                                         gcgcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGa                              <  1:1481299/150‑1 (MQ=255)
                                                                                                                         gcgcAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGa                              <  1:782993/150‑1 (MQ=255)
                                                                                                                              gATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGc                         >  1:582840/1‑150 (MQ=255)
                                                                                                                                   gCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTcag                    <  1:374067/150‑1 (MQ=255)
                                                                                                                                    cGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTcagc                   >  2:594262/1‑150 (MQ=255)
                                                                                                                                     gTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTcagca                  >  2:1561116/1‑150 (MQ=255)
                                                                                                                                       tgaccaccACGCGTTTCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCagcagc                >  1:1346609/7‑150 (MQ=255)
                                                                                                                                       tAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCagcagc                >  1:1367932/1‑150 (MQ=255)
                                                                                                                                           aCCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCgtgt            <  2:843552/150‑1 (MQ=255)
                                                                                                                                           aCCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCgtgt            >  1:775455/1‑150 (MQ=255)
                                                                                                                                            accACGCGTTTCGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTc           >  1:1626164/2‑150 (MQ=255)
                                                                                                                                             ccACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcg          >  1:323592/1‑150 (MQ=255)
                                                                                                                                                cgcgAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgacc       <  1:1404821/150‑1 (MQ=255)
                                                                                                                                                cgcgAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgacc       <  1:1609890/150‑1 (MQ=255)
                                                                                                                                                cgcgAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgacc       <  1:853424/150‑1 (MQ=255)
                                                                                                                                                 gcgAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgaccg      <  1:1254910/150‑1 (MQ=255)
                                                                                                                                                   gAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTcgaccgac    >  1:61248/1‑150 (MQ=255)
                                                                                                                                                     aaCGCGTGGTGGTCAAACTCCGACTGTTTTAGCTCACAAATCATAGCGTAAACCGCGTCTGTGTCGTACTGCGTGGCCGGGCGAAGATCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCATCGTCAGCAGCGGGTCGACCGACTg  >  2:804394/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCACCAACTCCTGAATTTCGCCGATCCAGTTGACATGATGCAGATGAAACTGCAAATGCAGGCCGATCATGCCGACAACTTCGCCATCAAGCAGCGCCAGATGGTAGCGCATGTTTGGGTCGCGCAGATTGGCGTTAAAACCCACGCGAAACGCGTGGTGGTCAAACTCCGCCTGTTTTAGCTCACAAATCAGCGCGTAAACCGCGTCGGTGTCGTACTGCGTGGCCGGGCGAAGCTCACAAGCAGGCATGGTGTTTCTCCTTCTGATGGATCAGCGTCAGCAGCGTGTCGACCGACTG  >  NC_000913/4315289‑4315587

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