Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,640,668 G→T 7.0% intergenic (‑93/+18) hisS ← / ← ispG histidyl tRNA synthetase/1‑hydroxy‑2‑methyl‑2‑(E)‑butenyl 4‑diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,640,6680GT 7.0% 395.6 / 14.5 118intergenic (‑93/+18)hisS/ispGhistidyl tRNA synthetase/1‑hydroxy‑2‑methyl‑2‑(E)‑butenyl 4‑diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (48/61);  new base T (4/4);  total (53/65)
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.72e-01

Only 100 of 125 total aligned reads displayed.

TTTCGCCAGGCAGGTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACGCCATCTTCAT  >  NC_000913/2640521‑2640816
                                                                                                                                                   |                                                                                                                                                    
tttCGCCAGGCAGGTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAgg                                                                                                                                                      >  3:879111/1‑148 (MQ=255)
tttCGCCAGGCAGGTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGgc                                                                                                                                                     <  2:988520/149‑1 (MQ=255)
tttCGCCAGGCAGGTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGgc                                                                                                                                                     >  2:142737/1‑149 (MQ=255)
    gCCAGGCAGGTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCt                                                                                                                                                  <  4:959727/148‑1 (MQ=255)
         gcaggTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCAt                                                                                                                                            <  1:125417/149‑1 (MQ=255)
         gcaggTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCAt                                                                                                                                            <  2:463934/149‑1 (MQ=255)
            ggTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCAc                                                                                                                                         >  2:1043684/1‑149 (MQ=255)
             gTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCAc                                                                                                                                         >  4:764736/1‑148 (MQ=255)
              tAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACg                                                                                                                                        >  3:875616/1‑148 (MQ=255)
                aTCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGtt                                                                                                                                      <  4:568815/148‑1 (MQ=255)
                 tCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGCCTCGCTTcccgt                                                                                                                                            <  1:371130/141‑3 (MQ=255)
                 tCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGCCTCGCTTcccgt                                                                                                                                            >  2:371130/1‑139 (MQ=255)
                   gTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGttatta                                                                                                                                  >  1:901974/1‑149 (MQ=255)
                    ttCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGttatta                                                                                                                                  >  3:648517/1‑148 (MQ=255)
                        tGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTAtttt                                                                                                                              >  3:606415/1‑148 (MQ=255)
                          ccGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTc                                                                                                                            >  3:924989/1‑148 (MQ=255)
                           cgcgAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCaa                                                                                                                          >  2:860171/1‑149 (MQ=255)
                           cgcgAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCa                                                                                                                           >  3:1035330/1‑148 (MQ=255)
                                aTGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAATCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACctg                                                                                                                      <  4:648517/148‑1 (MQ=255)
                                aTGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACctgc                                                                                                                     <  2:137510/149‑1 (MQ=255)
                                    cTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGaa                                                                                                                 >  1:89895/1‑149 (MQ=255)
                                     ttGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGaa                                                                                                                 <  4:864433/148‑1 (MQ=255)
                                     ttGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGaa                                                                                                                 >  4:716703/1‑148 (MQ=255)
                                     ttGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAAc                                                                                                                <  2:85110/149‑1 (MQ=255)
                                      tGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCaa                                                                                                                          <  4:941700/138‑1 (MQ=255)
                                       gAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACTCGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGt                                                                                                              <  2:126278/149‑1 (MQ=255)
                                       gAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGt                                                                                                              <  1:764823/149‑1 (MQ=255)
                                         aTGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTc                                                                                                             <  4:924989/148‑1 (MQ=255)
                                            ttttttGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGCAGCCATGCGGGTTCAATCATACACGGGAAGCGATGCTCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAAt                                                                                                          <  4:701037/148‑1 (MQ=255)
                                               tttGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCg                                                                                                       <  4:428346/148‑1 (MQ=255)
                                               tttGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCg                                                                                                       >  4:830301/1‑148 (MQ=255)
                                                  gCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGAcg                                                                                                    <  3:207348/148‑1 (MQ=255)
                                                    cACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACgcgc                                                                                                 >  1:714325/1‑149 (MQ=255)
                                                      cGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTc                                                                                                                                                >  3:990134/1‑100 (MQ=255)
                                                        tttATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTc                                                                                              >  3:195220/1‑148 (MQ=255)
                                                        tttATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCg                                                                                             >  1:368408/1‑149 (MQ=255)
                                                         ttATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCg                                                                                             >  3:379889/1‑148 (MQ=255)
                                                         ttATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCg                                                                                             >  4:151788/1‑148 (MQ=255)
                                                            ttCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCa                                                                                         >  1:220410/1‑149 (MQ=255)
                                                               tctTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCGTCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCt                                                                                      >  1:611742/1‑149 (MQ=255)
                                                                 tttCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGa                                                                                    <  2:529861/149‑1 (MQ=255)
                                                                   tCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACt                                                                                  >  1:289138/1‑149 (MQ=255)
                                                                    cTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGAGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACt                                                                                  <  3:924143/148‑1 (MQ=255)
                                                                     tGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGCCTCGCTTcccgt                                                                                                                                            <  2:707678/89‑3 (MQ=255)
                                                                      gAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGc                                                                               >  2:52039/1‑149 (MQ=255)
                                                                          ataAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTc                                                                           <  1:860171/149‑1 (MQ=255)
                                                                            aaaaaTGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTcc                                                                                                                                               <  4:303022/79‑1 (MQ=255)
                                                                            aaaaaTGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTcc                                                                                                                                               >  3:303022/1‑79 (MQ=255)
                                                                                    aCCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGc                                                                 <  2:368408/149‑1 (MQ=255)
                                                                                     cccTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGc                                                                 >  3:308040/1‑148 (MQ=255)
                                                                                       cTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTg                                                              >  1:805085/1‑149 (MQ=255)
                                                                                        tCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTg                                                              >  4:776333/1‑148 (MQ=255)
                                                                                        tCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGc                                                             >  2:624538/1‑149 (MQ=255)
                                                                                         cAACGCTTCCCTCAATGTTTCGGGAGCCATGCGTGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGTTTCGTCCAGCTGACTGGCTTTCGCACGAACGCGTGc                                                             >  3:705545/1‑148 (MQ=255)
                                                                                         cAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGc                                                             >  3:332929/1‑148 (MQ=255)
                                                                                            cGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCCGCTGACTGGCTTTCGCACGAATGCGTGCTTc                                                          <  4:904914/148‑1 (MQ=255)
                                                                                            cGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGCATTCGTGcgaaa                                                         <  1:465331/149‑5 (MQ=255)
                                                                                            cGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGCATTCGTGcgaaa                                                         <  2:1014990/149‑5 (MQ=255)
                                                                                              cTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCa                                                        <  4:620819/148‑1 (MQ=255)
                                                                                               ttCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGc                                                      >  2:1070454/1‑149 (MQ=255)
                                                                                                tCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCt                                                     <  1:142737/149‑1 (MQ=255)
                                                                                                  ccTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTg                                                    <  4:606415/148‑1 (MQ=255)
                                                                                                   cTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTgg                                                   <  3:658546/148‑1 (MQ=255)
                                                                                                   cTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTgg                                                   >  4:658546/1‑148 (MQ=255)
                                                                                                   cTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGt                                                  <  2:685071/149‑1 (MQ=255)
                                                                                                    tCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTc                                                 >  2:92745/1‑149 (MQ=255)
                                                                                                     cAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCttt                                                                            <  3:46097/121‑1 (MQ=255)
                                                                                                     cAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCttt                                                                            >  4:46097/1‑121 (MQ=255)
                                                                                                        tGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGAt                                              >  4:473268/1‑148 (MQ=255)
                                                                                                        tGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATc                                             <  1:1043684/149‑1 (MQ=255)
                                                                                                        tGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATc                                             <  1:92745/149‑1 (MQ=255)
                                                                                                        tGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATc                                             <  2:714325/149‑1 (MQ=255)
                                                                                                        tGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATc                                             <  2:797937/149‑1 (MQ=255)
                                                                                                           ttCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCata                                          <  2:998048/149‑1 (MQ=255)
                                                                                                               ggAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTg                                                                                 <  1:406575/106‑1 (MQ=255)
                                                                                                                   ccATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgttgtt                                  <  1:632113/149‑1 (MQ=255)
                                                                                                                   ccATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgttgtt                                  >  2:242961/1‑149 (MQ=255)
                                                                                                                   ccATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgttgt                                   <  3:335654/148‑1 (MQ=255)
                                                                                                                      tGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTc                               <  1:297688/149‑1 (MQ=255)
                                                                                                                       gCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTc                               <  3:764736/148‑1 (MQ=255)
                                                                                                                       gCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTc                               <  3:606219/148‑1 (MQ=255)
                                                                                                                          ggTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAg                            >  4:138557/1‑148 (MQ=255)
                                                                                                                           gTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGa                           <  4:308040/148‑1 (MQ=255)
                                                                                                                            ttCAATCATACACGGGAAGCGAGAGGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACACGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGAc                          >  4:295724/1‑148 (MQ=255)
                                                                                                                                aTCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgtt                                     <  2:13557/133‑1 (MQ=255)
                                                                                                                                aTCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgtt                                     >  1:13557/1‑133 (MQ=255)
                                                                                                                                 tCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGa                                                                                                      <  2:954885/67‑1 (MQ=255)
                                                                                                                                   aTACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgttgtt                                  <  1:97406/133‑1 (MQ=255)
                                                                                                                                   aTACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgttgtt                                  >  2:97406/1‑133 (MQ=255)
                                                                                                                                    tACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgt                                      <  4:152337/128‑1 (MQ=255)
                                                                                                                                       atgggAAGCGCCTCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgtt                                     >  1:537877/3‑126 (MQ=255)
                                                                                                                                          ggAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCgt                                      >  3:222386/1‑122 (MQ=255)
                                                                                                                                            aaGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACGCCa         <  1:809415/149‑1 (MQ=255)
                                                                                                                                              gCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTTCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACGCCAt        <  4:319410/148‑1 (MQ=255)
                                                                                                                                              gCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACGCCAt        <  4:742206/148‑1 (MQ=255)
                                                                                                                                               cGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACGCCATc       <  4:875616/148‑1 (MQ=255)
                                                                                                                                                  ggCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCACTTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACGCCATCTTc    >  3:605435/1‑148 (MQ=255)
                                                                                                                                                   gcgcTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACGCCATCTTCAt  >  1:266762/1‑149 (MQ=255)
                                                                                                                                                   gcgcTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACGCCATCTTCAt  <  1:624386/149‑1 (MQ=255)
                                                                                                                                                   gcgcTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACACCATCTTCAt  <  2:226666/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTCGCCAGGCAGGTAATCGTTCATGCCGCGAATGGCTTGAATGTTTTTTGCCACGTTTATTCTCTTTCTGAATATAAAAATGAACCCTCAACGCTTCCCTCAATGTTTCGGGAGCCATGCGGGTTCAATCATACACGGGAAGCGAGGCGCTTCCCATCACGTTATTATTTTTCAACCTGCTGAACGTCAATTCGACGCGCTTCGTCCAGCTGACTGGCTTTCGCACGAATGCGTGCTTCCAGCTGGTCGATCATATCGTTGTTGTCCAGACGGTCTTTGCGCACGCCATCTTCAT  >  NC_000913/2640521‑2640816

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: