Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,897,614 C→G 7.7% L37V (CTG→GTG) ‡ 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,6140CG 7.7% 480.6 / 26.0 156L37V (CTG→GTG) ‡cbrBPRK09823 family inner membrane protein, creBC regulon
Reads supporting (aligned to +/- strand):  ref base C (58/84);  new base G (7/5);  total (67/89)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.20e-02

Only 100 of 157 total aligned reads displayed.

GCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATA  >  NC_000913/3897466‑3897763
                                                                                                                                                    |                                                                                                                                                     
gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCt                                                                                                                                                      <  2:1434845/150‑1 (MQ=255)
gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCt                                                                                                                                                      <  1:471365/150‑1 (MQ=255)
gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCt                                                                                                                                                      <  1:1190655/150‑1 (MQ=255)
gCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCt                                                                                                                                                      >  2:827873/1‑150 (MQ=255)
 cTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTg                                                                                                                                                     >  2:32169/1‑150 (MQ=255)
  ttAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTgg                                                                                                                                                    <  2:1495093/150‑1 (MQ=255)
     aaTTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTc                                                                                                                                                 <  2:1593002/150‑1 (MQ=255)
         ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTg                                                                                                                                             <  2:1034930/150‑1 (MQ=255)
         ggCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGGGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTg                                                                                                                                             <  1:469654/150‑1 (MQ=255)
          gcgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGCGGTTCttttg                                                                                                                                            <  2:489900/150‑3 (MQ=255)
           cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGCGGTTCttttgc                                                                                                                                           <  2:1633279/150‑4 (MQ=255)
           cgcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTGCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCAGCGGTTCttttgc                                                                                                                                           <  1:1363940/150‑4 (MQ=255)
            tggcAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGt                                                                                                                                          <  1:1270914/147‑1 (MQ=255)
            gcTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGt                                                                                                                                          <  2:1574512/150‑1 (MQ=255)
              tCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGTACCGTTAATAGGTGTTCTGTTTCTTGTCCTCTACATATTATTCGCAAAAGAACCGCTGGTAATATGGGTGa                                                                                                                                        >  2:389092/1‑150 (MQ=255)
                aaaGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTTATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGata                                                                                                                                      <  1:1378535/150‑1 (MQ=255)
                   ggTAAGTAAAGGGAGTTTGATATGTCTGTTTCAGGTCGGGTAATACATCACGGACTTTATTTTGCAGTTGTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGataata                                                                                                                                   <  1:828171/150‑1 (MQ=255)
                   ggTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGataata                                                                                                                                   >  2:303981/1‑150 (MQ=255)
                   ggTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGataata                                                                                                                                   <  1:1157604/150‑1 (MQ=255)
                      aagtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACAt                                                                                                                                >  2:596317/1‑150 (MQ=255)
                       tgtatAGGGAGTTTGATATGTCTGTTTGACTTCTGGTAATACATCAATTACTTTATTTTGCAGTTTTATGACCGTTATTTGTTGTTATGTTTCTTGTCCTATACATATTCTTCGCATAAGATGCGCTGGTTATTTGGGTGATAATACATc                                                                                                                               <  2:1285622/149‑1 (MQ=255)
                       agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                                               <  1:440643/150‑1 (MQ=255)
                       agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                                               <  2:446696/150‑1 (MQ=255)
                       agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                                               <  1:1249232/150‑1 (MQ=255)
                       agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                                               <  1:1415056/150‑1 (MQ=255)
                       agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                                               <  2:711575/150‑1 (MQ=255)
                       agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                                               >  2:629794/1‑150 (MQ=255)
                       agtaaAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATc                                                                                                                               <  1:1593014/150‑1 (MQ=255)
                          aaaGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGGAGTTTAAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCaaa                                                                                                                            >  1:175420/1‑150 (MQ=255)
                               gAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATAAATCAAAttttt                                                                                                                       >  2:1604272/1‑150 (MQ=255)
                                aGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTtc                                                                                                                      <  1:516564/150‑1 (MQ=255)
                                      atatGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTACTCTACATATTCTTCGCAAAAGAACCGCTGGATCTTTGGGTGATAATACATCCAATTTTTCTCTTAt                                                                                                                >  1:1498436/1‑150 (MQ=255)
                                      atatCTCTGTTTCACGTCGGGTAATACATAACCGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTACTCTACATATTCTTCGCAAAAGAACCGCTTGTTCTTTGGTTTATAATACATCCAATTTTTCACTTAt                                                                                                                >  2:819309/1‑150 (MQ=255)
                                          gTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTc                                                                                                            <  1:1017987/150‑1 (MQ=255)
                                           tCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCg                                                                                                           >  1:635470/1‑150 (MQ=255)
                                             tGTTTCACGTCGGGTAATACTTCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGAt                                                                                                         <  1:985671/150‑1 (MQ=255)
                                              gTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATa                                                                                                        <  1:313535/150‑1 (MQ=255)
                                              gTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATa                                                                                                        <  2:1418086/150‑1 (MQ=255)
                                              gTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACAGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATa                                                                                                        >  1:171177/1‑150 (MQ=255)
                                               tttCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATaa                                                                                                       <  2:224319/150‑1 (MQ=255)
                                                 tCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACa                                                                                                                                 <  1:1705550/122‑1 (MQ=255)
                                                   aCGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTAc                                                                                                   <  2:416630/150‑1 (MQ=255)
                                                   aCGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTAc                                                                                                   <  2:372726/150‑1 (MQ=255)
                                                     gTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACgg                                                                                                 >  1:1293664/1‑150 (MQ=255)
                                                        gggTAATACCTCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                                              <  2:1374223/150‑1 (MQ=255)
                                                        gggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCTCAAAAGAACCGCTGGTTCTTTGTGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                                              <  2:30379/150‑1 (MQ=255)
                                                        gggTAATACATCACGGAATTTATTTTGCAGTTTTAGTACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAg                                                                                              <  1:716509/150‑1 (MQ=255)
                                                         ggTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGc                                                                                             <  2:1310562/150‑1 (MQ=255)
                                                               aCATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTcc                                                                                       >  2:862050/1‑150 (MQ=255)
                                                               aCATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTcc                                                                                       <  2:1424849/150‑1 (MQ=255)
                                                                cATCACTGACTTTATTTTGCATTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTATACATATTCTTCGCAAAAGAACCGCTGGTTCTTTTGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCt                                                                                      <  1:426728/150‑1 (MQ=255)
                                                                      ggACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCgttg                                                                                <  1:1021855/150‑1 (MQ=255)
                                                                         cTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTa                                                                             >  2:1089326/1‑150 (MQ=255)
                                                                           ttATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAc                                                                           >  2:551899/1‑150 (MQ=255)
                                                                            tATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAcc                                                                          <  1:745977/150‑1 (MQ=255)
                                                                            tATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAAcc                                                                          <  1:629794/150‑1 (MQ=255)
                                                                            tATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCg                                                                         >  2:753453/1‑150 (MQ=255)
                                                                             aTTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCgg                                                                        <  2:404227/150‑1 (MQ=255)
                                                                                ttGCAGTTTTATGAACGTTAATTGGTTTTCTGTTTATTGTCCTCTACATATTCTTCTCAAAAGAAACGCTGTTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgtg                                                                      <  1:325485/150‑1 (MQ=255)
                                                                                  gCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTa                                                                    >  1:385059/1‑150 (MQ=255)
                                                                                   cAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTaa                                                                   >  2:1709546/1‑150 (MQ=255)
                                                                                         tAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAAACGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTg                                                             <  2:883718/150‑1 (MQ=255)
                                                                                           ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgtg                                                                      >  1:1378931/1‑139 (MQ=255)
                                                                                           ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGgtg                                                                      <  2:1378931/139‑1 (MQ=255)
                                                                                           ggACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGcc                                                           <  1:698617/150‑1 (MQ=255)
                                                                                             aCCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTg                                                         >  2:1013740/1‑150 (MQ=255)
                                                                                                   aaTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcc                                                   <  1:32169/150‑1 (MQ=255)
                                                                                                   aaTTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGGTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGcc                                                   <  1:382906/150‑1 (MQ=255)
                                                                                                   aaTTGGTGTTATGTTTCTTTTACTCTACATATTCTTAGTAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGAAAACTACGGGAGCTTTTCCTGCGTTGTTAAGCGGTGTAATGGTTGCATGTCTGcc                                                   <  1:819309/150‑1 (MQ=255)
                                                                                                         tgtTATGTTTATTGTCCTCTACATATTCTTCGAAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaaa                                             <  2:1056114/150‑1 (MQ=255)
                                                                                                          gtTCTGTTTCTTGTCCTCTACATATTCTTCGCAAATGAACCGCTGGTTCTTTGGGTGATAATACAACCATTTTTTCTCTTATTGTCGATAACTACGGTAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAg                                            <  1:843999/150‑1 (MQ=255)
                                                                                                            tCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGAt                                          >  1:237443/1‑150 (MQ=255)
                                                                                                                  tCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTc                                    <  2:635470/150‑1 (MQ=255)
                                                                                                                    ttGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTcac                                  <  1:1192500/150‑1 (MQ=255)
                                                                                                                     tGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTGATGGTTGCCTGTCTGCCGGAAAAGATCGGGTcaca                                 >  2:109712/1‑150 (MQ=255)
                                                                                                                        ccTCTACATATTCTTCGCAAAAGAACCACTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaa                              >  1:246877/1‑150 (MQ=255)
                                                                                                                         ctctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGaaa                             <  2:668981/150‑1 (MQ=255)
                                                                                                                           ctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACg                           >  2:1240389/1‑150 (MQ=255)
                                                                                                                           ctACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACg                           >  2:162035/1‑150 (MQ=255)
                                                                                                                            tACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGt                          >  2:974884/1‑150 (MQ=255)
                                                                                                                                tatTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaaa                                             >  1:807771/1‑127 (MQ=255)
                                                                                                                                tatTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGaaaa                                             <  2:807771/127‑1 (MQ=255)
                                                                                                                                    cttcGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTg                  <  1:1240389/150‑1 (MQ=255)
                                                                                                                                    cttcGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTg                  >  1:1458840/1‑150 (MQ=255)
                                                                                                                                    attcGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTg                  <  1:162035/149‑1 (MQ=255)
                                                                                                                                      cacccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTc                >  1:1582811/6‑150 (MQ=255)
                                                                                                                                        gCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTc                <  2:355050/148‑1 (MQ=255)
                                                                                                                                        gCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGGTGCCTGGCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGGTGTctg              <  1:1359543/150‑1 (MQ=255)
                                                                                                                                        cccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCTATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTg                  >  2:919218/4‑146 (MQ=255)
                                                                                                                                        cccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCTATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTg                  <  1:919218/143‑1 (MQ=255)
                                                                                                                                        cccaaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctg              >  1:1683161/4‑150 (MQ=255)
                                                                                                                                         cAAAAGAAGCGCTGTTTATTTGGGTTATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTTCGTTGTTAACCGGTTTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctgg             <  2:175420/150‑1 (MQ=255)
                                                                                                                                         cAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctgg             <  2:1520416/150‑1 (MQ=255)
                                                                                                                                          caaaGAACCAGCGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctggc            >  1:696967/2‑150 (MQ=255)
                                                                                                                                          caaaGAACCAGCGGTTCTTTGGGTGATAATACATACAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctggc            >  1:703180/2‑150 (MQ=255)
                                                                                                                                          aaaaGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCTGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTctggc            >  2:1315208/1‑150 (MQ=255)
                                                                                                                                             aGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCGTAGTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTTACCGGTGTAATGGGTGCCTGTCTGCCGGAAAAGCTCGGGGCCCAGAATCGTTATCGGGGGctggctgg         >  2:1148385/1‑150 (MQ=255)
                                                                                                                                               aaCCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCTGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCtggtg       <  1:1315208/150‑1 (MQ=255)
                                                                                                                                                 ccGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGc     >  1:1086891/1‑150 (MQ=255)
                                                                                                                                                 ccGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGc     >  2:258577/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTTAAATTGGCGCTCAAAGGTAAGTAAAGGGAGTTTGATATGTCTGTTTCACGTCGGGTAATACATCACGGACTTTATTTTGCAGTTTTAGGACCGTTAATTGGTGTTCTGTTTCTTGTCCTCTACATATTCTTCGCAAAAGAACCGCTGGTTCTTTGGGTGATAATACATCCAATTTTTCTCTTATTGTCGATAACTACGGGAGCTATTCCTGCGTTGTTAACCGGTGTAATGGTTGCCTGTCTGCCGGAAAAGATCGGGTCACAGAAACGTTATCGTTGTCTGGCTGGTGGCATA  >  NC_000913/3897466‑3897763

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