Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,897,613 G→A 7.2% P36P (CCG→CCA cbrB → PRK09823 family inner membrane protein, creBC regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,897,6130GA 7.2% 286.7 / 9.7 110P36P (CCG→CCAcbrBPRK09823 family inner membrane protein, creBC regulon
Reads supporting (aligned to +/- strand):  ref base G (52/50);  new base A (5/3);  total (57/53)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGC  >  NC_000913/3897466‑3897760
                                                                                                                                                   |                                                                                                                                                   
gCTTAAATTGGCGCTCAAAGGTAAGTAAAGTGAGTTTGATATGTCTGTTTCACGTCGGGTAATAAATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGc                                                                                                                                                    <  1:1509095/149‑1 (MQ=255)
   tAAATTTGAGCTCTAAGTTAAGTAAATGTAGTTTTATATGTCTGTTTCACGTCGGGTAATACATCACGTACTTTATTTTGCTGTTTTATGACCTTTAATTGTTGTTCTGTTTCTTTTCCTCTACATATTCTTCGCAAAAGAACAGCTgg                                                                                                                                                 <  1:262852/149‑1 (MQ=255)
        tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGCGGTTCtt                                                                                                                                            <  2:1551047/149‑1 (MQ=255)
           ctcgcAAATGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGTTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTgg                                                                                                                                         <  1:494348/145‑1 (MQ=255)
           cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGCGGTTCttttg                                                                                                                                         <  1:188193/149‑3 (MQ=255)
              tCAAAGGGAAGTAAAGGGAGTTTGATATGTCTGTTTCACGGCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTg                                                                                                                                      <  1:998622/149‑1 (MQ=255)
               cAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGa                                                                                                                                     >  1:1314481/1‑149 (MQ=255)
               cAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGa                                                                                                                                     >  2:8925/1‑149 (MQ=255)
               cAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGGTCTTTGGGTGa                                                                                                                                     >  1:1020897/1‑149 (MQ=255)
                 aaGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGata                                                                                                                                   <  1:1355817/149‑1 (MQ=255)
                  aGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAAATGGTGTTCTGTATCTTGTGCTCTACATATTCTTCACAAAAGAACCACTGTTTCTTTGGTTGataa                                                                                                                                  >  1:389980/1‑149 (MQ=255)
                   ggTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGataat                                                                                                                                 <  2:1261284/149‑1 (MQ=255)
                    gtaagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGCCCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTATTCTCAAAAGACCCGCTGGTTCTTTGGGTGataata                                                                                                                                >  2:114797/1‑149 (MQ=255)
                    gtaagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGataata                                                                                                                                >  2:835706/1‑149 (MQ=255)
                    gtaagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGataata                                                                                                                                >  1:139096/1‑149 (MQ=255)
                     tatctaaAGGTAGTTAGATATGTCTGTCTCCCGTCTGGGAATACATCACGGTCATTATTTTGCAGTGTTAGGACCGTTAATTGGTGTTCTGTTTCTTTTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGTTGATAATAc                                                                                                                               <  1:304218/145‑1 (MQ=255)
                       agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACAt                                                                                                                             <  1:1449118/149‑1 (MQ=255)
                        gtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTAGTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                                            >  2:1352176/1‑149 (MQ=255)
                              ggAGTTTGATATTTCTGTTTCACGTCGTGTAATACATCACGGACTTTTTTTAGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCAACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAttt                                                                                                                      <  2:108183/149‑1 (MQ=255)
                              ggAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGCGGTTCttttg                                                                                                                                         <  2:1071526/130‑3 (MQ=255)
                              ggAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGCGGTTCttttg                                                                                                                                         >  1:1071526/1‑128 (MQ=255)
                                aggttGATATGTCTGTTGCACGGCGGGTTATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACAGCTGGTTCTTTGGGTGATAATACATCCAAttttt                                                                                                                    <  2:956071/146‑1 (MQ=255)
                                 gTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCAAAtttcta                                                                                                                   >  1:254566/1‑146 (MQ=255)
                                     gATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGGGTTGTGTTTCTTGTCCTCAACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCtt                                                                                                               >  2:1485020/1‑149 (MQ=255)
                                      atatGGCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAGGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTa                                                                                                              <  2:636642/149‑1 (MQ=255)
                                       tatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTTTTATGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTTATAATACATCCAATTTTTCTCTTAt                                                                                                             >  2:745669/1‑149 (MQ=255)
                                        atGTCTGTTTCACGTCGGTTAATACATGACGGACTTTATTTTGCAGTGTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAtt                                                                                                            <  2:305734/149‑1 (MQ=255)
                                        atGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAtt                                                                                                            >  1:719866/1‑149 (MQ=255)
                                           tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                         <  2:479673/149‑1 (MQ=255)
                                             tGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGa                                                                                                       <  1:515687/149‑1 (MQ=255)
                                               tttCACGTCGGGTAATACATAACGGACTTTATTATGCAGTTTTATGACCGATAATTGGTGATCTCTTTCTTGTCCTCTACATATTCTTCTCAAAAGAACCGCTGGTTCTTTGGGTCATAATACATCCAATTTTTCTATTATTGTCGATa                                                                                                     >  2:550258/1‑149 (MQ=255)
                                                   aCGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTa                                                                                                 >  1:1249739/1‑149 (MQ=255)
                                                     gTCGTGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTATTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACg                                                                                               <  1:628209/149‑1 (MQ=255)
                                                       cGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACggg                                                                                             <  2:762258/149‑1 (MQ=255)
                                                          gTAATACATCACGTAATTTATTTTGCAGTTTTAGTACCGTTTATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGTTAACTACGGGAGc                                                                                          <  1:878559/149‑1 (MQ=255)
                                                          gTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGc                                                                                          <  1:1109420/149‑1 (MQ=255)
                                                          gTAATAAATCACTGACTTTAAATTGCATTTTTAGTACAGTTAATTGGTTTTCTGTTTCTTGTCCTCTAAATATTCTTCGCAAAATAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGc                                                                                          <  1:1112085/149‑1 (MQ=255)
                                                                 aTCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCt                                                                                   >  1:506656/1‑149 (MQ=255)
                                                                 aTCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCATCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCt                                                                                   >  1:263680/1‑149 (MQ=255)
                                                                 aTCACGGACTTTAATTTGCAGTTTTAGGACCGTTACTAGGTGTTCTGTTTCTTGTCCGCTACATATTCTACGCCACCGACCCGCTGGATCTTTGGGTGATAATACATCCAATTGTTCTCTTATAGTCGAAAACTACGGGAGCTATTCCt                                                                                   <  2:326100/149‑1 (MQ=255)
                                                                  tCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTg                                                                                  >  2:837408/1‑149 (MQ=255)
                                                                    aCGGACTTTATTTTGCAGTTTTAGGACCGTTAATTTGTGTTCTGTTTATTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGTTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCg                                                                                <  2:930648/149‑1 (MQ=255)
                                                                    aCGGACCTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCg                                                                                >  2:544273/1‑149 (MQ=255)
                                                                    aCGGACCTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTGTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGGTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCg                                                                                >  2:245721/1‑149 (MQ=255)
                                                                       gACTTTATTTTGCAGTTTTAGGACCGTTAATTTGTGTTCTGTGTCTTGTCCTCTACATATTCTTCTCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttg                                                                             <  1:1534165/149‑1 (MQ=255)
                                                                         cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgtt                                                                           >  2:1264987/1‑149 (MQ=255)
                                                                              ttttGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCg                                                                      <  2:1371486/149‑1 (MQ=255)
                                                                              ttttGCAGTTTTAGGACCGGTAATTGGTGTGCTGTTTCTTGTACTATACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCg                                                                      <  2:417417/149‑1 (MQ=255)
                                                                                 tGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACAACGGGAGCTATTCCTGCGTTGTTAACCGgtg                                                                   >  1:185683/1‑149 (MQ=255)
                                                                                      ttttAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTATTATTGTCGATAACTACGGGAGCTATTCCTGCTTTGTTAAACGGTGTAAtt                                                              >  1:464410/1‑148 (MQ=255)
                                                                                           ggACCGTTAATTGGTGGTCTGTTTCTAGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAATTTTCTCTTATTGACTATAACTACGGTAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGc                                                         <  1:1106804/149‑1 (MQ=255)
                                                                                            gACCGTTAATTGGTGTTCTGTATCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTCATAATAAATCCAAATTATCTCTTATTGTCGATAACTACGGGATCTATTCCTGCGTTCTTAACCGGTTTAATGGTTGca                                                        >  2:791657/1‑148 (MQ=255)
                                                                                             aCCGTTATTTGGTGTTCTTTTTCTTGTCCTCTACATCTTCTTCGCAAAAGAACCGCTGGTTCTTTGTGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCt                                                       <  2:1220007/149‑1 (MQ=255)
                                                                                             aCCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCt                                                       <  2:1249739/149‑1 (MQ=255)
                                                                                                gTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTc                                                    >  2:364905/1‑149 (MQ=255)
                                                                                                 ttAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCt                                                   >  1:156184/1‑149 (MQ=255)
                                                                                                   aaTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTTCCTGTCTGc                                                 >  1:1540099/1‑149 (MQ=255)
                                                                                                    aTTGGTGTTCTGTTTCTTGTCATCTACATATTCTTCGAAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTATCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcc                                                <  2:292297/149‑1 (MQ=255)
                                                                                                    aTTGGTGTTCTGTTTCTTGTACTCTACATATTCTTCGCAAATGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcc                                                <  1:460995/149‑1 (MQ=255)
                                                                                                     ttGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCg                                               <  2:778459/149‑1 (MQ=255)
                                                                                                       ggTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGa                                             >  1:985748/1‑149 (MQ=255)
                                                                                                        gtgtTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTTACCGGTGTAATGGTTGCCTGTCTGCCGGaa                                            >  1:216525/1‑149 (MQ=255)
                                                                                                        gtgtTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGGGTAATGGTTGCCTGTCTGCCGGaa                                            >  2:504083/1‑149 (MQ=255)
                                                                                                         tgtTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaa                                           >  2:140208/1‑149 (MQ=255)
                                                                                                           ttCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAg                                         <  1:978730/149‑1 (MQ=255)
                                                                                                           ttCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGGGTAATGGTTGCCTGTCTGCCGGAAAAg                                         <  1:837459/149‑1 (MQ=255)
                                                                                                             cTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGAt                                       <  2:1314481/149‑1 (MQ=255)
                                                                                                               gTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATTACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCg                                     >  2:1266399/1‑149 (MQ=255)
                                                                                                               gTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCg                                     <  1:1204627/149‑1 (MQ=255)
                                                                                                                tttCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                    <  1:871367/149‑1 (MQ=255)
                                                                                                                tttCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCgg                                    <  1:243207/149‑1 (MQ=255)
                                                                                                                  tCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGt                                  <  2:409283/149‑1 (MQ=255)
                                                                                                                    ttGTCCTCTACATATTCTTCGCAAAAGACCCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTca                                >  1:244080/1‑149 (MQ=255)
                                                                                                                    ttGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGGGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTca                                >  1:102275/1‑149 (MQ=255)
                                                                                                                     tGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTCCGTTGTTAACCGGTGTAATGGTTGCCTGTCATCCGGAAAAGATCGGGTcac                               >  2:316865/1‑149 (MQ=255)
                                                                                                                         ctctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGGCGGAAAAGATCGGGTCACAGaa                           >  2:908139/1‑149 (MQ=255)
                                                                                                                         ctctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaa                           >  2:819968/1‑149 (MQ=255)
                                                                                                                         ctctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTACTGCGTTGTTAACCTGTGTAATGATTGCATGTCTGCCTGAAAAGATCGGGTCACAGaa                           >  1:1008186/1‑149 (MQ=255)
                                                                                                                         ctctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTAATCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaa                           <  1:835706/149‑1 (MQ=255)
                                                                                                                         ctctACATATTCTTCGCAAAAGAACCGCTGGTAATTTGGATGATAATACATCCAATTTTTCTCTTATTGTCTATAACTACGGGAGCTATTCATGCGTATTTAACAGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaa                           >  1:392519/1‑149 (MQ=255)
                                                                                                                           ctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATTATACATCCAATTTTTCTCTTATTGTCGATAACTACTGGAGCTATTACTGCGTTGTTAACCGCTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAAc                         >  1:419249/1‑149 (MQ=255)
                                                                                                                             aCATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATTACTACGGGAGCTATTCCTGAGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGt                       >  1:1538116/1‑149 (MQ=255)
                                                                                                                               atatTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATACCTACGGGAGCTATTCCTCCGTTTTTAACCGGTGTAATGGTTGCATGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa                     >  2:852919/1‑149 (MQ=255)
                                                                                                                               atatTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTa                     <  1:1282409/149‑1 (MQ=255)
                                                                                                                                 atTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTATTATTGTAGATAACTACGGGAGCTATTCCAGCGATGTTAACCGGTGTAATGTTTGCCTGTCTGCAGGAAAATATCTGGTCACAGAAACGTTAta                   >  2:858793/1‑148 (MQ=255)
                                                                                                                                   tcttcGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGt                 >  2:507820/1‑149 (MQ=255)
                                                                                                                                    cttcGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGtt                <  1:121694/149‑1 (MQ=255)
                                                                                                                                       cGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGTGAGCTATTCCTGCGTTGTTAACCAGTGTAATGGTTGCCTGTCTGCCGCAAAAGATCTGGTCACAGAAACGATATCGTTGta             >  2:1382692/1‑148 (MQ=255)
                                                                                                                                        gCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTACTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTct            >  2:1460081/1‑149 (MQ=255)
                                                                                                                                         ccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg           >  1:770405/3‑149 (MQ=255)
                                                                                                                                          caaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCTCAGAAACGTTATCGTTGTctgg          >  2:824497/2‑149 (MQ=255)
                                                                                                                                            aaGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctggct        <  1:775686/149‑1 (MQ=255)
                                                                                                                                             aGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctggctg       <  2:407150/149‑1 (MQ=255)
                                                                                                                                              gAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctggctgg      >  2:1520245/1‑149 (MQ=255)
                                                                                                                                               aaCCGCTGGTTCTTTGGGTGATAATAGATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCtggt     <  2:185683/149‑1 (MQ=255)
                                                                                                                                               aaCCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCtggt     <  2:156184/149‑1 (MQ=255)
                                                                                                                                                 ccGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCtggtgg   >  2:325764/1‑149 (MQ=255)
                                                                                                                                                 ccGCGGGTTCTTTGGGTGATAATTCTTCGTATTTTTGTCTTATTGTCGATATCTTCGGGAGCTATTCCTGCGTTGTTATCCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCtggtgg   <  2:1082541/149‑1 (MQ=255)
                                                                                                                                                  cGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCTGTGTAATGGTTGCCTGTATGACGGAAAAGATCGGGTCACATAAACGTTATCGTTGTCTGGCTTGTGGc  >  1:1377776/1‑149 (MQ=255)
                                                                                                                                                  cGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGc  <  1:8925/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGC  >  NC_000913/3897466‑3897760

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