Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,671,548 A→G 5.2% Q58R (CAG→CGG)  asr → acid shock‑inducible periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,671,5480AG 5.2% 245.0 / 3.0 77Q58R (CAG→CGG) asracid shock‑inducible periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (39/34);  new base G (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

GGTCTGTCTTCTGCCGCCTTTGCTGCAGAGACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACA  >  NC_000913/1671415‑1671694
                                                                                                                                     |                                                                                                                                                  
ggTCTGTCTTCTGCCGCCTTTGCTGCAGAGACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAg                                                                                                                                       <  3:386571/147‑1 (MQ=255)
 gtctgtctTCTGCCGCCTTTGCTGCAGAGACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGca                                                                                                                                     <  1:153703/148‑1 (MQ=255)
  tctgtctTCTGCCGCCTTTGCTGCAGAGACTACGACCACACCTGCTCCGACTGCGACGCCCACCAAAGCCGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGcat                                                                                                                                    >  1:346162/1‑148 (MQ=255)
      tcttctGCCGCCTTTGCTGCAGAGACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGcatca                                                                                                                                  <  5:95895/146‑1 (MQ=255)
       cttctGCCGCCTTTGCTGCAGAGACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATa                                                                                                                                <  3:263445/147‑1 (MQ=255)
              cgccTTTGCTGCAGAGACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATAc                                                                                                                         >  7:378786/1‑147 (MQ=255)
                   ttGCTGCAGAGACTACGACCCCACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAACCTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGc                                                                                                                   >  2:255000/1‑148 (MQ=255)
                      ctgcAGAGACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTg                                                                                                                 >  7:104243/1‑147 (MQ=255)
                         cAGAGACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAAc                                                                                                              <  1:314287/147‑1 (MQ=255)
                          agagACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGaa                                                                                                               >  3:158683/1‑145 (MQ=255)
                            agaCTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGaa                                                                                                          <  8:378786/148‑1 (MQ=255)
                              aCTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCACCATAAAAATACGAAAGCTGAACAGAAAg                                                                                                        <  1:340039/148‑1 (MQ=255)
                              aCTACGACCACACCTGATCCGACTGCGACGACCAACAAAGCAGCGCCGTAGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAg                                                                                                        <  7:366668/148‑1 (MQ=255)
                                tACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGcc                                                                                                      >  3:326229/1‑148 (MQ=255)
                                 aCGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGccc                                                                                                     >  1:26901/1‑148 (MQ=255)
                                 aCGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGccc                                                                                                     >  7:252465/1‑148 (MQ=255)
                                  ggACCACACCTGCTCCGACTGCGACGAACACCAAAGCAGCTCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGcc                                                                                                      <  2:26901/145‑1 (MQ=255)
                                  cGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGcc                                                                                                      <  8:252465/146‑1 (MQ=255)
                                    aCCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCt                                                                                                   >  5:288594/1‑147 (MQ=255)
                                     ccACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACCCATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGa                                                                                                 <  1:432288/148‑1 (MQ=255)
                                     ccACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGa                                                                                                 <  1:200686/148‑1 (MQ=255)
                                        cacCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAACCAGCATAAAGCAGAACCTGCCCAGAAAGCGCAGGCGGCTAAAAACCAGCATAAAAATACGAAAGCTGAACAGAAAGCCCCAGAAc                                                                                               >  8:226335/1‑147 (MQ=255)
                                         acCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACa                                                                                              >  3:234921/1‑147 (MQ=255)
                                           cTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAACGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACaaaa                                                                                           <  8:99782/148‑1 (MQ=255)
                                            tGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACaaa                                                                                            >  5:7502/1‑146 (MQ=255)
                                              cTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAgc                                                                                        <  8:148223/148‑1 (MQ=255)
                                                 cGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAgcgc                                                                                      <  4:326229/147‑1 (MQ=255)
                                                   acagCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCCGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGCACAGAAAGCCCCTCAACAAAAAGCgcagg                                                                                   >  8:182948/4‑148 (MQ=255)
                                                        gacgacCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCGGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGc                                                                               >  1:346444/1‑147 (MQ=255)
                                                         acgacCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCga                                                                             <  6:288594/148‑1 (MQ=255)
                                                               accaAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGaaa                                                                        >  6:2104/1‑147 (MQ=255)
                                                                ccaAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAaca                                                                      >  1:55793/1‑148 (MQ=255)
                                                                  aaaGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACg                                                                    <  8:104243/148‑1 (MQ=255)
                                                                         cgcCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCTTAATAATACGAAAGCTGAACAGAAAGCGCCTGAACAAAAAGCGCCGGCAGCGAAGAAACACGCCAAgga                                                             >  2:313528/1‑146 (MQ=255)
                                                                          gcCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCGGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAac                                                                       <  4:343855/137‑1 (MQ=255)
                                                                          gcCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCGGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAac                                                                       >  3:343855/1‑137 (MQ=255)
                                                                              gCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAaca                                                                      >  5:388878/1‑134 (MQ=255)
                                                                              gCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAaca                                                                      <  6:388878/134‑1 (MQ=255)
                                                                              gCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAaca                                                          >  7:380666/1‑146 (MQ=255)
                                                                                   aaCTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCAt                                                   <  6:7502/148‑1 (MQ=255)
                                                                                    aCTACACATCATAAAAAACAGCATAAAGCATCACCTGCCCAGAAAGCGCAGGCGTCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCTGGCAGCTAAGAAACACGCCAAGAAACACAGCCATc                                                  <  1:355842/148‑1 (MQ=255)
                                                                                    aCTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCGGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATc                                                  <  2:346444/148‑1 (MQ=255)
                                                                                       acacATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAgcgc                                                                                      <  8:320474/109‑1 (MQ=255)
                                                                                       acacATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAgcgc                                                                                      >  7:320474/1‑109 (MQ=255)
                                                                                       acacATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCgaa                                                                            >  7:152659/1‑119 (MQ=255)
                                                                                          catcatAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGATAAACGCCAAGAAACACAGCCagacgcaa                                             >  6:228107/1‑140 (MQ=255)
                                                                                            tcatAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCg                                          <  3:311219/148‑1 (MQ=255)
                                                                                             catAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAgcgc                                                                                      >  3:186519/1‑103 (MQ=255)
                                                                                             catAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAgcgc                                                                                      <  4:186519/103‑1 (MQ=255)
                                                                                             catAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCgaagaa                                                                         >  7:123860/1‑116 (MQ=255)
                                                                                             catAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCgaagaa                                                                         <  8:123860/116‑1 (MQ=255)
                                                                                             catAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCgg                                         >  4:302544/1‑148 (MQ=255)
                                                                                              atAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGaaa                                                                                                                     >  5:242061/1‑71 (MQ=255)
                                                                                              atAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGaaa                                                                                                                     <  6:242061/71‑1 (MQ=255)
                                                                                               tAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATc                                                  <  8:352564/137‑1 (MQ=255)
                                                                                               tAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATc                                                  >  7:352564/1‑137 (MQ=255)
                                                                                                 aaaaaCAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCaa                                      >  3:290168/1‑147 (MQ=255)
                                                                                                 aaaaaCAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCaa                                      >  3:3639/1‑147 (MQ=255)
                                                                                                 aaaaaCAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCaa                                      >  3:295762/1‑147 (MQ=255)
                                                                                                  aaaaCAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAAGACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACCCCAAGAAACACAGCCATAAGCAACCGGCaaaa                                    >  5:174194/1‑148 (MQ=255)
                                                                                                       aGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACctgc                               <  1:255000/148‑1 (MQ=255)
                                                                                                        gCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACctgct                              <  8:380666/148‑1 (MQ=255)
                                                                                                             aaGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGcac                          >  4:251671/1‑147 (MQ=255)
                                                                                                                 agcaCCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAAccc                     >  4:123031/1‑148 (MQ=255)
                                                                                                                  gcaCCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAacac                                                         >  3:22133/1‑111 (MQ=255)
                                                                                                                  gcaCCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAacac                                                         <  4:22133/111‑1 (MQ=255)
                                                                                                                      cTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGc                >  8:434928/1‑148 (MQ=255)
                                                                                                                        gCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGa                                                                                                           >  6:221895/1‑55 (MQ=255)
                                                                                                                        gCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGa                                                                                                           <  5:221895/55‑1 (MQ=255)
                                                                                                                         cccAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTa             <  5:2104/148‑1 (MQ=255)
                                                                                                                            aGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGt          >  3:166841/1‑148 (MQ=255)
                                                                                                                            aGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGt          >  5:218678/1‑148 (MQ=255)
                                                                                                                             gAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCTCAACCCGCAGCGTAAGtt         >  4:234177/1‑148 (MQ=255)
                                                                                                                               aaGCGCAGGCGGCTAAAAAGCATCTTAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGtttt       >  4:407622/1‑148 (MQ=255)
                                                                                                                                aGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTc      <  3:251671/148‑1 (MQ=255)
                                                                                                                                  cgcAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCaa    <  5:228107/148‑1 (MQ=255)
                                                                                                                                    cAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACa  <  4:234921/148‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                  
GGTCTGTCTTCTGCCGCCTTTGCTGCAGAGACTACGACCACACCTGCTCCGACTGCGACGACCACCAAAGCAGCGCCGGCGAAAACTACACATCATAAAAAACAGCATAAAGCAGCACCTGCCCAGAAAGCGCAGGCGGCTAAAAAGCATCATAAAAATACGAAAGCTGAACAGAAAGCCCCTGAACAAAAAGCGCAGGCAGCGAAGAAACACGCCAAGAAACACAGCCATCAGCAACCGGCAAAACCTGCTGCACAACCCGCAGCGTAAGTTTTCAACA  >  NC_000913/1671415‑1671694

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: ‑