Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,672,639 A→C 5.9% Q198P (CAG→CCG) ‡ 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,6390AC 5.9% 262.7 / 7.5 85Q198P (CAG→CCG) ‡hcaDputative 3‑phenylpropionate/cinnamate dioxygenase ferredoxin reductase subunit
Reads supporting (aligned to +/- strand):  ref base A (42/38);  new base C (3/2);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

GATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTC  >  NC_000913/2672490‑2672785
                                                                                                                                                      |                                                                                                                                                  
gATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACcagca                                                                                                                                                    <  1:400898/150‑1 (MQ=255)
        gTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGgtg                                                                                                                                            <  1:456542/150‑1 (MQ=255)
          gCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCCTGCTGgtggc                                                                                                                                          <  2:1249991/150‑3 (MQ=255)
           ccGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCACACCAGCCTGCTGgtggc                                                                                                                                          <  1:1153248/150‑3 (MQ=255)
            cGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTgcg                                                                                                                                        <  1:1971119/150‑1 (MQ=255)
              gAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCa                                                                                                                                      >  2:799424/1‑150 (MQ=255)
                   aTTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTg                                                                                                                                 >  2:1826333/1‑150 (MQ=255)
                         cTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaat                                                                                                                           >  2:1289447/1‑150 (MQ=255)
                          tGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaata                                                                                                                          >  2:1707733/1‑150 (MQ=255)
                          tGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaata                                                                                                                          >  1:1598189/1‑150 (MQ=255)
                          tGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCaata                                                                                                                          >  1:1204751/1‑150 (MQ=255)
                              aCTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGc                                                                                                                      <  2:1808322/150‑1 (MQ=255)
                              aCTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGc                                                                                                                      <  2:1473074/150‑1 (MQ=255)
                                 ggCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCAt                                                                                                                   >  1:1399900/1‑150 (MQ=255)
                                   cTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                                 >  2:1633994/1‑150 (MQ=255)
                                   cTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTg                                                                                                                 <  2:388984/150‑1 (MQ=255)
                                      ccAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAAc                                                                                                              <  2:1070497/150‑1 (MQ=255)
                                       cAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACa                                                                                                             <  2:754846/150‑1 (MQ=255)
                                       cAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACa                                                                                                             <  1:1981097/150‑1 (MQ=255)
                                         gcgcCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtg                                                                                                           >  1:1530898/1‑150 (MQ=255)
                                           gcCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgtg                                                                                                         >  1:703216/1‑150 (MQ=255)
                                           gcCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACAtgtg                                                                                                         <  2:19698/150‑1 (MQ=255)
                                              aCGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTc                                                                                                      <  2:1204751/150‑1 (MQ=255)
                                              aCGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTc                                                                                                      <  2:1530898/150‑1 (MQ=255)
                                                gcagcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGa                                                                                                    >  1:475759/1‑150 (MQ=255)
                                                 cagcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGcc                                                                                                                     >  1:1182539/1‑132 (MQ=255)
                                                 cagcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGcc                                                                                                                     <  2:1182539/132‑1 (MQ=255)
                                                  agcgcagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATg                                                                                                  <  2:843474/150‑1 (MQ=255)
                                                      cagATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGa                                                                                              <  1:799424/150‑1 (MQ=255)
                                                          tGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaaa                                                                                          >  2:407156/1‑150 (MQ=255)
                                                          tGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaaa                                                                                          >  2:546895/1‑150 (MQ=255)
                                                          tGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGaaaaa                                                                                          <  1:1707733/150‑1 (MQ=255)
                                                           gTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAg                                                                                         <  2:1525903/150‑1 (MQ=255)
                                                             aaGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTa                                                                                       >  2:605098/1‑150 (MQ=255)
                                                             aaGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTa                                                                                       >  2:642736/1‑150 (MQ=255)
                                                                 tGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAAc                                                                                   >  1:1926725/1‑150 (MQ=255)
                                                                   aCAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTg                                                                                 >  2:1993863/1‑150 (MQ=255)
                                                                    cAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGa                                                                                >  1:1915101/1‑150 (MQ=255)
                                                                     aGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGAc                                                                               >  2:745996/1‑150 (MQ=255)
                                                                      gTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACg                                                                              >  1:747204/1‑150 (MQ=255)
                                                                       tGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGACCTGACGc                                                                             >  2:1028482/1‑150 (MQ=255)
                                                                        gATTGAACTGGCCGCAACCGTCATGGGCCGTAATGCACCAACGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCt                                                                            >  2:892137/1‑150 (MQ=255)
                                                                         aTTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTg                                                                           <  1:52962/150‑1 (MQ=255)
                                                                          ttGAACTGGCGGCACCCGTCATGGGCCGTAATGCACCACCGACCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGc                                                                          >  1:1217568/1‑150 (MQ=255)
                                                                          ttGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGc                                                                          <  1:1683617/150‑1 (MQ=255)
                                                                            gAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCaa                                                                        >  2:716695/1‑150 (MQ=255)
                                                                               cTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGt                                                                     >  1:1702353/1‑150 (MQ=255)
                                                                                 ggcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgg                                                                   >  2:1587253/1‑150 (MQ=255)
                                                                                  gcggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTggg                                                                  <  2:36439/150‑1 (MQ=255)
                                                                                   cggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTgggg                                                                 <  1:1289447/150‑1 (MQ=255)
                                                                                    ggcAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGGCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGa                                                                <  1:1353168/150‑1 (MQ=255)
                                                                                        aCCGTCATGGGCCGTAATGCACCACCGCCCGTCCAACGCTATCTTTTACAGCGCCACCAGCAGGCTCGTGTGCGCATTCTGCTCAAAAATGCCATTGAACATGTGGTCGATGCTGAAAAAGTAGAACTGACGCTGCAAAGAGGGGAAAcc                                                            >  2:1363386/1‑149 (MQ=255)
                                                                                         ccGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGc                                                           >  1:230973/1‑150 (MQ=255)
                                                                                         ccGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGc                                                           >  1:1165157/1‑150 (MQ=255)
                                                                                         ccGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGc                                                           >  1:626359/1‑150 (MQ=255)
                                                                                                ggtccctAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGAGGTCGATGGTGAAAAAGTAGAACTCACGCTGCAAAGTGGGGAAAAGCTTCACGc                                                    >  2:70716/7‑150 (MQ=255)
                                                                                                  gCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTg                                                  <  2:1399900/150‑1 (MQ=255)
                                                                                                  gCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTg                                                  >  1:1367772/1‑150 (MQ=255)
                                                                                                      tAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGAtgt                                              >  1:1349831/1‑150 (MQ=255)
                                                                                                      tAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGAtgt                                              <  1:443801/150‑1 (MQ=255)
                                                                                                        aTGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtgg                                            >  1:1000993/1‑150 (MQ=255)
                                                                                                          gCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTATAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATgtggtg                                          <  2:747204/150‑1 (MQ=255)
                                                                                                                ccGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTAc                                    <  1:663876/150‑1 (MQ=255)
                                                                                                                 cGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACg                                   <  2:1165157/150‑1 (MQ=255)
                                                                                                                  gCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACgg                                  <  2:878565/150‑1 (MQ=255)
                                                                                                                       tGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTg                             <  2:556185/150‑1 (MQ=255)
                                                                                                                       tGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTg                             <  2:1915101/150‑1 (MQ=255)
                                                                                                                       tGCAACGCTATATTATACCGCGCCACCAGCAGGCTTGTGTGCGCATTCTGCTCTATAATGCCATTGAAATTGTGTTCGATGGTGAAAAAGTATAACTGACGCTTCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTTATATACGGTATTg                             <  2:512470/150‑1 (MQ=255)
                                                                                                                        gCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTgg                            <  1:1826333/150‑1 (MQ=255)
                                                                                                                        gCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTgg                            <  1:70716/150‑1 (MQ=255)
                                                                                                                        gCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTgg                            <  1:1608537/150‑1 (MQ=255)
                                                                                                                        gCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTgg                            <  2:1217568/150‑1 (MQ=255)
                                                                                                                           aCGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTAt                         >  2:552784/1‑150 (MQ=255)
                                                                                                                                aTCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAgcg                    >  1:756856/1‑150 (MQ=255)
                                                                                                                                   ttttACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCa                 >  1:772575/1‑150 (MQ=255)
                                                                                                                                     ttACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAAc               <  1:745996/150‑1 (MQ=255)
                                                                                                                                      tACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACg              >  2:304656/1‑150 (MQ=255)
                                                                                                                                      tACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACg              >  2:1949560/1‑150 (MQ=255)
                                                                                                                                      tACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACg              <  2:475759/150‑1 (MQ=255)
                                                                                                                                        cAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGGGGGCGCTGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAg            <  2:970399/150‑1 (MQ=255)
                                                                                                                                           cgcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCa          >  1:1899991/5‑150 (MQ=255)
                                                                                                                                           cgcacaccaGCCTGCTGGTGGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCaa         >  2:9251/5‑150 (MQ=255)
                                                                                                                                            gcacaccaGCCTGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCaa         >  1:1789402/4‑150 (MQ=255)
                                                                                                                                             ccaccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTggg                                                                  >  2:1947025/1‑91 (MQ=255)
                                                                                                                                             ccaccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTggg                                                                  <  1:1947025/91‑1 (MQ=255)
                                                                                                                                             ccaccaGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACt       >  1:813887/1‑150 (MQ=255)
                                                                                                                                                   agcagGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTc  >  1:709331/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GATTATCGGTGCCGGAACTATTGGTCTGGAACTGGCTGCCAGCGCCACGCAGCGCAGATGTAAGGTGACAGTGATTGAACTGGCGGCAACCGTCATGGGCCGTAATGCACCACCGCCCGTGCAACGCTATCTTTTACAGCGCCACCAGCAGGCTGGTGTGCGCATTCTGCTCAATAATGCCATTGAACATGTGGTCGATGGTGAAAAAGTAGAACTGACGCTGCAAAGTGGGGAAACGCTTCAGGCTGATGTGGTGATTTACGGTATTGGTATCAGCGCCAACGAGCAACTGGCTC  >  NC_000913/2672490‑2672785

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: