Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,698,203 G→T 7.4% A60S (GCA→TCA)  bcsF → DUF2636 family cellulose production small membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,698,2030GT 7.4% 222.2 / 3.8 71A60S (GCA→TCA) bcsFDUF2636 family cellulose production small membrane protein
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base T (2/4);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCT  >  NC_000913/3698056‑3698347
                                                                                                                                                   |                                                                                                                                                
attGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAggg                                                                                                                                                  >  7:103351/1‑148 (MQ=255)
attGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAggg                                                                                                                                                  <  5:367233/148‑1 (MQ=255)
            gcACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATg                                                                                                                                     <  4:67760/149‑1 (MQ=255)
             cACTGATATTTTTCCCTCTGGTCTATCTGTAGCGGCTCTCTTTGCGACGCATGCGCGACTCCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTCAGCCGCACGGAAAAAGCCAGGGCAACCTAAACATGa                                                                                                                                    <  4:116792/149‑1 (MQ=255)
             cACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGa                                                                                                                                    <  1:193953/149‑1 (MQ=255)
               cTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACt                                                                                                                                  <  1:14494/149‑1 (MQ=255)
                tGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAACCCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACt                                                                                                                                  >  5:110874/1‑148 (MQ=255)
                tGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACt                                                                                                                                  >  4:141142/1‑148 (MQ=255)
                  atatTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCaa                                                                                                                               >  1:13614/1‑149 (MQ=255)
                     tttttCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAAtt                                                                                                                             <  2:167993/148‑1 (MQ=255)
                               tGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAggg                                                                                                                                                  <  1:294852/117‑1 (MQ=255)
                               tGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAggg                                                                                                                                                  >  2:294852/1‑117 (MQ=255)
                                   gtATCTTGCGCGGCACTATTTGCGAGGCATTCGCGATACCTTACGTTTTTTCTTTGCTAAACCTCGTTATTTTAATCCGGCCGGGACGTTACTCCTCACGGAAAAAGCCAGTTCAACCAAAAAATGACTCAATTTTCGCAAATTATCGc                                                                                                              <  2:273311/148‑1 (MQ=255)
                                   cTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTAAGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCTAAAAATGACTCAATTTACGCAAAATACCGc                                                                                                              <  5:108690/149‑1 (MQ=255)
                                      tCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGcc                                                                                                             <  5:70119/147‑1 (MQ=255)
                                      tCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCa                                                                                                            <  3:373633/148‑1 (MQ=255)
                                      tCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCAt                                                                                                           >  5:320934/1‑149 (MQ=255)
                                            gcgGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCcttc                                                                                                     >  1:219603/1‑149 (MQ=255)
                                              ggCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCcttctt                                                                                                   <  2:13614/149‑1 (MQ=255)
                                                  ctctTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCaaaa                                                                                                                    <  5:171277/128‑1 (MQ=255)
                                                  ctctTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCaaaa                                                                                                                    >  6:171277/1‑128 (MQ=255)
                                                   tctTTGCGACGCATTCGCGAAACCTTACGTTTGTTCTTTGCTAAAACTCGTTATGTTAACCCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTAACGCAAAATACCGCCATGACTTCTTACctc                                                                                              >  6:255580/1‑149 (MQ=255)
                                                     tttGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCctc                                                                                              <  7:181108/147‑1 (MQ=255)
                                                         cGACGCATTCGCGACAGCTTACGTTTGTTCTTTGCTACAGCTCGTTATGTTAAACCGGCGGGGACGTTACGCCGCAGGGATATAGCCAGGGCAACCAAATATTGACTCAATTTAGGCTTAATACCGGCATGCCTTCTTCCCTCTGGCaa                                                                                        <  7:170761/149‑1 (MQ=255)
                                                         cGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCaa                                                                                        <  8:103351/149‑1 (MQ=255)
                                                          gACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTATGCAAAATACCGCCATGCCTTCTTCCCTCTGGCaa                                                                                        >  8:386733/1‑148 (MQ=255)
                                                            cGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATAc                                                                                     >  6:347332/1‑149 (MQ=255)
                                                             gCATTCGCGACACCTTCCTGTTGTTCTTTGCTTAACCTCTTTATGTTAAACCGGCCGTGACGTTACGCCGGACGGAAAAAGCCATGGCAACCAAAAACTGACTCAATTTACTCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACt                                                                                    <  2:64670/149‑1 (MQ=255)
                                                             gCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAACCCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTAAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACt                                                                                    >  1:366695/1‑149 (MQ=255)
                                                               aTTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGTCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACt                                                                                    >  8:7529/1‑147 (MQ=255)
                                                               aTTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGTCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACt                                                                                    <  7:7529/147‑1 (MQ=255)
                                                               aTTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTgg                                                                                  >  8:151873/1‑149 (MQ=255)
                                                                           tttcGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCAGGGACGTTACGTCGCACGGAAAAAGCCAGGTCAACCAAAAATTGACTAAATTTTCGCAAAATTCAGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTcc                                                                      <  6:290316/146‑1 (MQ=255)
                                                                           ttACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTcc                                                                                                 >  7:418644/1‑122 (MQ=255)
                                                                           ttACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTcc                                                                                                 <  8:418644/122‑1 (MQ=255)
                                                                             aCGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCgg                                                                    >  1:114681/1‑149 (MQ=255)
                                                                                 ttGCTCTTTGCTCAACCTCGTTATGTTAAACCGGCCGGTACGGTACGCCGCACGTAAAAATCCAGCGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTgg                                                                <  6:110874/149‑1 (MQ=255)
                                                                                    ttCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAAc                                                             >  1:43147/1‑149 (MQ=255)
                                                                                     tCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAACTACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACt                                                            >  6:20179/1‑149 (MQ=255)
                                                                                     tCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACt                                                            >  1:66304/1‑149 (MQ=255)
                                                                                       tttGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGTCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTc                                                          >  7:291281/1‑149 (MQ=255)
                                                                                       tttGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCAAGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTc                                                          >  1:155063/1‑149 (MQ=255)
                                                                                       tttGCTAAACCTCGTGATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCATTACTGGCGCGGCCTTTCCGGCTGGAGCTTc                                                          >  6:209883/1‑149 (MQ=255)
                                                                                          gCTAAACCTCGATATCTTAAACCCGCCGCGACGTTACCCCGCACGGAAAAAGCCAGGCCCACCAAAAAATGACTCAATTTACGCAAAATCCCGCCATGCCTTCTTCCCTCTTGCAATACTGGCACGGCCTTTCCGGCTGGAACATCAAt                                                       >  4:209295/1‑149 (MQ=255)
                                                                                               aCCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCt                                                  <  2:219603/149‑1 (MQ=255)
                                                                                                ccTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTAttt                                                     >  3:27229/1‑145 (MQ=255)
                                                                                                ccTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTAttt                                                     <  4:27229/145‑1 (MQ=255)
                                                                                                ccTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTg                                                 >  5:83993/1‑149 (MQ=255)
                                                                                                    gTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGtt                                              <  7:151873/148‑1 (MQ=255)
                                                                                                        tGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCgg                                                                             >  6:268130/1‑113 (MQ=255)
                                                                                                        tGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCgg                                                                             <  5:268130/113‑1 (MQ=255)
                                                                                                          ttAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCg                                       >  1:418490/1‑149 (MQ=255)
                                                                                                          ttAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCg                                       >  6:318079/1‑149 (MQ=255)
                                                                                                             aaCCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGcc                                    >  1:347504/1‑149 (MQ=255)
                                                                                                             aaCCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGc                                     <  7:351779/148‑1 (MQ=255)
                                                                                                                 ggccggGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGAGTCAATTTACGCAAAATAGCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCtgtt                                <  5:255580/149‑1 (MQ=255)
                                                                                                                 ggccggGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCtgtt                                <  5:347332/149‑1 (MQ=255)
                                                                                                                   ccggGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCtgttgt                              >  1:235211/1‑149 (MQ=255)
                                                                                                                      ggACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTggg                           <  6:83993/149‑1 (MQ=255)
                                                                                                                          gTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCggg                       <  2:418490/149‑1 (MQ=255)
                                                                                                                           ttACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGa                      >  7:179990/1‑149 (MQ=255)
                                                                                                                           ttACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGa                      >  6:164716/1‑149 (MQ=255)
                                                                                                                            tACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCACTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGGTGTGGGCGGGa                      >  8:255635/1‑148 (MQ=255)
                                                                                                                                 cgcACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCtt                <  2:43147/149‑1 (MQ=255)
                                                                                                                                  gcACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGTCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTa               <  2:347504/149‑1 (MQ=255)
                                                                                                                                  gcACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGGACTCCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCtt                >  7:33518/1‑148 (MQ=255)
                                                                                                                                  gcACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGGACTCCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCtt                >  7:27119/1‑148 (MQ=255)
                                                                                                                                    aCGGAAAAAGCCAGGTCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAAc             <  8:291281/149‑1 (MQ=255)
                                                                                                                                        aaaaaGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAAGTTc          >  6:47002/1‑148 (MQ=255)
                                                                                                                                           aaGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATc      <  6:320934/149‑1 (MQ=255)
                                                                                                                                           aaGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATc      <  8:179990/149‑1 (MQ=255)
                                                                                                                                           aaGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATc      >  1:245336/1‑149 (MQ=255)
                                                                                                                                              ccAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGc   >  6:90547/1‑149 (MQ=255)
                                                                                                                                              ccAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGc   >  8:206338/1‑149 (MQ=255)
                                                                                                                                               cAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGct  >  8:164981/1‑149 (MQ=255)
                                                                                                                                                  ggCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGc   <  6:162773/145‑1 (MQ=255)
                                                                                                                                                  ggCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGc   >  5:162773/1‑145 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
ATTGTCGTTTGCGCACTGATATTTTTCCCGCTGGGCTATCTGGCGCGGCACTCTTTGCGACGCATTCGCGACACCTTACGTTTGTTCTTTGCTAAACCTCGTTATGTTAAACCGGCCGGGACGTTACGCCGCACGGAAAAAGCCAGGGCAACCAAAAAATGACTCAATTTACGCAAAATACCGCCATGCCTTCTTCCCTCTGGCAATACTGGCGCGGCCTTTCCGGCTGGAACTTCTATTTTCTGGTTAAGTTCGGCCTGTTGTGGGCGGGATATCTTAACTTCCATCCGCT  >  NC_000913/3698056‑3698347

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: