Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 194,272 A→G 5.8% H251R (CAC→CGC)  dxr → 1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913194,2720AG 5.8% 212.6 / 3.3 70H251R (CAC→CGC) dxr1‑deoxy‑D‑xylulose 5‑phosphate reductoisomerase
Reads supporting (aligned to +/- strand):  ref base A (40/25);  new base G (2/2);  total (43/27)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCGTGAAGCCGCTCGATTTTTGC  >  NC_000913/194132‑194420
                                                                                                                                            |                                                                                                                                                    
tCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGcagt                                                                                                                                              >  8:234334/1‑149 (MQ=255)
    gCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGcagtcagt                                                                                                                                          >  4:74728/1‑149 (MQ=255)
        gCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCGCCCGCAGTCAGTGAtt                                                                                                                                      >  4:374922/1‑149 (MQ=255)
        gCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCGCCCGCAGTCAGTGAtt                                                                                                                                      >  4:369055/1‑149 (MQ=255)
          agTCATCCGAATTGGTCGATGTGTCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCTCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCa                                                                                                                                    <  6:392341/148‑1 (MQ=255)
           gTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCAc                                                                                                                                   >  2:91649/1‑149 (MQ=255)
           gTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGGCTGCAATACATAGAAGCGCGTTGCCGGTTTAACGCCACCGACAGCCAGATGGAAGTGCAGATTCACCCGCAGTCAGTGATTCa                                                                                                                                    >  3:279769/1‑148 (MQ=255)
            tCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACt                                                                                                                                  <  7:328190/149‑1 (MQ=255)
             cATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACt                                                                                                                                  >  3:69786/1‑148 (MQ=255)
              aTCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTc                                                                                                                                 >  8:3553/1‑148 (MQ=255)
               tCCGCACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCaa                                                                                                                               <  1:174756/149‑1 (MQ=255)
               tCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGAtt                                                                                                                                      >  5:170338/1‑142 (MQ=255)
               tCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGAtt                                                                                                                                      <  6:170338/142‑1 (MQ=255)
               tCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCaa                                                                                                                               <  3:74728/149‑1 (MQ=255)
                  gAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATgg                                                                                                                            >  8:186426/1‑149 (MQ=255)
                    aCTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTg                                                                                                                          <  3:286834/149‑1 (MQ=255)
                       ggTCGATTGGGCGAAAAATTTCTGTCGATTCGGCTACAATGATGAACAAAGGTATGGAATAAATTGAAGATAATTCGCTGTTTAACGCAAGAGCCAGCCACATGGAATGGCATAATCAACCGCAGAAAGAGATTAACTCAATGGTtcg                                                                                                                        >  8:192900/1‑148 (MQ=255)
                        gTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCCGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTgcgc                                                                                                                       >  6:117965/1‑148 (MQ=255)
                         tCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTa                                                                                                                     <  7:186426/149‑1 (MQ=255)
                           gATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTAt                                                                                                                    >  7:65227/1‑148 (MQ=255)
                            aTGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCa                                                                                                                  <  1:91649/149‑1 (MQ=255)
                                     aaaaTTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTTAATGGTGCGCTATCAGGACGGCAg                                                                                                         >  4:109558/1‑149 (MQ=255)
                                               tCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGgcg                                                                                               <  8:35058/149‑1 (MQ=255)
                                               tCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTAGCAGGACGGCAGTGTTCTGgcg                                                                                               >  2:288374/1‑149 (MQ=255)
                                               tCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTCTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGgcg                                                                                               >  8:146222/1‑149 (MQ=255)
                                                cGATTCGGCGACCATGATGAACAAAGGTCTGGTAGACATTGAAGCGGGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGTAGTGCTGAGTCACCCGCAGGCAGTGTTTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGgcgc                                                                                              <  8:316262/149‑1 (MQ=255)
                                                 gATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCa                                                                                             <  7:234334/149‑1 (MQ=255)
                                                  aTTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGgcg                                                                                               >  7:290689/1‑146 (MQ=255)
                                                            cATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAAc                                                                                  >  5:314054/1‑149 (MQ=255)
                                                            cATGATGAACAAAGGGCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAAc                                                                                  <  7:239083/149‑1 (MQ=255)
                                                                 tgaACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATc                                                                                                                   <  5:189215/111‑1 (MQ=255)
                                                                 tgaACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATc                                                                                                                   >  6:189215/1‑111 (MQ=255)
                                                                   aaCAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGatat                                                                           >  5:202746/1‑149 (MQ=255)
                                                                     cAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGc                                                                         >  8:224320/1‑149 (MQ=255)
                                                                      aaaGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCg                                                                        >  7:44044/1‑149 (MQ=255)
                                                                        aGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTa                                                                      <  6:202746/149‑1 (MQ=255)
                                                                           tcgggACTACATTGAAGCGCGGTGGCTGGTTAACGCCTGCGCCAGCCAGTTGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCTCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGc                                                                   <  4:279769/146‑1 (MQ=255)
                                                                             tGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGcc                                                                  >  4:134038/1‑148 (MQ=255)
                                                                             tGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGcc                                                                  >  2:325816/1‑148 (MQ=255)
                                                                              ggAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCaa                                                                >  2:68948/1‑149 (MQ=255)
                                                                               gAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCaa                                                                <  1:325816/148‑1 (MQ=255)
                                                                                    cATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGcc                                                           <  6:314054/148‑1 (MQ=255)
                                                                                     aTTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCca                                                         >  6:287465/1‑149 (MQ=255)
                                                                                      ttGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAAtt                                                              <  7:151841/143‑1 (MQ=255)
                                                                                      ttGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAAtt                                                              >  8:151841/1‑143 (MQ=255)
                                                                                      ttGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGcc                                                           >  1:219265/1‑146 (MQ=255)
                                                                                      ttGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGcc                                                           <  2:219265/146‑1 (MQ=255)
                                                                                         aaGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCacac                                                      >  8:267724/1‑148 (MQ=255)
                                                                                           gcgcGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACcat                                                   >  5:329528/1‑149 (MQ=255)
                                                                                           gcgcGGTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTTAATGGTGCGCTATCAGGACGGCAGTGGTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACcat                                                   <  3:109558/149‑1 (MQ=255)
                                                                                                     tGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGcc                                          >  2:365936/1‑148 (MQ=255)
                                                                                                       tttAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGaa                                       >  2:300473/1‑149 (MQ=255)
                                                                                                          aaCGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATc                                     <  4:69786/148‑1 (MQ=255)
                                                                                                             gccagcgccagcCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTgggg                                                                                      >  7:362518/1‑96 (MQ=255)
                                                                                                             gccagcgccagcCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTgggg                                                                                      <  8:362518/96‑1 (MQ=255)
                                                                                                              ccagcgccagcCAGATGGAAGTGCTGATTCGCCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTg                                <  3:369055/149‑1 (MQ=255)
                                                                                                              ccagcgccagcCAGATGGAAGTGCTGATTCGCCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTg                                <  3:374922/149‑1 (MQ=255)
                                                                                                                   gccagccaGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAAt                                      >  4:161938/1‑138 (MQ=255)
                                                                                                                   gccagccaGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAAt                                      <  3:161938/138‑1 (MQ=255)
                                                                                                                    ccagccagATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAActct                          >  2:65865/1‑149 (MQ=255)
                                                                                                                        ccagATGGAAGTGCTGATTCCCCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCg                      >  4:63562/1‑149 (MQ=255)
                                                                                                                        ccagATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCg                      >  5:360546/1‑149 (MQ=255)
                                                                                                                         cagATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACAGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCGt                     >  8:412212/1‑149 (MQ=255)
                                                                                                                               gAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCGTGAAGc                <  7:146222/148‑1 (MQ=255)
                                                                                                                                     cTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGa                                        <  5:216905/118‑1 (MQ=255)
                                                                                                                                     cTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGa                                        >  6:216905/1‑118 (MQ=255)
                                                                                                                                     cTGATTAACCCGCAGTCAGTGATTCACTAAATGGTGCGATATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACTCCAATTGCACACACCATGTCATGTACGAATCGCGTGAACTCTTGCGTTAAGACTCTCg          >  6:22498/1‑148 (MQ=255)
                                                                                                                                     atGATTCACCCGCAGTCAGTGAATCACTCAATGGTGCGCTATCAGGACGGAAGTCTTCTAGAGAAGCTGGGGGAACCGGATATGCGTACGCCAATTCCCCACACCATCGCATGGCAGAATCGCGTGAACTCTGGCGTGAAGCCTCTCGa         >  4:98175/2‑149 (MQ=255)
                                                                                                                                       gATTCACCCGCAGTCAGTGATTCCCTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCGTGAAGCCGCTCGAtt       >  5:243829/1‑149 (MQ=255)
                                                                                                                                        aTTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCGTGAAGCCGCTCGAttt      >  2:199379/1‑149 (MQ=255)
                                                                                                                                            aCCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTg                                                             <  4:179607/90‑1 (MQ=255)
                                                                                                                                            aCCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTg                                                             >  3:179607/1‑90 (MQ=255)
                                                                                                                                            aCCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGc                       >  6:156058/1‑128 (MQ=255)
                                                                                                                                            aCCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGc                       <  5:156058/128‑1 (MQ=255)
                                                                                                                                            aCCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCGTGAAGCCGCTCGATTTTTg   <  8:65227/148‑1 (MQ=255)
                                                                                                                                            aCCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCGTGAAGCCGCTCGATTTTTGc  <  8:44044/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
TCAAGCCTGCCGTCATCCGAACTGGTCGATGGGGCGTAAAATTTCTGTCGATTCGGCTACCATGATGAACAAAGGTCTGGAATACATTGAAGCGCGTTGGCTGTTTAACGCCAGCGCCAGCCAGATGGAAGTGCTGATTCACCCGCAGTCAGTGATTCACTCAATGGTGCGCTATCAGGACGGCAGTGTTCTGGCGCAGCTGGGGGAACCGGATATGCGTACGCCAATTGCCCACACCATGGCATGGCCGAATCGCGTGAACTCTGGCGTGAAGCCGCTCGATTTTTGC  >  NC_000913/194132‑194420

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: ‑