Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,358,010 C→T 17.2% P1096S (CCG→TCG)  gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,358,0100CT17.2% 151.6 / 22.2 64P1096S (CCG→TCG) gltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base C (29/24);  new base T (5/6);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

CTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCC  >  NC_000913/3357869‑3358157
                                                                                                                                             |                                                                                                                                                   
cTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCCCTGGCCCGAtgg                                                                                                                                               >  7:255929/1‑148 (MQ=255)
      gAAACCCAGCAGGCGCTGGTTGCTACCGGCTTGCGTCACAAGCTCCGTTTGCAGGTCGAGGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGAGCATAAAGATTCGGCTTAGTCACTGGCCCGATGGTGgcgc                                                                                                                                         >  6:180722/1‑148 (MQ=255)
       aaaCCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCt                                                                                                                                        >  2:126148/1‑148 (MQ=255)
       aaaCCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCt                                                                                                                                        >  8:22908/1‑148 (MQ=255)
        aaCCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTc                                                                                                                                       >  6:234620/1‑148 (MQ=255)
           ccAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGc                                                                                                                                    <  2:88821/148‑1 (MQ=255)
           ccAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGc                                                                                                                                    <  1:289206/148‑1 (MQ=255)
            cagcagGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCt                                                                                                                                   >  2:340060/1‑148 (MQ=255)
              gcagGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCTGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTg                                                                                                                                  >  4:92810/1‑147 (MQ=255)
                agGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGCTGTaa                                                                                                                               <  1:8312/148‑1 (MQ=255)
                     cTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATc                                                                                                                          >  3:286190/1‑148 (MQ=255)
                      tGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATc                                                                                                                          >  1:265548/1‑147 (MQ=255)
                        gTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTa                                                                                                                        <  5:186480/147‑1 (MQ=255)
                              aaCGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTaa                                                                                                                               >  2:30185/1‑134 (MQ=255)
                              aaCGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTaa                                                                                                                               <  1:30185/134‑1 (MQ=255)
                              aaCGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTAttt                                                                                                                 >  7:315896/1‑148 (MQ=255)
                              aaCGGCTTGCGTCACAAGATCCGTTAGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTAttt                                                                                                                 >  5:180616/1‑148 (MQ=255)
                               aCGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTAtt                                                                                                                  <  7:83322/146‑1 (MQ=255)
                               aCGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTAtt                                                                                                                  <  5:113947/146‑1 (MQ=255)
                                   cTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCAt                                                                                                            <  4:286190/148‑1 (MQ=255)
                                     tGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCt                                                                                                          >  7:35287/1‑148 (MQ=255)
                                       cGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGa                                                                                                        <  7:116479/148‑1 (MQ=255)
                                            caAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGACCGATGGTGGCGCTCGGCTGTAAATATCTa                                                                                                                        >  6:316151/1‑127 (MQ=255)
                                            caAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGACCGATGGTGGCGCTCGGCTGTAAATATCTa                                                                                                                        <  5:316151/127‑1 (MQ=255)
                                                     tttGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACg                                                                                          <  5:234620/148‑1 (MQ=255)
                                                      ttGCAGGTCGATGGCGGCCTGAAAAAGGGTGTCGATATCAGCACGGCGGCGTTTCTCGGCGCATAAAGCTTCGGCTTCGGCACTGGCCAGATGGTGGCGCGCGGCGGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACgg                                                                                         <  5:180722/148‑1 (MQ=255)
                                                       tGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACgg                                                                                         <  2:287564/147‑1 (MQ=255)
                                                          aGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGgtgt                                                                                     <  3:126370/148‑1 (MQ=255)
                                                            gTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAg                                                                                   >  2:266163/1‑148 (MQ=255)
                                                                 tGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGcc                                                                                                              <  5:171681/116‑1 (MQ=255)
                                                                 tGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGcc                                                                                                              >  6:171681/1‑116 (MQ=255)
                                                                 tGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAAc                                                                               >  7:36777/1‑147 (MQ=255)
                                                                       ccTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGAt                                                                        >  4:345260/1‑148 (MQ=255)
                                                                        cTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGGGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATg                                                                       <  5:286292/148‑1 (MQ=255)
                                                                            aaaCGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGc                                                                                  >  8:250031/1‑133 (MQ=255)
                                                                            aaaCGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGc                                                                                  <  7:250031/133‑1 (MQ=255)
                                                                                     tCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGt                                                            >  2:37536/1‑146 (MQ=255)
                                                                                      cGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGt                                                            >  1:184509/1‑145 (MQ=255)
                                                                                        atatCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGaa                                                       >  3:81212/1‑148 (MQ=255)
                                                                                        atatCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGaa                                                       <  7:241704/148‑1 (MQ=255)
                                                                                          atcatcAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAAc                                                      >  3:258620/1‑147 (MQ=255)
                                                                                           tcatcaAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCa                                                    >  6:317968/1‑148 (MQ=255)
                                                                                            catcaAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCAc                                                   >  2:244360/1‑148 (MQ=255)
                                                                                              tcaAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACt                                                  >  4:296642/1‑147 (MQ=255)
                                                                                               caAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTAt                                                >  1:164316/1‑148 (MQ=255)
                                                                                                acggcggcGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATc                                               <  3:228884/146‑1 (MQ=255)
                                                                                                aaGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCGGGATGACAAACTGCGTAAGAACCACTAt                                                <  1:126148/147‑1 (MQ=255)
                                                                                                 aGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCa                                              >  2:287100/1‑148 (MQ=255)
                                                                                                  ggcggcGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGCCAAACTGCGTAAGAACCACTATCAc                                             >  8:285709/1‑148 (MQ=255)
                                                                                                   gcggcgATTCTCGGCGCAGAAAGCTTCGGCTTCAGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACg                                            >  1:145771/1‑148 (MQ=255)
                                                                                                    cggcgATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACgg                                           <  1:340060/148‑1 (MQ=255)
                                                                                                       cgATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCt                                        >  3:307777/1‑148 (MQ=255)
                                                                                                            ctcGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCa                                    <  3:296642/147‑1 (MQ=255)
                                                                                                              cGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAATTATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCAtt                                  >  7:14515/1‑147 (MQ=255)
                                                                                                                  gcAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAATTATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAgg                             <  8:14515/148‑1 (MQ=255)
                                                                                                                  gcAGAAAGCTGCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAgg                             <  2:374564/148‑1 (MQ=255)
                                                                                                                     gAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACTACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGa                          >  5:255521/1‑148 (MQ=255)
                                                                                                                     gAAAGCTTCGGCTTCGGCACTGGCTCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGa                          <  1:266163/148‑1 (MQ=255)
                                                                                                                     gAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGAGGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGa                          <  6:180616/148‑1 (MQ=255)
                                                                                                                                   cggcACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAg            >  7:212327/1‑148 (MQ=255)
                                                                                                                                     gcACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGtt          <  8:36777/148‑1 (MQ=255)
                                                                                                                                      cACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGttt         <  8:315896/148‑1 (MQ=255)
                                                                                                                                       aCTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTa        <  1:287100/148‑1 (MQ=255)
                                                                                                                                        cTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTAt       <  2:145771/148‑1 (MQ=255)
                                                                                                                                             ccGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGccc  <  4:307777/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CTTGTTGAAACCCAGCAGGCGCTGGTTGCTAACGGCTTGCGTCACAAGATCCGTTTGCAGGTCGATGGCGGCCTGAAAACGGGTGTCGATATCATCAAGGCGGCGATTCTCGGCGCAGAAAGCTTCGGCTTCGGCACTGGCCCGATGGTGGCGCTCGGCTGTAAATATCTACGTATTTGCCATCTGAACAACTGCGCAACGGGTGTAGCAACTCAGGATGACAAACTGCGTAAGAACCACTATCACGGCCTGCCATTCAAGGTGACGAATTACTTTGAGTTTATCGCCC  >  NC_000913/3357869‑3358157

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