Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,314,886 T→C 5.2% S113S (TCA→TCG)  yciE ← putative rubrerythrin/ferritin‑like metal‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,314,8860TC 5.2% 259.5 / 3.5 77S113S (TCA→TCG) yciEputative rubrerythrin/ferritin‑like metal‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base C (2/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATG  >  NC_000913/1314738‑1315033
                                                                                                                                                    |                                                                                                                                                   
cATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGt                                                                                                                                                     >  3:133312/1‑149 (MQ=255)
 aTCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTg                                                                                                                                                    >  1:392513/1‑149 (MQ=255)
     gTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGgt                                                                                                                                                >  2:89929/1‑149 (MQ=255)
        tCAGAGCTAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTa                                                                                                                                             <  4:321350/149‑1 (MQ=255)
        tCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTa                                                                                                                                             >  2:508380/1‑149 (MQ=255)
          agagCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAAc                                                                                                                                           <  1:506934/149‑1 (MQ=255)
           gagCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACa                                                                                                                                          >  4:156496/1‑149 (MQ=255)
             gCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAgg                                                                                                                                        >  2:173101/1‑149 (MQ=255)
                aaTTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGGGGTGTAACAGGCGa                                                                                                                                     <  3:142366/149‑1 (MQ=255)
                aaTTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGa                                                                                                                                     <  3:280074/149‑1 (MQ=255)
                aaTTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGa                                                                                                                                     <  3:229374/149‑1 (MQ=255)
                aaTTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGa                                                                                                                                     <  4:133312/149‑1 (MQ=255)
                aaTTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGa                                                                                                                                     >  4:609079/1‑149 (MQ=255)
                 aTTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGAt                                                                                                                                    <  4:60273/149‑1 (MQ=255)
                   tAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGAttt                                                                                                                                  <  2:615814/149‑1 (MQ=255)
                    aaaaaTTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTc                                                                                                                                 >  4:536740/1‑149 (MQ=255)
                     aaaaTTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCa                                                                                                                                >  2:123930/1‑149 (MQ=255)
                      aaaTTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGgt                                                                                                                                                >  3:611437/1‑132 (MQ=255)
                      aaaTTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGgt                                                                                                                                                <  4:611437/132‑1 (MQ=255)
                            ctcAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCt                                                                                                                         >  3:412976/1‑149 (MQ=255)
                             tcAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGgt                                                                                                                                                <  3:206725/125‑1 (MQ=255)
                             tcAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGgt                                                                                                                                                >  4:206725/1‑125 (MQ=255)
                             tcAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTc                                                                                                                        <  2:45806/149‑1 (MQ=255)
                                         ggAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATcc                                                                                                            >  1:687168/1‑149 (MQ=255)
                                          gAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCg                                                                                                           >  2:180574/1‑149 (MQ=255)
                                                  tGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATaga                                                                                                   <  4:412976/149‑1 (MQ=255)
                                                         gCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATa                                                                                                     >  2:491429/1‑140 (MQ=255)
                                                         gCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATa                                                                                                     <  1:491429/140‑1 (MQ=255)
                                                         gCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTg                                                                                            >  1:214943/1‑149 (MQ=255)
                                                          ccAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGa                                                                                           <  2:346010/149‑1 (MQ=255)
                                                           cAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGAc                                                                                          >  2:27658/1‑149 (MQ=255)
                                                               cGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTAtt                                                                                      <  4:79695/149‑1 (MQ=255)
                                                                ggCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGCGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTAttt                                                                                     >  1:350418/1‑149 (MQ=255)
                                                                ggCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGCGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTAttt                                                                                     >  3:149203/1‑149 (MQ=255)
                                                                  ccATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTtca                                                                                   >  4:13655/1‑149 (MQ=255)
                                                                   cATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTtcat                                                                                  <  2:540147/149‑1 (MQ=255)
                                                                     tttGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTtcatca                                                                                >  3:557980/1‑149 (MQ=255)
                                                                      ttGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAg                                                                               >  1:187213/1‑149 (MQ=255)
                                                                        gCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGaa                                                                             >  4:126174/1‑149 (MQ=255)
                                                                        gCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGaa                                                                             <  3:13655/149‑1 (MQ=255)
                                                                        gCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGaa                                                                             <  1:180574/149‑1 (MQ=255)
                                                                        gCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGaa                                                                             <  4:670178/149‑1 (MQ=255)
                                                                         cTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAg                                                                            >  2:20899/1‑149 (MQ=255)
                                                                         cTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAg                                                                            <  1:173101/149‑1 (MQ=255)
                                                                             tCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGaa                                                                        <  1:648616/149‑1 (MQ=255)
                                                                                tCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAatat                                                                     >  4:653669/1‑149 (MQ=255)
                                                                                tCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAatat                                                                     >  3:373733/1‑149 (MQ=255)
                                                                                  aTTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATac                                                                   >  4:120201/1‑149 (MQ=255)
                                                                                   tttAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAg                                                                                                                                                  >  3:472635/1‑69 (MQ=255)
                                                                                   tttAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAg                                                                                                                                                  <  4:472635/69‑1 (MQ=255)
                                                                                    ttAAAATCTCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATacca                                                                 <  2:563590/149‑1 (MQ=255)
                                                                                    ttAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATacca                                                                 <  3:536740/149‑1 (MQ=255)
                                                                                          tCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGAtt                                                           >  3:607843/1‑149 (MQ=255)
                                                                                              ttCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATc                                                                                                             >  2:12801/1‑95 (MQ=255)
                                                                                              ttCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATc                                                                                                             <  1:12801/95‑1 (MQ=255)
                                                                                              ttCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACt                                                       <  3:609079/149‑1 (MQ=255)
                                                                                              ttCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACt                                                       >  3:167904/1‑149 (MQ=255)
                                                                                                    ggTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGCGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCaa                                                 <  2:350418/149‑1 (MQ=255)
                                                                                                    ggTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGCGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCaa                                                 <  4:149203/149‑1 (MQ=255)
                                                                                                           aTTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGc                                          <  2:392513/149‑1 (MQ=255)
                                                                                                                  cTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTa                                   >  2:129039/1‑149 (MQ=255)
                                                                                                                   tGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTAc                                  <  4:557980/149‑1 (MQ=255)
                                                                                                                   tGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTAc                                  >  4:642081/1‑149 (MQ=255)
                                                                                                                    gTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGccc                                                   <  4:127801/131‑1 (MQ=255)
                                                                                                                    gTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGccc                                                   >  3:127801/1‑131 (MQ=255)
                                                                                                                     tATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTc                                <  1:123930/149‑1 (MQ=255)
                                                                                                                       tCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCAt                              <  3:653669/149‑1 (MQ=255)
                                                                                                                       tCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCAt                              <  1:508380/149‑1 (MQ=255)
                                                                                                                       tCACCGGCATTGTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCAt                              >  3:403920/1‑149 (MQ=255)
                                                                                                                                ttttttGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCttt                     <  1:20899/149‑1 (MQ=255)
                                                                                                                                  ttttGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTaa                   <  3:156496/149‑1 (MQ=255)
                                                                                                                                      gctgctgctAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGAc               >  4:297656/1‑149 (MQ=255)
                                                                                                                                           tgctAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGAtt                                                           <  2:104919/100‑1 (MQ=255)
                                                                                                                                           tgctAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGAtt                                                           >  1:104919/1‑100 (MQ=255)
                                                                                                                                            gctAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACg         >  4:442468/1‑149 (MQ=255)
                                                                                                                                             ctAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGt        <  4:384286/149‑1 (MQ=255)
                                                                                                                                                 tAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGaaa    <  1:89929/149‑1 (MQ=255)
                                                                                                                                                   gTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATg  >  4:387745/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATCAGTTTCAGAGCGAATTAAAAATTTCTCAGTTGTTTGCGGAATATTCTGAATCAGCCAGTCGGCCATTTGCTTTTCCTCATTTAAAATCGCTTCGATGGTTGGAATTGAAGCTGTATCACCGGCATTTTTTGCTGCTGCTAATAGTGAGGTGTAACAGGCGATTTCAAATTGCTCGAAGACATATCCGCTAATAGAGCCTTTGACTATTTCATCAGAAGGGAATATACCACCGATTGACTGCCCAAGCGCAGCCATTTTACTCATGGAATCTTTAATGACTGAACGTGAAATG  >  NC_000913/1314738‑1315033

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑