Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,897,615 T→C 6.4% L37P (CTG→CCG) ‡ 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,6150TC 6.4% 304.3 / 9.7 109L37P (CTG→CCG) ‡cbrBPRK09823 family inner membrane protein, creBC regulon
Reads supporting (aligned to +/- strand):  ref base T (53/48);  new base C (4/3);  total (57/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 109 total aligned reads displayed.

TAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGC  >  NC_000913/3897469‑3897760
                                                                                                                                                  |                                                                                                                                                 
tAAATTTGAGCTCTAAGTTAAGTAAATGTAGTTTTATATGTCTGTTTCACGTCGGGTAATACATCACGTACTTTATTTTGCTGTTTTATGACCTTTAATTGTTGTTCTGTTTCTTTTCCTCTACATATTCTTCGCAAAAGAACAGCTgg                                                                                                                                                 <  1:262852/149‑1 (MQ=255)
     tGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCtt                                                                                                                                            >  1:305734/1‑149 (MQ=255)
        ctcgcAAATGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGTTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTgg                                                                                                                                         <  1:494348/145‑1 (MQ=255)
        cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGCGGTTCttttg                                                                                                                                         <  1:188193/149‑3 (MQ=255)
           tCAAAGGGAAGTAAAGGGAGTTTGATATGTCTGTTTCACGGCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTg                                                                                                                                      <  1:998622/149‑1 (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)
                 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)
                       aaaGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCa                                                                                                                          >  2:243207/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)
                                                       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)
                                                                   ggACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgtt                                                                              <  1:1544044/149‑1 (MQ=255)
                                                                    gACTTTATTTTGCAGTTTTAGGACCGTTAATTTGTGTTCTGTGTCTTGTCCTCTACATATTCTTCTCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttg                                                                             <  1:1534165/149‑1 (MQ=255)
                                                                      cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttgtt                                                                           >  2:1264987/1‑149 (MQ=255)
                                                                      cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTAATGTCGATAACTACGGGAGCTATTCCTGCgttgtt                                                                           >  1:767205/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)
                                                                                        ggACCGTTAATTGGTGGTCTGTTTCTAGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAAATTTTCTCTTATTGACTATAACTACGGTAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGc                                                         <  1:1106804/149‑1 (MQ=255)
                                                                                          aCCGTTATTTGGTGTTCTTTTTCTTGTCCTCTACATCTTCTTCGCAAAAGAACCGCTGGTTCTTTGTGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCt                                                       <  2:1220007/149‑1 (MQ=255)
                                                                                          aCCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCt                                                       <  2:1249739/149‑1 (MQ=255)
                                                                                            cGTTCATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGt                                                     <  2:132949/149‑1 (MQ=255)
                                                                                             gTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTc                                                    >  2:440924/1‑149 (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)
                                                                                                     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)
                                                                                                                 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)
                                                                                                                                     gCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTct            <  2:797813/149‑1 (MQ=255)
                                                                                                                                     gCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTACTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTct            >  2:1460081/1‑149 (MQ=255)
                                                                                                                                     gCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTGATAGTCTATAACTACGGGAGCTATTCCTGCGTTGTTAAACGGTGTAAGGGTTGACTGTCTGACGTAAAAGATCGGGTCACAGAAACGTGAGCGTTGTct            >  1:1256008/1‑149 (MQ=255)
                                                                                                                                      ccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg           >  1:770405/3‑149 (MQ=255)
                                                                                                                                      ccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg           >  2:1052500/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)
                                                                                                                                                  |                                                                                                                                                 
TAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGC  >  NC_000913/3897469‑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: