Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 441,069 A→T 9.7% V11V (GTT→GTA) ‡ ispA ← geranyl diphosphate/farnesyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913441,0690AT 9.7% 172.9 / 6.8 61V11V (GTT→GTA) ‡ispAgeranyl diphosphate/farnesyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base T (3/3);  total (29/32)
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.44e-01

GCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTG  >  NC_000913/440922‑441214
                                                                                                                                                   |                                                                                                                                                 
gCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAAc                                                                                                                                                 <  2:462062/150‑1 (MQ=255)
        tatGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCtt                                                                                                                                         >  1:838779/1‑150 (MQ=255)
           gACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGGTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCGTTAAGCAGGccatcc                                                                                                                                      <  1:1026745/150‑6 (MQ=255)
            aCCGGTGGCATAAAGCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAg                                                                                                                                     <  2:1312925/150‑1 (MQ=255)
            aCCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCGTTAAGCAGGccaacca                                                                                                                                     <  1:750815/150‑7 (MQ=255)
            aCCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGGTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCGGCTTAACGCAGGCTTCGAg                                                                                                                                     <  2:685911/150‑1 (MQ=255)
               ggTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtg                                                                                                                                  >  1:7828/1‑150 (MQ=255)
               ggTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtg                                                                                                                                  >  1:166546/1‑150 (MQ=255)
                gTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTtgc                                                                                                                                 <  1:955773/150‑1 (MQ=255)
                    cATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGTGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCg                                                                                                                             >  1:356801/1‑150 (MQ=255)
                       aaaCCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGaa                                                                                                                          <  2:395631/150‑1 (MQ=255)
                           cAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAATTc                                                                                                                      <  1:348320/150‑1 (MQ=255)
                             ggAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCa                                                                                                                    <  1:251214/150‑1 (MQ=255)
                                  ggTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACt                                                                                                               >  1:1408252/1‑150 (MQ=255)
                                     cGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGTCCTGCTTAACGCAGGCTTCGAGATGCTGCGGAAAGTCCATTACTCAt                                                                                                            >  2:668964/1‑150 (MQ=255)
                                     cGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCAt                                                                                                            >  2:62384/1‑150 (MQ=255)
                                           cgcTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCg                                                                                                      <  2:154426/150‑1 (MQ=255)
                                             cTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGgt                                                                                                    <  2:7828/150‑1 (MQ=255)
                                                    ccTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAg                                                                                             <  2:838779/150‑1 (MQ=255)
                                                           atGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTag                                                                                      >  2:1021925/1‑150 (MQ=255)
                                                                  aTACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGTGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGATGCg                                                                               <  1:1310049/150‑1 (MQ=255)
                                                                  aTACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCg                                                                               <  1:536315/150‑1 (MQ=255)
                                                                   tACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGt                                                                              >  2:39314/1‑150 (MQ=255)
                                                                     atGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCTGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCt                                                                            <  1:590465/149‑1 (MQ=255)
                                                                      tGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGTGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTATAGAGGcttctt                                                                           <  2:356801/150‑1 (MQ=255)
                                                                        cATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCa                                                                         >  2:835402/1‑150 (MQ=255)
                                                                            gTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGt                                                                     <  1:1021925/150‑1 (MQ=255)
                                                                             tttCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTc                                                                    <  2:166546/150‑1 (MQ=255)
                                                                                 gACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGAc                                                                >  2:241959/1‑150 (MQ=255)
                                                                                   ccACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGAcag                                                              >  1:640304/1‑150 (MQ=255)
                                                                                    cACGGGAGTGTTCTGATAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGAcagc                                                             <  2:180642/150‑1 (MQ=255)
                                                                                          aGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAAtt                                                       <  2:33467/150‑1 (MQ=255)
                                                                                            tgtTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTg                                                     <  2:1140194/150‑1 (MQ=255)
                                                                                             gtTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGt                                                    >  1:34314/1‑150 (MQ=255)
                                                                                             gtTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGt                                                    >  1:128284/1‑150 (MQ=255)
                                                                                                cTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCATGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAcg                                                 >  2:1271691/1‑150 (MQ=255)
                                                                                                 tGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTAcgc                                                >  1:436797/1‑150 (MQ=255)
                                                                                                  gAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAAtt                                                       >  2:1214229/1‑142 (MQ=255)
                                                                                                  gAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAAtt                                                       <  1:1214229/142‑1 (MQ=255)
                                                                                                   aaaGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACgcgc                                              >  1:949291/1‑150 (MQ=255)
                                                                                                      gggCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGt                                           >  1:441333/1‑150 (MQ=255)
                                                                                                           gTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGc                                      >  1:1078576/1‑150 (MQ=255)
                                                                                                            tGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCt                                     <  2:1408252/150‑1 (MQ=255)
                                                                                                               ggCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTGATTGTCAGACAGCAGCATTTGTACGCGCTGTTCGGCttgt                                  >  1:1430554/1‑150 (MQ=255)
                                                                                                                gCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCttgtt                                 <  2:640304/150‑1 (MQ=255)
                                                                                                                 cGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCttgttg                                >  2:105297/1‑150 (MQ=255)
                                                                                                                      aaaCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAAtt                           >  2:395645/1‑150 (MQ=255)
                                                                                                                       aaCGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAAttt                          >  1:1178235/1‑150 (MQ=255)
                                                                                                                           gCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGcc                      <  1:105297/150‑1 (MQ=255)
                                                                                                                           gCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGcc                      <  2:128284/150‑1 (MQ=255)
                                                                                                                               aGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGcc                  <  1:835402/150‑1 (MQ=255)
                                                                                                                                 cgcCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGcccc                <  2:186103/150‑1 (MQ=255)
                                                                                                                                    cTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCTCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGa             <  2:1363106/150‑1 (MQ=255)
                                                                                                                                     tGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGAc            <  1:39314/150‑1 (MQ=255)
                                                                                                                                     ctcgaagCCTGCGTTAAGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGAc            >  1:424990/7‑150 (MQ=255)
                                                                                                                                     ctcgaagCCTGCGTTAAGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGAc            >  1:228149/7‑150 (MQ=255)
                                                                                                                                      ggTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACg           <  2:844201/150‑1 (MQ=255)
                                                                                                                                         tGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGc        >  2:371313/1‑150 (MQ=255)
                                                                                                                                         aagCCTGCGTTAAGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGc        >  1:1257926/3‑150 (MQ=255)
                                                                                                                                          agCCTGCGTTAAGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGcc       <  2:228149/149‑1 (MQ=255)
                                                                                                                                               gCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTg  >  1:149993/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCGAACATATGACCGGTGGCATAAACCAGGAAAGGTCGCAGGCGCTTACCACCTAATAATGCGCCATACTGCATGGTTTCGACCACGGGAGTGTTCTGAAAGGGCAGTGGGGCGATAAAACGGCTCAGCGCCTGGTTGGCCTGCTTAACGCAGGCTTCGAGTTGCTGCGGAAAGTCCATTACTCATTGTCCGGTGTAAAAGGGGTTAGAGAGGCGTCTTCATTGTCAGACAGCAGAATTTGTACGCGCTGTTCGGCTTGTTGTAATTTGGCCTGCCCCTGACGTGCCAGCTG  >  NC_000913/440922‑441214

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