Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 911,159 G→T 5.2% intergenic (‑110/+23) poxB ← / ← hcr pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding/HCP oxidoreductase, NADH‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913911,1590GT 5.2% 267.0 / 2.3 77intergenic (‑110/+23)poxB/hcrpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding/HCP oxidoreductase, NADH‑dependent
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base T (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCGGAAACCACTTTTGTC  >  NC_000913/911011‑911305
                                                                                                                                                    |                                                                                                                                                  
gAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAgg                                                                                                                                                    <  7:71431/149‑1 (MQ=255)
          gCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATAtcatcg                                                                                                                                          >  7:213279‑M2/1‑143 (MQ=255)
           cGATATAAGCTGCACCCGTTTGTTTCATTGTTCTCCATCTCCTGAATGTGTTAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGTGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATAtcatcgg                                                                                                                                         >  7:395237‑M2/1‑142 (MQ=255)
            gATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATAtcatcggg                                                                                                                                        <  1:136769‑M2/149‑9 (MQ=255)
                   gCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAgg                                                                                                                                 >  6:205622/1‑149 (MQ=255)
                   gCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTGCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAgg                                                                                                                                 >  4:16951/1‑149 (MQ=255)
                     tGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGc                                                                                                                                >  2:315449/1‑148 (MQ=255)
                         aCCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGaa                                                                                                                                                  <  4:280530/126‑1 (MQ=255)
                         aCCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGaa                                                                                                                                                  >  3:280530/1‑126 (MQ=255)
                            gtttgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGc                                                                                                                        <  6:381429/149‑1 (MQ=255)
                               tgtttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCga                                                                                                                      <  1:200288/148‑1 (MQ=255)
                                 tttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATg                                                                                                                         >  6:40752/1‑143 (MQ=255)
                                 tttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATg                                                                                                                         <  5:40752/143‑1 (MQ=255)
                                 tttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGaa                                                                                                                   >  4:344725/1‑149 (MQ=255)
                                 tttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGaa                                                                                                                   >  4:256186/1‑149 (MQ=255)
                                 tttCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGaa                                                                                                                   >  8:32171/1‑149 (MQ=255)
                                          tctcCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGATATTCCTTtcatcgggctatttaaccgttagtgcctcctttctctc                                                                                                          >  4:93412‑M2/1‑111 (MQ=255)
                                              cATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATccccc                                                                                                        <  2:157737/147‑1 (MQ=255)
                                                  tcCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGTTTAAATAGCCCGATGAAAGGAATATCATCGGTCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGa                                                                                                  <  5:105607/149‑1 (MQ=255)
                                                  tcCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGATATTCCTTtcatcgggctatttaaccgttagtgcctccttt                                                                                                               <  6:221589‑M2/136‑34 (MQ=255)
                                                  tcCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGATATTCCTTtcatcgggctatttaaccgttagtgcctccttt                                                                                                               >  5:221589‑M2/1‑103 (MQ=255)
                                                   ccTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAATGGATGGGAGAGAAATGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTTCGGAt                                                                                                 >  5:317157/1‑149 (MQ=255)
                                                    cTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATg                                                                                                <  8:134555/149‑1 (MQ=255)
                                                    cTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATg                                                                                                <  8:128766/149‑1 (MQ=255)
                                                       aaTGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATgg                                                                                               >  2:327365/1‑147 (MQ=255)
                                                                  ggTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGc                                                                                    >  3:15194/1‑147 (MQ=255)
                                                                  ggTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGc                                                                                    <  4:15194/147‑1 (MQ=255)
                                                                  ggTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCa                                                                                  <  7:10341/149‑1 (MQ=255)
                                                                   gTAACAAGTTTAGTTCATATGACTTTGGGGGAAGGGTTGGGAGAGAAATGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAATGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGTAGCAGGCCAg                                                                                 <  7:356692/149‑1 (MQ=255)
                                                                     aaCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTa                                                                               >  2:118579/1‑149 (MQ=255)
                                                                      aCAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTAc                                                                              <  3:344725/149‑1 (MQ=255)
                                                                         aGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGACAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTa                                                                           <  3:16951/149‑1 (MQ=255)
                                                                           tttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCga                                                                                                                      >  7:197237/1‑104 (MQ=255)
                                                                           tttAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCga                                                                                                                      <  8:197237/104‑1 (MQ=255)
                                                                               gTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGTAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCtt                                                                      <  6:300238/148‑1 (MQ=255)
                                                                                ttCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCtt                                                                      >  6:364993/1‑147 (MQ=255)
                                                                                   aTCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGc                                                                  >  3:180616/1‑148 (MQ=255)
                                                                                    tCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGTAGGAGCAGGCCAGTACGTAACCTTCAGCGa                                                                >  8:296397/1‑149 (MQ=255)
                                                                                      tGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGAtt                                                              >  6:170019/1‑149 (MQ=255)
                                                                                       gACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGAttt                                                             >  7:180700/1‑149 (MQ=255)
                                                                                         cGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCg                                                           <  7:32171/149‑1 (MQ=255)
                                                                                          ggAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCgg                                                          >  8:122170/1‑149 (MQ=255)
                                                                                           gAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGc                                                         <  1:315449/149‑1 (MQ=255)
                                                                                              ggggAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGTGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTc                                                      >  6:367193/1‑149 (MQ=255)
                                                                                               gggAAGTGATGGGAGAGAAAGGAGGCACTTACGGTTAAATAGCCCGATGAAAGGAATATCATCGGTCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGATCTGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCg                                                     <  4:81574/149‑1 (MQ=255)
                                                                                                 gAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGAGAGAACCAAATCCCCCTGCGGATGGCAGGAGAAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGt                                                   >  1:397613/1‑149 (MQ=255)
                                                                                                      gATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAACCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCg                                              >  8:312848/1‑149 (MQ=255)
                                                                                                           gagagaAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAACCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCAtt                                         <  7:312848/149‑1 (MQ=255)
                                                                                                                aaaGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGtgct                                    >  7:1301/1‑149 (MQ=255)
                                                                                                                     aggCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCa                               <  5:364993/149‑1 (MQ=255)
                                                                                                                      ggCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCAc                              >  7:174388/1‑149 (MQ=255)
                                                                                                                         aCTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGt                           <  5:170019/149‑1 (MQ=255)
                                                                                                                          cTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGta                          <  3:256186/149‑1 (MQ=255)
                                                                                                                           tAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGtgct                                    >  2:393314/1‑138 (MQ=255)
                                                                                                                           tAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGtgct                                    <  1:393314/138‑1 (MQ=255)
                                                                                                                             aCGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATgg                                                                                               >  1:153609/1‑77 (MQ=255)
                                                                                                                             aCGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATgg                                                                                               <  2:153609/77‑1 (MQ=255)
                                                                                                                             aCGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGc                                                                                    <  6:32958/88‑1 (MQ=255)
                                                                                                                             aCGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGc                                                                                    >  5:32958/1‑88 (MQ=255)
                                                                                                                             aCGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCg                                                           >  1:70896/1‑113 (MQ=255)
                                                                                                                             aCGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCg                                                           <  2:70896/113‑1 (MQ=255)
                                                                                                                             aCGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTatat                       >  5:143229/1‑149 (MQ=255)
                                                                                                                               ggTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATAtt                      >  5:183981/1‑148 (MQ=255)
                                                                                                                                 ttAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCAc                   <  7:122170/149‑1 (MQ=255)
                                                                                                                                 ttAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCAc                   >  1:305129/1‑149 (MQ=255)
                                                                                                                                 ttAAATAGCACGATGAAAGGAATCTCATCGGGCATAAGGCGATTATGCGAGAACCAAATACCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCaa                   <  4:194536/149‑2 (MQ=255)
                                                                                                                                    aaTAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCg                 <  4:180616/148‑1 (MQ=255)
                                                                                                                                     aTAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCGGa               <  1:198625/149‑1 (MQ=255)
                                                                                                                                         cccGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCAtt                                         >  3:174097/1‑119 (MQ=255)
                                                                                                                                         cccGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCAtt                                         <  4:174097/119‑1 (MQ=255)
                                                                                                                                           cGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCGGAAACCAc         >  1:332831/1‑149 (MQ=255)
                                                                                                                                            gATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCGGAAACCACt        >  3:119736/1‑149 (MQ=255)
                                                                                                                                             aTGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCGGAAACCACtt       <  2:305129/149‑1 (MQ=255)
                                                                                                                                               gAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGTGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCGGAAACCACtttt     >  1:157194/1‑149 (MQ=255)
                                                                                                                                                 aaGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGc                                                                                    >  3:22129/1‑68 (MQ=255)
                                                                                                                                                 aaGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGc                                                                                    <  4:22129/68‑1 (MQ=255)
                                                                                                                                                  aGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGTAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCGGAAACCACTTTTGTc  <  7:296397/149‑1 (MQ=255)
                                                                                                                                                  aGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCGGAAACCACTTTTGTc  >  1:293515/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GAGTGTTTTGGCGATATAAGCTGCAACCGTTTGTTTCATGGTTCTCCATCTCCTGAATGTGATAACGGTAACAAGTTTAGTTCATCTGACGGAGGGGGAAGGGATGGGAGAGAAAGGAGGCACTAACGGTTAAATAGCCCGATGAAAGGAATATCATCGGGCATAAGGCGATTATGCGAGAACCAAATCCCCCTGCGGATGGCAGGAGCAGGCCAGTACGTAACCTTCAGCGATTTCGGCGTCGGTCAGCGTCATTGTGCTGCTCACCGTATATTCACCGGAAACCACTTTTGTC  >  NC_000913/911011‑911305

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