Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,453 A→C 7.4% 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 7.4% 380.2 / 17.3 134intergenic (‑97/+19)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (52/72);  new base C (4/6);  total (56/78)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

Only 100 of 138 total aligned reads displayed.

GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGTA  >  NC_000913/4365304‑4365602
                                                                                                                                                      |                                                                                                                                                     
gCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGa                                                                                                                                                       >  1:1195561/1‑150 (MQ=255)
 ccGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACCACGGACGGGGGAAGGATGGTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                      <  2:1539087/150‑1 (MQ=255)
 ccGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGTGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                      <  2:692082/150‑1 (MQ=255)
    tACATAGCACCACGACAGACGCGGGTTTCTAACTTTTTTCATCAAGCCTAAACACCTCGCAAGAACCGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGc                                                                                                                                                   <  2:923293/150‑1 (MQ=255)
         aGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCAAAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgc                                                                                                                                              <  1:1125049/150‑1 (MQ=255)
            accacGACAGACGCTGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                                           <  2:1305804/150‑4 (MQ=255)
            accacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                                           <  1:1210969/150‑4 (MQ=255)
             ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                          >  2:626369/1‑150 (MQ=255)
             ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGCAACCGCAAGAGAGAGGGCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                          <  1:1016897/150‑1 (MQ=255)
              cacGACAGACTCGGTATTCTAACTTTTATCATCAAGCATAAAGACCTCGCAAGAACCGATCGAACCCCAAGAGCGAGGTCACTATACAACGGACCGGGGACGCATGTGAGGCGGCGGGAAAGAAAGCGGGGCGCTCACGCCCcgctttct                                                                                                                                         <  2:300333/150‑6 (MQ=255)
              cacGACAGACGCGGTATTCGAACTTTTTTCATCACGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGTGAGAGGTCAATATACAACGGACGGGGGATGTATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttct                                                                                                                                         <  1:1731769/150‑6 (MQ=255)
              cacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGAGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTg                                                                                                                                         >  1:808673/1‑150 (MQ=255)
               acGACAGACGCTGTATTCGAACTTTTTTCATGAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttct                                                                                                                                         <  1:1018560/150‑6 (MQ=255)
               acGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGa                                                                                                                                        <  1:872414/150‑1 (MQ=255)
                 gACAGACGCGGTATTCGAACTTTTTTCAGCAAGCATAACCACCTCGCAAGAACAGATGGAACCGAAAGAGATAGGGCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATc                                                                                                                                      <  2:1416064/150‑1 (MQ=255)
                  aCAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCa                                                                                                                                     <  2:415890/150‑1 (MQ=255)
                   cAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGCGAGAGGTCACTCTACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                                                                    <  1:812505/150‑1 (MQ=255)
                   cAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGATAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                                                                    >  1:1311463/1‑150 (MQ=255)
                    aGACGCGTTATTCGAAATTTTTTCATCAAGCATAAACCCCTCGCAAGAACAGATGGCACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                                   <  1:120230/150‑1 (MQ=255)
                    aGACGCGGTATTCGAACTTTTATCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                                   <  1:1723075/150‑1 (MQ=255)
                    aGACCCGGTATTCGAACTTTATTCATCACGCATAAACACCTCGCACGAACAGCTGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGGCGGCGCGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                                   <  1:299604/150‑1 (MQ=255)
                     gACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTa                                                                                                                                  <  1:668450/150‑1 (MQ=255)
                       cgcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACa                                                                                                                                >  1:414406/1‑150 (MQ=255)
                        gcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                               <  1:1035401/150‑1 (MQ=255)
                        gcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                               <  1:1818508/150‑1 (MQ=255)
                         cgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGc                                                                                                                              <  2:758065/150‑1 (MQ=255)
                          ggTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAACAACAGCTGGAACCGCAAGAGAGAGGTCACAATACAACGGACGGGGGAAGGATGTTAGTCGGCGGTAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCa                                                                                                                             <  2:1395006/150‑1 (MQ=255)
                              ttCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGaa                                                                                                                         <  1:1118518/150‑1 (MQ=255)
                               tCGAACTTTTTTCATCAAGCAAAAACACCTCGCAAGAACAGATGTAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                                                                        <  1:1567607/150‑1 (MQ=255)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                       >  2:198318/1‑150 (MQ=255)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                       <  1:795028/150‑1 (MQ=255)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                       <  1:228356/150‑1 (MQ=255)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGAATGAAGc                                                                                                                       >  1:1641011/1‑150 (MQ=255)
                                cGAAATTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                       <  1:634315/150‑1 (MQ=255)
                                      tttttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCa                                                                                                                 <  1:793195/150‑1 (MQ=255)
                                        tttCATCAAGCATAAACACCTCGCAAGAACAGATGTAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                                                                               <  1:888103/150‑1 (MQ=255)
                                           catcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                                                                            <  1:246734/150‑1 (MQ=255)
                                            atcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcga                                                                                                           >  2:1714790/1‑150 (MQ=255)
                                                aGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCCGCAcgcagcc                                                                                                       >  2:712992/1‑150 (MQ=255)
                                                    tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                   >  1:273722/1‑150 (MQ=255)
                                                    tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                   <  2:210930/150‑1 (MQ=255)
                                                    tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACAATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                   <  1:128969/150‑1 (MQ=255)
                                                     aaaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaagc                                                                                                  <  2:1690070/150‑1 (MQ=255)
                                                      aaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCa                                                                                                 >  1:1192798/1‑150 (MQ=255)
                                                      aaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCa                                                                                                 >  2:1344821/1‑150 (MQ=255)
                                                        cacCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCaaa                                                                                               >  1:1349134/1‑150 (MQ=255)
                                                         acCTCGCAAGAACAGATGGAACCGCAAGATAGAGGTCACTATACTACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAAc                                                                                              <  2:1169026/150‑1 (MQ=255)
                                                          ccTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACg                                                                                             >  1:1160178/1‑150 (MQ=255)
                                                           cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGCTTGCTCATTACAGCATGCCGCTACCCAGCCCGCCGCCGAAGCCAACgg                                                                                            >  1:986195/1‑150 (MQ=255)
                                                                  taaCAGATGGAACGGCAAGAGAGAGGTCACTATACAACGGAGGGGGGAAGGATGTAAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCTGGGc                                                                                     <  1:169971/149‑1 (MQ=255)
                                                                  gAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                                     <  2:763422/150‑1 (MQ=255)
                                                                   aaCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCa                                                                                    <  2:564566/150‑1 (MQ=255)
                                                                      aGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAaca                                                                                 >  2:41598/1‑150 (MQ=255)
                                                                       gATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAacac                                                                                >  2:506732/1‑150 (MQ=255)
                                                                         tGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                                              <  1:1175719/150‑1 (MQ=255)
                                                                          ggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                                                                             >  2:952166/1‑150 (MQ=255)
                                                                          ggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                                                                             >  2:1224996/1‑150 (MQ=255)
                                                                          ggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                                                                             >  1:801222/1‑150 (MQ=255)
                                                                           gAACCGCAAGAGTGAGGTCACTCTACACCGGACGGGTTAAGGATGTTAGTCGGCGGGAATGAAAGCGGGGCGTGAACGCCCCGAGATTGATCATTACATCATTAAGCTACCCAGCACTTAGCCGAAGCAAACGGCCAGGGCAACACCCag                                                                            <  1:560329/150‑1 (MQ=255)
                                                                             aCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                                                                          <  2:949044/150‑1 (MQ=255)
                                                                              ccGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAATACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGt                                                                         >  1:1738389/1‑150 (MQ=255)
                                                                              ccGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATCTTAGTCGGCGGGAAAGAAAGCGGCGCGTGACCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCCCGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGt                                                                         >  1:360869/1‑150 (MQ=255)
                                                                               cGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTa                                                                        >  1:903137/1‑150 (MQ=255)
                                                                               cGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTa                                                                        <  2:1257603/150‑1 (MQ=255)
                                                                                gCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGACGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAc                                                                       >  1:975587/1‑150 (MQ=255)
                                                                                  aaGAGAGAGGTCACTATACAACGGACGGGGGTAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAccc                                                                     >  1:1450121/1‑150 (MQ=255)
                                                                                     agagagGTCACTATACAACGGACGGGGGAAGGATGTAAGTCGGCGGGACAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACATCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCggg                                                                  <  2:1652387/150‑1 (MQ=255)
                                                                                            tCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCCATAGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaaga                                                           >  1:1518461/1‑150 (MQ=255)
                                                                                             cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaagaa                                                          >  2:1578896/1‑150 (MQ=255)
                                                                                              aCTATCCAACGGACGGGGGAAGGAGGTTAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAAc                                                         <  2:1434006/150‑1 (MQ=255)
                                                                                               cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg                                                        >  1:449154/1‑150 (MQ=255)
                                                                                                tataCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACgg                                                       <  2:425497/150‑1 (MQ=255)
                                                                                                 ataCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGa                                                      <  1:1235459/150‑1 (MQ=255)
                                                                                                    cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCACCACCCAGAGTACCCGGGATGAAGAACGGATgg                                                   >  2:123987/1‑150 (MQ=255)
                                                                                                    cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATgg                                                   >  2:1361503/1‑150 (MQ=255)
                                                                                                     aaCGGACGGGGGAAGGATGATAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGt                                                  <  1:143551/150‑1 (MQ=255)
                                                                                                      acggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCACGGCAACACCCAGAGTACCCTGCATCAAGCACGGATGGtt                                                 >  2:597762/1‑150 (MQ=255)
                                                                                                       cggacggGGGATGGATGTTTGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTTCAGCATGTATCTACCCAGCACGAAGGCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                                <  1:1564480/150‑1 (MQ=255)
                                                                                                       cggacggGGGAAGGATGTTAGTCGGCGGGAAAGACAGCGGTGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAGGGTTg                                                <  1:1692006/150‑1 (MQ=255)
                                                                                                       cggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                                <  1:190314/150‑1 (MQ=255)
                                                                                                          acggGGGAAGTATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAg                                             <  1:164582/150‑1 (MQ=255)
                                                                                                              ggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACCCCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACg                                          >  1:1449159/1‑150 (MQ=255)
                                                                                                               gggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGa                                         >  2:170303/1‑150 (MQ=255)
                                                                                                                      aTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAcc                                  <  2:1417078/150‑1 (MQ=255)
                                                                                                                            gTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACg                            >  1:585337/1‑150 (MQ=255)
                                                                                                                              cggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACggg                          <  2:1761098/150‑1 (MQ=255)
                                                                                                                               ggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgt                         <  2:819992/150‑1 (MQ=255)
                                                                                                                                       agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg                 >  2:254164/1‑150 (MQ=255)
                                                                                                                                       agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg                 >  1:420419/1‑150 (MQ=255)
                                                                                                                                       agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg                 >  2:662611/1‑150 (MQ=255)
                                                                                                                                       agaaagCGGGGCGTGAACGCCCAGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg                 <  2:117934/150‑1 (MQ=255)
                                                                                                                                        gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt                >  2:1564368/1‑150 (MQ=255)
                                                                                                                                        gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt                <  1:606337/150‑1 (MQ=255)
                                                                                                                                        gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCAATACGGGTAGTACCCgt                >  2:800995/1‑150 (MQ=255)
                                                                                                                                         atcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGTTGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg               >  2:703653/4‑150 (MQ=255)
                                                                                                                                         aaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg               <  1:950209/150‑1 (MQ=255)
                                                                                                                                                 gCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatccat       <  2:273722/150‑1 (MQ=255)
                                                                                                                                                  cGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatccatc      >  1:1760004/1‑150 (MQ=255)
                                                                                                                                                    tGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTg    <  2:1234969/150‑1 (MQ=255)
                                                                                                                                                      aaCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGTa  >  1:1363839/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGTA  >  NC_000913/4365304‑4365602

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: