Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,528,908 T→A 12.9% I81N (ATC→AAC)  yrfG → GMP/IMP nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,528,9080TA12.9% 191.0 / 4.5 88I81N (ATC→AAC) yrfGGMP/IMP nucleotidase
Reads supporting (aligned to +/- strand):  ref base T (39/35);  new base A (2/10);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 2.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAAAGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGAG  >  NC_000913/3528760‑3529050
                                                                                                                                                     |                                                                                                                                              
cAAAAGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGatat                                                                                                                                                >  1:839891/1‑149 (MQ=255)
cAAAAGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGAGGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGatat                                                                                                                                                >  1:265042/1‑149 (MQ=255)
 aaaaGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGAGGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATc                                                                                                                                               >  2:607212/1‑149 (MQ=255)
 aaaaGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGCTTACTGGAGTGAGCAACTGGGTCTGGATATc                                                                                                                                               >  2:380965/1‑149 (MQ=255)
 aaaaGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATc                                                                                                                                               >  1:1328767/1‑149 (MQ=255)
 aaaaGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGGGAGCAACTGGGTCTGGATGTc                                                                                                                                               >  1:101103/1‑149 (MQ=255)
 aaaaGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGGCTTGATTACTGGAGTGAGCAACTGGGTCTGGAAATc                                                                                                                                               >  1:724733/1‑149 (MQ=255)
  aaaGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCt                                                                                                                                              <  1:761921/149‑1 (MQ=255)
     gCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTACAAGTAGT‑ATTACTGGAGTGAGCAACTGGGTCTGGATATCtgtg                                                                                                                                           <  2:1240001/149‑1 (MQ=255)
       tGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAG‑TGTCTTGATTACGTGAGTGAGCCACTGGGTCTGGATATCTGTGCg                                                                                                                                         >  1:587535/1‑149 (MQ=255)
         gTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCCAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATAACTGTGCGAt                                                                                                                                       >  2:383403/1‑149 (MQ=255)
             cTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGGGCAACTGGGTCTGGATAGCTGTGCGATGgcc                                                                                                                                   >  1:1012371/1‑149 (MQ=255)
                aaaCATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGaccacc                                                                                                                                >  1:227010/1‑149 (MQ=255)
                aaaCATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACCTACAGCATACGCTAAACCAGTAC‑CGGCTTGAATACTGGGGTGAGCAACCGGGTCTCGAGAACTGTGCGAttctaccc                                                                                                                                >  2:506532/1‑142 (MQ=255)
                       gggCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATAGCTGTGCGATGACCACCGAGATgg                                                                                                                         >  1:147898/1‑149 (MQ=255)
                        ggCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGAGCCAGCAATGTGACAACGTACAGCATACACTAACGTGGTAC‑TGTCTAAATTACTGGAGAGAACGACTGGGTCTGGACATGTGGACGATGACGACCGACATgtg                                                                                                                        >  2:671460/1‑147 (MQ=255)
                         gcgcGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCACCCATAGTATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGa                                                                                                                       <  1:1049908/149‑1 (MQ=255)
                          ggcgAAATATGTGGTTAGGCCACAGGAAGCGAGGGAATATGTGCGCCAGCATTATCGCGACGTACAGCATACGCTGTTCTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGTTTTTTTTTAATGATGCGGCGACCACCGAGATGGGAc                                                                                                                      <  2:1712/148‑1 (MQ=255)
                            cgAAAAACTGGATTACGCCACAGGAAGCGATGGAATATATGCGCCATCAATATCACGACTTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTTTTTTTTTTTATTTATACGACGACCACCGAGATGGGACcg                                                                                                                    <  2:249676/149‑1 (MQ=255)
                            cgAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACAGGGACTGGATATCTGTGCGATGACCACCGAGATGGGACcg                                                                                                                    >  2:495502/1‑149 (MQ=255)
                             gAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACcgc                                                                                                                   <  1:243792/149‑1 (MQ=255)
                              aaaaaCGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACcgcg                                                                                                                  >  1:930654/1‑149 (MQ=255)
                              aaaaaCGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACAGGGGATGGATGTCTGTGCGATGACGACCGGGATGGGACcgag                                                                                                                  >  2:1106290/1‑147 (MQ=255)
                               aaaaCGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACCGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGt                                                                                                                 <  1:1195410/149‑1 (MQ=255)
                                  aCGGTGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGTGCAATTGTTTTTTTTTAATGATGCGGCGACCACCGAGATGGGACCGCGTGcc                                                                                                              <  2:1032326/149‑1 (MQ=255)
                                  aCGGGGTTACCCCACAGGAAGCGATGGAATATATGCGCCAGCCATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTTTGTTTTTATAACGATACGGCGACCACCGAGATGGGACCGCGTGcc                                                                                                              <  2:1306451/149‑1 (MQ=255)
                                   cGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTGC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACTACCGAGATGGGACCGCGTGCCg                                                                                                             >  1:1208975/1‑149 (MQ=255)
                                     gggTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTa                                                                                                           <  1:1035660/149‑1 (MQ=255)
                                        ttACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAGACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTg                                                                                                        >  2:1130554/1‑149 (MQ=255)
                                         tACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATTTGTGCGGTGACCACCGAGATGGGACCGCGTGCCGTACTGc                                                                                                       <  2:839891/149‑1 (MQ=255)
                                           cGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGt                                                                                                     >  1:1346912/1‑149 (MQ=255)
                                             ccACGGGAAGCGATGCAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGTTTTTTTTATTGGTGCGGCGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGa                                                                                                   <  2:1342238/149‑1 (MQ=255)
                                             ccACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGTGTTTTTTTAATGATACGGTGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGa                                                                                                   <  2:663717/149‑1 (MQ=255)
                                             ccACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGa                                                                                                   >  1:732983/1‑149 (MQ=255)
                                              cacaGGAAGCGATGGAATATATGCGCCAGCAATATCACGCCGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGGGGGAGCAACTGGGCCTGGATATCTGTGCGGTGACCACCGAGCTGGGACCGCGTACCGTTCTACGTCaa                                                                                                  >  2:92611/1‑149 (MQ=255)
                                                     aGCGATGGAATATATGCGCCAGCCCTATCGCGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGGGAGCAACTGGGACTGGATATCTGTGCGCTGACCACCGAGAAGGGACCGCGTGCCGTACTGCCTGAAGATACCa                                                                                           >  2:491915/1‑149 (MQ=255)
                                                      gCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCAt                                                                                          >  1:1211104/1‑149 (MQ=255)
                                                        taTGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGTTTCTTTATAACGATACGACGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTc                                                                                        <  2:1223422/148‑1 (MQ=255)
                                                         aTGGAATATATGCGCCAGCGATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTTTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTcc                                                                                       >  1:982926/1‑149 (MQ=255)
                                                          tGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGGCCGCGGGCCGTACTGCGTGAAGATACCATTCCg                                                                                      >  1:438568/1‑149 (MQ=255)
                                                          tGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGGGAAGCTTCCATTCCg                                                                                      >  1:239232/1‑149 (MQ=255)
                                                          tGGAATATATGCGCCAGCAATATCACGACATACAGCATACGCTAAACTGGTAC‑TGTCTTGAATACTGAAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGTCGTACGGCGGGAAGATACCATTCCg                                                                                      >  1:828149/1‑149 (MQ=255)
                                                           ggAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGGCCACCAAGATGGGACCGCGGGCCGTCCTGCGCGAAGAGACCATTCCGt                                                                                     >  2:569753/1‑149 (MQ=255)
                                                            gagcatctGCGCCAGCAATAGCACGACGTACAGCATACGCTAACCTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTTTGGTTAACGGTGCGGCGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGtt                                                                                    <  2:1006989/142‑1 (MQ=255)
                                                            gAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGtt                                                                                    >  1:87549/1‑149 (MQ=255)
                                                              atatatGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTAGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGGGCCGTACTGCGTGAACATACAATTCGGTTTc                                                                                  >  2:290731/1‑149 (MQ=255)
                                                               tatatGCGCCAGCAATATCACGACGTACGGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGGGAGCAACTGGGTCTGGATATCTGTGCGATGACAACCGAGATGGGAGCGCGTGCCGTACTGCGTGAAGATACCATTCCGTGTCt                                                                                 >  1:882407/1‑149 (MQ=255)
                                                                aactGCGGCAGCAATATCACGATGTACAGCATACGCCAAACTTGTAC‑TGTCTTGATTACTGGAGTGAGCAACTTTTTTTTTTTATTGGTGCGGCGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCtt                                                                                <  2:151694/146‑1 (MQ=255)
                                                                 tatGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTg                                                                               >  1:1128033/1‑149 (MQ=255)
                                                                 tatGCGCCAGCAATATCACGACGTACAGCATACGCGAAACTGGGAC‑TGTCTTGATTACTGCAGTGAGCAACTGGGTCTGGAGATCTGTGCGATCACCACCGAGATGGGACCGCGTGCCGTACTGCGGGAAGAGACCGTTCCGTTTCTTg                                                                               >  1:633741/1‑149 (MQ=255)
                                                                 catGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTg                                                                               <  1:1229943/148‑1 (MQ=255)
                                                                  atGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGGTTCTTGa                                                                              >  1:286609/1‑149 (MQ=255)
                                                                  atGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGAAGTGAGCAACTGGGTCTGGCTATATGTTCGATGCCCACCCGGGGGGGACCCGGTGCGGTACTGCGGGACGATACCACTACGTTTCtggg                                                                              >  1:182540/1‑146 (MQ=255)
                                                                   tGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAg                                                                             <  1:1147428/149‑1 (MQ=255)
                                                                   tGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAg                                                                             <  2:1246712/149‑1 (MQ=255)
                                                                       ccAGCAATAGCACGTCGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCAc                                                                         <  1:103351/149‑1 (MQ=255)
                                                                        cAGCAATATCACTACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTTTTTAACGGTGCGGCGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACt                                                                        <  2:1230933/149‑1 (MQ=255)
                                                                          gCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGCGTGAGCCATTGGGTCTGGATGTCAATTTCATGGCCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGa                                                                      >  1:1331472/1‑149 (MQ=255)
                                                                           cAATATCCCGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGaa                                                                     <  1:491915/149‑1 (MQ=255)
                                                                            aaTATCAGGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCGGTGCGACGGCCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGaaa                                                                    <  2:198456/149‑1 (MQ=255)
                                                                            aaTATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTCAGCAACTGTGTCTTGTTATTGATGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGaaa                                                                    <  2:578346/149‑1 (MQ=255)
                                                                            aaTATCACGACGCACCGCGTACGCTAAACTGGTAC‑TGTCTTGAGTACTGGAGGGAGCAACTGTGTCTTGTTAATGGTGCGGCGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGaaa                                                                    <  2:667655/149‑1 (MQ=255)
                                                                                 caccacGTACTGCATACGCTAAACTGGTAC‑TGTCGTGATTACTGGAGTGAGCAACTGGGTCTTGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAg                                                               <  2:930654/149‑1 (MQ=255)
                                                                                 cACGACGTACCGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGTGTTTGTTTAACGATACGGCGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAg                                                               <  2:488969/149‑1 (MQ=255)
                                                                                     acgTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTCGAGGCACTGAAAGCCAGCGGt                                                           <  1:707442/149‑1 (MQ=255)
                                                                                       gTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCAAGGAATCGTTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTaa                                                         <  1:984548/149‑1 (MQ=255)
                                                                                         aCAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTTGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAAATACCATTCCGTTTCTTGAGGCACCGAAAGCCGGCGGTAagc                                                       >  2:369657/1‑149 (MQ=255)
                                                                                              aTACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAAATGCGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGa                                                  <  1:1261826/149‑1 (MQ=255)
                                                                                                   cTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAAtttt                                             >  1:389893/1‑149 (MQ=255)
                                                                                                    tAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGGTGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTg                                            <  2:87549/149‑1 (MQ=255)
                                                                                                    gaaaCTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTg                                            <  2:1331472/148‑1 (MQ=255)
                                                                                                      aaCTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCt                                          <  1:607212/149‑1 (MQ=255)
                                                                                                        cTGGTCC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCa                                        <  1:290440/149‑1 (MQ=255)
                                                                                                             gc‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGTTATCTGTGCGACGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAAt                                   <  2:147898/148‑1 (MQ=255)
                                                                                                             aC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAAt                                   >  1:28722/1‑149 (MQ=255)
                                                                                                                tGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATgc                                 >  1:1205713/1‑149 (MQ=255)
                                                                                                                tGTCGTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATgc                                 <  1:404047/149‑1 (MQ=255)
                                                                                                                ggTCTTGATTACTGGAGTGAGCAACTGGGTCTGGTTATCGGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCCTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATgc                                 <  2:724733/148‑1 (MQ=255)
                                                                                                                                  taGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGc               <  1:1106290/148‑1 (MQ=255)
                                                                                                                                   aGCAACTGGGTCTTGATAACTGTGCGACGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCg              <  2:227010/149‑1 (MQ=255)
                                                                                                                                   aGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGGAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCg              <  1:335318/149‑1 (MQ=255)
                                                                                                                                   aGCAACTGGGTCCGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCg              <  1:495502/149‑1 (MQ=255)
                                                                                                                                        cTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTaaa         <  2:982926/149‑1 (MQ=255)
                                                                                                                                         tGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTaaaa        <  1:1091946/149‑1 (MQ=255)
                                                                                                                                         tGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTaaaa        <  1:1130554/149‑1 (MQ=255)
                                                                                                                                             tCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTg    >  1:76459/1‑149 (MQ=255)
                                                                                                                                             tCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTg    >  1:671000/1‑149 (MQ=255)
                                                                                                                                             tCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTg    >  1:28212/1‑149 (MQ=255)
                                                                                                                                              ctttatatCTGGTCGACGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAATCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGa   <  2:438568/145‑1 (MQ=255)
                                                                                                                                              cTGGTTATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGa   <  2:1001635/149‑1 (MQ=255)
                                                                                                                                               tGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGATTCTTGAGGCCCTGAAAGCCAGCGGTAAGCAGCGGATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGGAGAACTTCAg  >  2:1211505/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
CAAAAGCTGGTGCCTGAAACATGGGGCGCGAAAAACGGGGTTACGCCACAGGAAGCGATGGAATATATGCGCCAGCAATATCACGACGTACAGCATACGCTAAACTGGTAC‑TGTCTTGATTACTGGAGTGAGCAACTGGGTCTGGATATCTGTGCGATGACCACCGAGATGGGACCGCGTGCCGTACTGCGTGAAGATACCATTCCGTTTCTTGAGGCACTGAAAGCCAGCGGTAAGCAGCGAATTTTGCTCACCAATGCGCATCCGCACAACCTGGCGGTAAAACTTGAG  >  NC_000913/3528760‑3529050

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑