Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,453 A→C 9.3% intergenic (‑97/+19) cutA ← / ← dcuA divalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,4530AC 9.3% 196.8 / 10.7 83intergenic (‑97/+19)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (35/40);  new base C (2/6);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

CCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCT  >  NC_000913/4365305‑4365599
                                                                                                                                                    |                                                                                                                                                  
ccGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACTCCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGa                                                                                                                                                    >  2:276236/1‑149 (MQ=255)
 cGTACATAGCACCACGACAGACGCGGAATTCGACCTTTATTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACCGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                   <  1:1151693/149‑1 (MQ=255)
    aCATAGCACCACTACAGACGCGCTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTAAACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGc                                                                                                                                                <  1:109965/149‑1 (MQ=255)
        aGCCCCACGACAGACGCGGTATTCGAACTTTTATCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGGTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGccccc                                                                                                                                            >  1:753545/1‑148 (MQ=255)
        aGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcg                                                                                                                                            >  2:419850/1‑149 (MQ=255)
         gCACCACGACAGACGCGGTATTCGAACTTTTATCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAATAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGGCCcgc                                                                                                                                           >  1:82665/1‑149 (MQ=255)
          caccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGCACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgcg                                                                                                                                          >  1:12793/1‑149 (MQ=255)
          caccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgcg                                                                                                                                          >  1:8058/1‑149 (MQ=255)
          caccacGACAGACGCGGTAATCGAACTTTTTTCATCAAGCATACACACCTCGCAAGAACAGCTGGACCCGCAAGAGAGAGGGCCCTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgct                                                                                                                                          <  1:502658/149‑2 (MQ=255)
           accacGACAGACGCGGTATTCGAACTTTTTTCATCCAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGCCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGa                                                                                                                                         <  1:448048/149‑1 (MQ=255)
           accacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctt                                                                                                                                         <  2:867154/149‑3 (MQ=255)
            ccacGACAGACGCGGTATTCGAAGTTTTTTCATCATGCATAAACAGCTCCCAAGAACAGATCGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                                        <  1:573008/149‑4 (MQ=255)
            ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                                        <  1:340134/149‑4 (MQ=255)
             cacGACAGCCGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTTGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                       >  2:183884/1‑149 (MQ=255)
             cacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttc                                                                                                                                       <  1:682664/149‑5 (MQ=255)
             cacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                       >  2:730298/1‑149 (MQ=255)
              acGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATTCAACGGACGGGGGAAGGTTGTTAGTAGGGGGGATAGATAGCGGGGCGTGATCGCGCCGCGATTg                                                                                                                                      >  2:564193/1‑149 (MQ=255)
               cGACAGACGCGGTATTCGAACTTTTTTCATCCAGCATAAACACCTCGCAAGAACAGATGGATCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGTAAAGAAAGCGGGGCGTTCACGCCCcgctttctt                                                                                                                                     <  1:510142/149‑7 (MQ=255)
                gACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGAt                                                                                                                                    <  2:41072/149‑1 (MQ=255)
                  cAGACGCGGTATTCGAACTTTTTTCATCAAGCTTAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCa                                                                                                                                  <  2:824828/149‑1 (MQ=255)
                  cAGACGCGGTAATCGAACTTATTTCATCAAGCATAAACCCCTCGCAAGAACAGAGGGAACCGCAAGAGGGAGGTCACAATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCa                                                                                                                                  <  1:457128/149‑1 (MQ=255)
                     aCGCGGTATTCCAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGCAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTa                                                                                                                               <  1:467679/149‑1 (MQ=255)
                        cgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGAAGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGGTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                            <  2:92032/149‑1 (MQ=255)
                         ggTATTCGACCTTTTTTCATCAAGCATAAACACCTCGCAAGACCAGATGGAACCGCAAGAGAGAGGTGACTATACAACGGACGGGGGAAGGAGGGTAGTCGGCGGGTAAGACAGCGGGGCGTGAAAGCCGCGGGAGTGATCAGTAGAGc                                                                                                                           >  2:788860/1‑149 (MQ=255)
                            aTTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTAACGCCCCGCGAATGATCATTACAGCATg                                                                                                                        >  2:83198/1‑149 (MQ=255)
                              tCGAACTTTTTTCCTCAAGCATCAACTCCTCGCAAGAACAGATGGAACAGCAATAGAGCGGTCACTCTACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGaa                                                                                                                      <  2:809561/149‑1 (MQ=255)
                                gAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                    >  2:625696/1‑149 (MQ=255)
                                 aaCTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACCACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCTTGAACGCCCCGCGATTGATCAATACAGCATTAAGCt                                                                                                                   >  2:1026304/1‑149 (MQ=255)
                                 aaCTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGTGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCt                                                                                                                   <  1:92217/149‑1 (MQ=255)
                                    ttttttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACAGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTAcc                                                                                                                <  2:1141025/149‑1 (MQ=255)
                                           atcaAGCATAAACACCTCGCAAGCACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTGAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGTGCATTTCAGCATGAAGCTACCCAGCAcg                                                                                                         >  1:671789/1‑149 (MQ=255)
                                            tcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTCCAGCATGAAGCTACCCAGCAcga                                                                                                        >  2:96372/1‑149 (MQ=255)
                                            tcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcga                                                                                                        >  1:573185/1‑149 (MQ=255)
                                               atcATAAACACCTAGCAAGAACAGAAGGAACCGCATGAGAGAGGTCACTTTACTTCGGAGGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagc                                                                                                     <  1:989756/147‑1 (MQ=255)
                                                      acacCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGAGGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAGTGATCATTACAGCATGAAGCGACCCAGCACGAAGCCGAAGCa                                                                                              >  2:569155/1‑149 (MQ=255)
                                                             gCAAGAACAGATGGAACCGAAAGAGAGAGGTCACTATAAAACGGACGTGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCTATTGATCATTACAGCATGAAGCTACCCATCACGAAGCCGAAGCAAACGGcc                                                                                       >  2:868551/1‑149 (MQ=255)
                                                                aGAACACAGGGAACCGCAAGAGAGAGGGCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgg                                                                                    <  1:638689/149‑1 (MQ=255)
                                                                 gAACAGATGGCACCGCCAGAGCGAGGGCACTCTCCACCGGCGGGGGGATGGATGTTAGGCGGCGGGAATGAAAGCGGGGCGGGAACGCCCCGCGATTGATCCTTACTGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAggg                                                                                   <  1:722491/149‑1 (MQ=255)
                                                                  aaCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGTGGGATGGATGGTAGGCTGCGGGAAAGAAAGCGGGGCGTGATCGCCGCGCGATTGATCATTACCTCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                                  <  1:802658/149‑1 (MQ=255)
                                                                  aaCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                                  >  1:962095/1‑149 (MQ=255)
                                                                    cctaTGGACCAGCCAGCGCGAGGTCACTATACACCGGACGGGGGAAGGATGTTAGAAGGCGGGAAAGAAAGCCGCGCGTGAACGCCCAGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCaa                                                                                <  2:786621/146‑1 (MQ=255)
                                                                    cAGATCGATCCGCAAGAGATAGGTAACTATACACGGGACGGGGGAAGGATGTTAGTCGGCGTGAAAGAAAGCGGGGCGTGAACGCCGCGCGATTGATCATTACAGCATGATGCTACCCAGCACGAAGCCTAAGCAAACGGCCAGGGCat                                                                                <  2:985209/149‑2 (MQ=255)
                                                                         ggAACCGCAAGTGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGGGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                                           <  2:131666/149‑1 (MQ=255)
                                                                         ggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCCAACGGCCAGGGCAACAccc                                                                           >  1:752134/1‑149 (MQ=255)
                                                                         ggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                                           >  2:492380/1‑149 (MQ=255)
                                                                            aCCGCAAGAGAGAGGTCACTATACACCGGACGGGGGAAGGAAGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCaga                                                                        <  1:1094453/149‑1 (MQ=255)
                                                                            aCCGCAAGAGAGAGGGCACTATCCAACGGACGGGGGAAGGACGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCACGGCAACACCCaga                                                                        <  1:725086/149‑1 (MQ=255)
                                                                              cGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCATTTCAGCATGAAGCTCCCCAGCACGCAGCCGCAGCAAACGGCCAGGGCGACAccaggagg                                                                      >  2:1101219/1‑143 (MQ=255)
                                                                              cGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGt                                                                      >  1:846225/1‑149 (MQ=255)
                                                                               gCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTa                                                                     >  2:506549/1‑149 (MQ=255)
                                                                                cAAGAGAGAGGTCACTTTTCTACGGACGGGGGAAGGATGTTAGTCGGAGGGAATGAAAGCGGGGCGTGAACGCCCCGCTATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACTGCCAGGGCAACACCCAGAGTag                                                                    <  2:1182066/149‑2 (MQ=255)
                                                                                   gagagagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCg                                                                 >  2:839416/1‑149 (MQ=255)
                                                                                    agagagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCgg                                                                <  1:758356/149‑1 (MQ=255)
                                                                                     cagagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCggg                                                               <  1:513386/148‑1 (MQ=255)
                                                                                              cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAAc                                                      <  1:276236/149‑1 (MQ=255)
                                                                                               aataCACCGGCCGGGGGAACGATGTTATTCGGCCGCAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg                                                     <  1:765374/148‑1 (MQ=255)
                                                                                                 taCAACGGACGGGGGAAGGATGTTAGACGGCGGGAAAGAACGCCGGGCGTGAACGCCCCGCGATTGATCATTAAAGCATGAAGCTACCCAGCACGAAGCCGAACCAAACGCCCAGCGCAACACCCAGAGTACCCGGGAAGAAGACCGGa                                                   >  1:296990/1‑149 (MQ=255)
                                                                                                    aaCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCCCGAAGCCGAAGCAAACGGCCAGGGCAACCCCCAGAGTACCCGGGATGAAGAACGGATgg                                                >  2:297048/1‑149 (MQ=255)
                                                                                                    aaCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATgg                                                >  1:759030/1‑149 (MQ=255)
                                                                                                     acggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGt                                               <  1:43721/149‑1 (MQ=255)
                                                                                                     acggacggGGGAAGGATGTTAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGt                                               <  1:208679/149‑1 (MQ=255)
                                                                                                       ggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                             <  1:466449/149‑1 (MQ=255)
                                                                                                          cggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAg                                          <  1:394759/149‑1 (MQ=255)
                                                                                                             gggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACCGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGTTGAAGCCCGGATCGTTGCAGACg                                       >  1:1054615/1‑149 (MQ=255)
                                                                                                                 aGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAAtt                                   <  2:963036/149‑1 (MQ=255)
                                                                                                                    aTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                                >  1:376017/1‑149 (MQ=255)
                                                                                                                         aGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAACCGGCCAGGGCACCACCCAGAGTACCCGGGAAGAAGAACGGATCGTTGAAGACGAATATACCGATa                           >  1:266389/1‑149 (MQ=255)
                                                                                                                             ggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGTTTGAAGACGAATTTACCGATACggg                       >  2:919380/1‑149 (MQ=255)
                                                                                                                             ggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACggg                       <  2:486719/149‑1 (MQ=255)
                                                                                                                             ggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACggg                       >  1:486719/1‑149 (MQ=255)
                                                                                                                              gcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATAATTACTCCATGAAGATACCCAGCAGGAAGCCGAAGCAAAAGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgt                      <  2:713604/149‑1 (MQ=255)
                                                                                                                                  gaaagaaaGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                                <  2:422183/135‑1 (MQ=255)
                                                                                                                                  gaaagaaaGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                                >  1:422183/1‑135 (MQ=255)
                                                                                                                                  gaaagaaaGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgtagta                  <  1:183884/149‑1 (MQ=255)
                                                                                                                                     agaaagCGGGGAGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAACCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTAccc               <  2:753545/149‑1 (MQ=255)
                                                                                                                                       atagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt             <  1:189693/147‑1 (MQ=255)
                                                                                                                                       aaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAATCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt             <  1:730298/149‑1 (MQ=255)
                                                                                                                                       aaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAGGAAGCACGGAGGGTTGAAGCCGAATTTACCGATACGGGTAGTACCCgg             >  1:557641/1‑148 (MQ=255)
                                                                                                                                         agCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGCGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtgg           >  1:1101850/1‑148 (MQ=255)
                                                                                                                                          gCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTc          >  2:231175/1‑149 (MQ=255)
                                                                                                                                           cGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGACGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGCCGAATTTACCGATACGGGTAGTCCCCGTGTca         >  1:827124/1‑149 (MQ=255)
                                                                                                                                                agTGACCGCTCCGTGATTGATCATTACAGCATGAAGCTACCCCGCATTAAGCCGAAGCAAACGGCCAGGGCAACACACTTAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcttccat    <  2:82665/148‑1 (MQ=255)
                                                                                                                                                 gTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatccatc   <  2:394356/149‑1 (MQ=255)
                                                                                                                                                  tGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCt  <  2:12793/149‑1 (MQ=255)
                                                                                                                                                  ggAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACTGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCt  <  2:8058/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCT  >  NC_000913/4365305‑4365599

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