Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,898,322 G→A 6.2% intergenic (+26/‑105) sdaA → / → pdeD L‑serine deaminase I/putative c‑di‑GMP phosphodiesterase PdeD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,898,3220GA 6.2% 300.4 / 6.3 97intergenic (+26/‑105)sdaA/pdeDL‑serine deaminase I/putative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base G (54/37);  new base A (2/4);  total (56/41)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

CAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTCGTAATCGA  >  NC_000913/1898173‑1898471
                                                                                                                                                     |                                                                                                                                                     
cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACg                                                                                                                                                       <  2:366895/150‑1 (MQ=255)
cAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACg                                                                                                                                                       <  1:215709/150‑1 (MQ=255)
    gCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATg                                                                                                                                                   >  1:1395450/1‑150 (MQ=255)
     cACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATgg                                                                                                                                                  >  1:1266287/1‑150 (MQ=255)
        cgcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCg                                                                                                                                               >  2:38763‑M2/7‑150 (MQ=255)
         gcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGa                                                                                                                                              >  1:987700‑M2/6‑150 (MQ=255)
         gcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGa                                                                                                                                              >  1:1104227‑M2/6‑150 (MQ=255)
         gcgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACAGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGa                                                                                                                                              >  2:595998‑M2/6‑150 (MQ=255)
          cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCgag                                                                                                                                             <  1:797433‑M2/146‑2 (MQ=255)
          cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCgag                                                                                                                                             <  2:535679‑M2/146‑2 (MQ=255)
          cgcgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGGGCgag                                                                                                                                             <  1:258629‑M2/146‑2 (MQ=255)
            cgTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAAtt                                                                                                                                           <  2:879686‑M2/148‑1 (MQ=255)
             gTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCCGTTGCAGATGggcgagtaa                                                                                                                                          <  1:767563‑M2/149‑5 (MQ=255)
              tctcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTta                                                                                                                                         <  1:1328229/150‑1 (MQ=255)
                tcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCAAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTtata                                                                                                                                       <  2:69869/150‑1 (MQ=255)
                gcGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTtata                                                                                                                                       <  1:289649/149‑1 (MQ=255)
                 cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATAc                                                                                                                                      >  1:1412417/1‑150 (MQ=255)
                 cGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCCCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATAc                                                                                                                                      >  2:1158384/1‑150 (MQ=255)
                   cTGGATAAGGTCATCGAACCGATGGCCGAAACCGGTAAGGCCATGAACGCCAAATACGGCGAAACCTCCCGCGGTGGTCTGGCACTCGCTGTCCAGTGTGAGGTATAGTTCTTACTCGCCCATCTGCACCGGAGGGGCGAATGTATGccc                                                                                                                                    >  2:1572230/1‑150 (MQ=255)
                    tGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCg                                                                                                                                   <  2:913252/150‑1 (MQ=255)
                     ggATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGc                                                                                                                                  >  2:389817/1‑150 (MQ=255)
                          aGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtct                                                                                                                             <  1:1316178/150‑1 (MQ=255)
                           ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtctc                                                                                                                            >  2:1040736/1‑150 (MQ=255)
                           ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtctc                                                                                                                            >  1:941964/1‑150 (MQ=255)
                           ggTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTtctc                                                                                                                            <  1:771071/150‑1 (MQ=255)
                            gTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCACATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCg                                                                                                                           >  1:1266160/1‑150 (MQ=255)
                               aTCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTct                                                                                                                        <  2:927307/150‑1 (MQ=255)
                                tCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctg                                                                                                                       >  2:1121910/1‑150 (MQ=255)
                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTta                                                                                                                                         <  1:1603232/130‑1 (MQ=255)
                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTta                                                                                                                                         >  2:1603232/1‑130 (MQ=255)
                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgct                                                                                                                     >  1:728381/1‑150 (MQ=255)
                                  gAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgct                                                                                                                     >  1:1580551/1‑150 (MQ=255)
                                   aaaCGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTctgctg                                                                                                                    <  2:1104227/150‑1 (MQ=255)
                                          gTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCAGGCGGTGGTGTGGCAATCAAAGTCCAGTGTGGCTAAAACTTCTTACTCGCCCAGCTGCAACGGATGGGCGAATTTATACCCGCTTTCGCGTCTGCTGTAATAtt                                                                                                             >  2:1506959/1‑150 (MQ=255)
                                           tACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATAATTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATAtta                                                                                                            >  2:1411256/1‑149 (MQ=255)
                                             cGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTccc                                                                                                          >  2:1544522/1‑150 (MQ=255)
                                              gAAACCTGTAAGGAGATGATCGCAAATTACCGCTTAACCTCACGGGGTGGTTTGGCAATCAAAGTCCATTGTTACTAATACTTGTTACGCGCCCTTGTGCAACGGATTTGCGAATTTATACGCGCTTTCTCGTCTGCTGTAATATTgccc                                                                                                         <  2:978437/150‑1 (MQ=255)
                                              gAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTcccc                                                                                                         <  2:737701/150‑1 (MQ=255)
                                                aaCCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCAc                                                                                                       <  1:1495297/150‑1 (MQ=255)
                                                   cGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTac                                                                                                    <  1:1473168/150‑1 (MQ=255)
                                                    ggTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTaca                                                                                                   >  2:1259388/1‑150 (MQ=255)
                                                       aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTac                                                                                                    <  2:1566426/146‑1 (MQ=255)
                                                       aaGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTac                                                                                                    >  1:1566426/1‑146 (MQ=255)
                                                           aCATGAACGACAAATACCGCGAAACCTCACGCGGTGGTCTGGAAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATAACCGCTTTCTCGTCTGCTGTAATATTACCAACTACACTTCCaa                                                                                            >  2:799050/1‑149 (MQ=255)
                                                             aTGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTg                                                                                          >  2:1021776/1‑150 (MQ=255)
                                                              tGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGt                                                                                         >  1:914641/1‑150 (MQ=255)
                                                                aaCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTg                                                                                       >  2:621946/1‑150 (MQ=255)
                                                                 aCGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGc                                                                                      >  1:376603/1‑150 (MQ=255)
                                                                     cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCCTCTGCAACGCATGGGCGAATTTATACCCGCTGTCTCGTCGGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCa                                                                                  <  2:983091/150‑1 (MQ=255)
                                                                     cAAATACCGCGAAACCTCACGCGGTGGTCTGGCAACCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCa                                                                                  >  2:356419/1‑150 (MQ=255)
                                                                       aaTACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAgg                                                                                <  2:1312033/150‑1 (MQ=255)
                                                                             gcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTg                                                                          >  2:293556/1‑150 (MQ=255)
                                                                             gcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTg                                                                          <  2:1412417/150‑1 (MQ=255)
                                                                             gcgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTg                                                                          <  1:1081247/150‑1 (MQ=255)
                                                                              cgAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGt                                                                         <  2:961288/150‑1 (MQ=255)
                                                                                  aCCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGcc                                                                     <  1:293556/150‑1 (MQ=255)
                                                                                    cTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCAt                                                                   >  1:621423/1‑150 (MQ=255)
                                                                                    cTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCAt                                                                   >  2:15632/1‑150 (MQ=255)
                                                                                          cgGTGGTCTCCCAATCAAAGTAAACTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCCGGCGTTTGTCGCCATACGcgt                                                             >  2:1426027/1‑148 (MQ=255)
                                                                                             tggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACa                                                          <  1:1040736/150‑1 (MQ=255)
                                                                                             tggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACa                                                          <  2:376603/150‑1 (MQ=255)
                                                                                              ggtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAg                                                         >  1:440142/1‑150 (MQ=255)
                                                                                               gtCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAgg                                                        <  1:1169810/150‑1 (MQ=255)
                                                                                                tCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAggg                                                       <  1:738146/150‑1 (MQ=255)
                                                                                                tCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAggg                                                       >  2:147487/1‑150 (MQ=255)
                                                                                                  tGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTg                                                     >  2:237511/1‑150 (MQ=255)
                                                                                                  tGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTg                                                     <  1:1474898/150‑1 (MQ=255)
                                                                                                   ggCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTgg                                                    >  2:538645/1‑150 (MQ=255)
                                                                                                      aaTCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGccc                                                 >  2:267152/1‑150 (MQ=255)
                                                                                                        tCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTAAAGGGTGGCCCGc                                               >  2:639218/1‑150 (MQ=255)
                                                                                                          aaaGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCAt                                             <  1:38763/150‑1 (MQ=255)
                                                                                                          aaaGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCAt                                             >  1:347608/1‑150 (MQ=255)
                                                                                                          aaaGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTACCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCAt                                             >  2:65585/1‑150 (MQ=255)
                                                                                                          aaaGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCATTCACGTCTGCTGTAATATTCCACACTACACTTACACTGTTTCCTCAGGCGTTTGTCGCCAAACGCTTACAGGATGCACCGCAt                                             >  2:5707/1‑150 (MQ=255)
                                                                                                           aaGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATg                                            <  2:368647/150‑1 (MQ=255)
                                                                                                             gTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCa                                          <  1:268324/150‑1 (MQ=255)
                                                                                                                cAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGGTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCaaaa                                       >  2:1432828/1‑150 (MQ=255)
                                                                                                                 tgtgtGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCTCTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCaaaat                                      <  2:1266160/149‑2 (MQ=255)
                                                                                                                 aGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCaaaaa                                      >  1:1283545/1‑150 (MQ=255)
                                                                                                                  gtgtGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAg                                     >  1:57/1‑150 (MQ=255)
                                                                                                                  gtgtGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAg                                     >  1:517302/1‑150 (MQ=255)
                                                                                                                   tgtgACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGc                                    <  1:65585/150‑1 (MQ=255)
                                                                                                                   tgtgACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGc                                    <  1:561514/150‑1 (MQ=255)
                                                                                                                     tgACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGcac                                  >  1:1213074/1‑150 (MQ=255)
                                                                                                                      gACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGcaca                                 >  1:782115/1‑150 (MQ=255)
                                                                                                                      gACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGcaca                                 >  2:405290/1‑150 (MQ=255)
                                                                                                                      gACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCAATGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGcaca                                 >  1:453083/1‑150 (MQ=255)
                                                                                                                      gACTAATACATCTAACTCACACATCTCCAACGGATAATCGAATTTAAACAAGATTTCAAGTCTGCTGTAATATTCCCCACTACAATTCCACTCTTGCGTCATTCGTTTGTCGCCATACGCTTACAGTCTGGCCCGCATGCAAAAAGcaca                                 >  2:1488151/1‑150 (MQ=255)
                                                                                                                        cTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAAc                               >  1:941775/1‑150 (MQ=255)
                                                                                                                                    attCGCCCATCCGTTGCAGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTaaaa                   >  1:853768/3‑150 (MQ=255)
                                                                                                                                    attCGCCCATCCGTTGCAGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTaaaa                   >  1:669062/3‑150 (MQ=255)
                                                                                                                                     cTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAc                  <  1:1121910/150‑1 (MQ=255)
                                                                                                                                     cTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAAc                  <  1:621946/150‑1 (MQ=255)
                                                                                                                                        gCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTa               >  1:1240225/1‑150 (MQ=255)
                                                                                                                                        gCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTa               >  2:1586963/1‑150 (MQ=255)
                                                                                                                                         cccATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTAt              <  2:897826/150‑1 (MQ=255)
                                                                                                                                          ccATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCAGAACGGATCATTAAAACCTAtc             <  2:1580551/150‑1 (MQ=255)
                                                                                                                                               tGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATcgtcgt        >  2:334776/1‑150 (MQ=255)
                                                                                                                                                  aaCGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTCGTAAt     <  1:267152/150‑1 (MQ=255)
                                                                                                                                                     ggATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGACGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTCATAATCGa  <  2:728381/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGTGCACCGCGCGTCTCGCTGGATAAGGTCATCGAAACGATGTACGAAACCGGTAAGGACATGAACGCCAAATACCGCGAAACCTCACGCGGTGGTCTGGCAATCAAAGTCCAGTGTGACTAATACTTCTTACTCGCCCATCTGCAACGGATGGGCGAATTTATACCCGCTTTCTCGTCTGCTGTAATATTCCCCACTACACTTCCACTGTTGCGTCAGGCGTTTGTCGCCATACGCTTACAGGGTGGCCCGCATGCAAAAAGCACAACGGATCATTAAAACCTATCGTCGTAATCGA  >  NC_000913/1898173‑1898471

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