Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,672,640 G→T 6.6% Q198H (CAG→CAT) ‡ hcaD → putative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,672,6400GT 6.6% 252.5 / 4.9 75Q198H (CAG→CAT) ‡hcaDputative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base G (27/43);  new base T (3/2);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

ATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGCG  >  NC_000913/2672491‑2672788
                                                                                                                                                      |                                                                                                                                                    
aTTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACcagcag                                                                                                                                                      <  1:159672/150‑1 (MQ=255)
      ggTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGgt                                                                                                                                                >  2:451982/1‑150 (MQ=255)
         gCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGtgtg                                                                                                                                             <  2:1233099/150‑1 (MQ=255)
            ggAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCACACCAGCCTGCTGgtggcgc                                                                                                                                           <  1:1605392/150‑5 (MQ=255)
            ggAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTgcgc                                                                                                                                          >  1:672694/1‑150 (MQ=255)
             gAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCACACCAGCCTGCTGgtggcgct                                                                                                                                          <  2:575073/150‑6 (MQ=255)
             gAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCa                                                                                                                                         <  1:1089851/150‑1 (MQ=255)
                cTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTc                                                                                                                                      <  1:359745/150‑1 (MQ=255)
                 tATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCt                                                                                                                                     <  2:973369/150‑1 (MQ=255)
                 tATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCt                                                                                                                                     <  2:918803/150‑1 (MQ=255)
                   ttGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGc                                                                                                                                   <  1:480350/150‑1 (MQ=255)
                     ggTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTc                                                                                                                                 <  1:912388/150‑1 (MQ=255)
                     ggTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTc                                                                                                                                 <  2:1606443/150‑1 (MQ=255)
                     ggTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTc                                                                                                                                 <  1:1550969/150‑1 (MQ=255)
                                  cTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGTCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                                    <  2:1619313/150‑1 (MQ=255)
                                  cTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                                    <  1:1277808/150‑1 (MQ=255)
                                    gCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGTTTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCATGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGaa                                                                                                                  >  2:246686/1‑150 (MQ=255)
                                     ccAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAAc                                                                                                                 <  2:672694/150‑1 (MQ=255)
                                      cAGCGCCACGCAGCGCAGATGTAAGGGGACATTGATTGAACTTGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACa                                                                                                                <  2:1176801/150‑1 (MQ=255)
                                        gcgcCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtg                                                                                                              <  2:1592864/150‑1 (MQ=255)
                                             aCGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTc                                                                                                         <  2:159858/150‑1 (MQ=255)
                                                     cagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACTCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGa                                                                                                 <  1:235939/150‑1 (MQ=255)
                                                     cagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGa                                                                                                 <  2:181780/150‑1 (MQ=255)
                                                     cagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGa                                                                                                 <  2:565833/150‑1 (MQ=255)
                                                           tAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGt                                                                                           >  1:503391/1‑150 (MQ=255)
                                                           tAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGt                                                                                           <  1:947935/150‑1 (MQ=255)
                                                           tAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGt                                                                                           <  1:73417/150‑1 (MQ=255)
                                                             aGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGtgt                                                                                         <  2:131252/150‑3 (MQ=255)
                                                                   cAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGa                                                                                   >  1:1405308/1‑150 (MQ=255)
                                                                       gATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCt                                                                               >  2:134365/1‑150 (MQ=255)
                                                                        tttGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCATGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                              <  1:246686/149‑1 (MQ=255)
                                                                           gAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCaa                                                                           <  2:634985/150‑1 (MQ=255)
                                                                            aaCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                                    >  1:576018/1‑108 (MQ=255)
                                                                            aaCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                                    <  2:576018/108‑1 (MQ=255)
                                                                            aaCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCaaa                                                                          >  1:964054/1‑150 (MQ=255)
                                                                             aCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAg                                                                         >  1:222422/1‑150 (MQ=255)
                                                                             aCTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAg                                                                         >  2:1566785/1‑150 (MQ=255)
                                                                                ggcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgg                                                                      >  1:439231/1‑150 (MQ=255)
                                                                                 gcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTggg                                                                     <  2:998799/150‑1 (MQ=255)
                                                                                  cggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                                    >  2:944403/1‑150 (MQ=255)
                                                                                  cggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                                    <  1:505614/150‑1 (MQ=255)
                                                                                  cggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                                    <  2:1186699/150‑1 (MQ=255)
                                                                                   ggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGa                                                                   <  1:906405/150‑1 (MQ=255)
                                                                                   ggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGa                                                                   >  1:708275/1‑150 (MQ=255)
                                                                                       aCCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACg                                                               <  1:187810/150‑1 (MQ=255)
                                                                                        ccGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGa                                                                                   <  2:432160/129‑1 (MQ=255)
                                                                                        ccGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGa                                                                                   >  1:432160/1‑129 (MQ=255)
                                                                                        ccGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGc                                                              >  2:1492945/1‑150 (MQ=255)
                                                                                              tGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAgg                                                        >  2:75186/1‑150 (MQ=255)
                                                                                              tGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAgg                                                        >  2:387294/1‑150 (MQ=255)
                                                                                               gggCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGc                                                       <  1:734635/150‑1 (MQ=255)
                                                                                                 gCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTg                                                     <  1:451982/150‑1 (MQ=255)
                                                                                                  ccGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGa                                                    <  1:1352688/150‑1 (MQ=255)
                                                                                                     tAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGAtgt                                                 >  2:1336190/1‑150 (MQ=255)
                                                                                                        tGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggt                                              >  1:1194778/1‑150 (MQ=255)
                                                                                                         gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAGGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgctg                                             <  2:193194/150‑1 (MQ=255)
                                                                                                         gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggtg                                             >  1:1619192/1‑150 (MQ=255)
                                                                                                          caccacCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGa                                            <  2:1096979/150‑1 (MQ=255)
                                                                                                                cGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACg                                      <  1:1566785/150‑1 (MQ=255)
                                                                                                                  cccGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGt                                    <  1:134365/150‑1 (MQ=255)
                                                                                                                      tGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTg                                >  2:105839/1‑150 (MQ=255)
                                                                                                                         aaCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTa                             <  1:1305559/150‑1 (MQ=255)
                                                                                                                         aaCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTAAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTa                             <  2:684548/150‑1 (MQ=255)
                                                                                                                           cGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATc                           <  2:1619192/150‑1 (MQ=255)
                                                                                                                             cTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAg                         >  2:504690/1‑150 (MQ=255)
                                                                                                                                tCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgcgc                      >  2:1500872/1‑150 (MQ=255)
                                                                                                                                tCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgcgc                      >  1:360673/1‑150 (MQ=255)
                                                                                                                                 cTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGcc                     <  1:105839/150‑1 (MQ=255)
                                                                                                                                  ttttACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCa                    <  2:503391/150‑1 (MQ=255)
                                                                                                                                  ttttACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCa                    <  1:778491/150‑1 (MQ=255)
                                                                                                                                   tttACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCaa                   >  2:251812/1‑150 (MQ=255)
                                                                                                                                       cAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAg               >  1:391244/1‑150 (MQ=255)
                                                                                                                                          cgcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCa             >  2:952526/5‑150 (MQ=255)
                                                                                                                                          cgcacaccaGCCTGCTGGTGGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCaa            >  1:1419102/5‑150 (MQ=255)
                                                                                                                                          cgcCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCaa            >  1:481293/1‑150 (MQ=255)
                                                                                                                                            cacaccaGCCTGCTGGTGGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACt          >  2:380400/3‑150 (MQ=255)
                                                                                                                                                     agGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTcgcg  >  1:912323/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTCGCG  >  NC_000913/2672491‑2672788

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: