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.

GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGTA  >  NC_000913/4365304‑4365602
                                                                                                                                                      |                                                                                                                                                     
gCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGa                                                                                                                                                       >  1:1195561/1‑150 (MQ=255)
 ccGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACCACGGAC‑GGGGGAAGGATGGTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                      <  2:1539087/150‑1 (MQ=255)
 ccGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGTGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                      <  2:692082/150‑1 (MQ=255)
    tACATAGCACCACGACAGACGCGGGTTTCTAACTTTTTTCATCAAGCCTAAACACCTCGCAAGAACCGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGc                                                                                                                                                   <  2:923293/150‑1 (MQ=255)
         aGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCAAAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgc                                                                                                                                              <  1:1125049/150‑1 (MQ=255)
            accacGACAGACGCTGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                                           <  2:1305804/150‑4 (MQ=255)
            accacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                                           <  1:1210969/150‑4 (MQ=255)
             ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                          >  2:626369/1‑150 (MQ=255)
             ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGCAACCGCAAGAGAGAGGGCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                          <  1:1016897/150‑1 (MQ=255)
              cacGACAGACTCGGTATTCTAACTTTTATCATCAAGCATAAAGACCTCGCAAGAACCGATCGAACCCCAAGAGCGAGGTCACTATACAACGGAC‑CGGGGACGCATGTGAGGCGGCGGGAAAGAAAGCGGGGCGCTCACGCCCcgctttct                                                                                                                                         <  2:300333/150‑6 (MQ=255)
              cacGACAGACGCGGTATTCGAACTTTTTTCATCACGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGTGAGAGGTCAATATACAACGGAC‑GGGGGATGTATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttct                                                                                                                                         <  1:1731769/150‑6 (MQ=255)
              cacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGAGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTg                                                                                                                                         >  1:808673/1‑150 (MQ=255)
               acGACAGACGCTGTATTCGAACTTTTTTCATGAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttct                                                                                                                                         <  1:1018560/150‑6 (MQ=255)
               acGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGa                                                                                                                                        <  1:872414/150‑1 (MQ=255)
                 gACAGACGCGGTATTCGAACTTTTTTCAGCAAGCATAACCACCTCGCAAGAACAGATGGAACCGAAAGAGATAGGGCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATc                                                                                                                                      <  2:1416064/150‑1 (MQ=255)
                  aCAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCa                                                                                                                                     <  2:415890/150‑1 (MQ=255)
                   cAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGCGAGAGGTCACTCTACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                                                                    <  1:812505/150‑1 (MQ=255)
                   cAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGATAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                                                                    >  1:1311463/1‑150 (MQ=255)
                    aGACGCGTTATTCGAAATTTTTTCATCAAGCATAAACCCCTCGCAAGAACAGATGGCACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                                   <  1:120230/150‑1 (MQ=255)
                    aGACGCGGTATTCGAACTTTTATCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                                   <  1:1723075/150‑1 (MQ=255)
                    aGACCCGGTATTCGAACTTTATTCATCACGCATAAACACCTCGCACGAACAGCTGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGGCGGCGCGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                                   <  1:299604/150‑1 (MQ=255)
                     gACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTa                                                                                                                                  <  1:668450/150‑1 (MQ=255)
                       cgcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACa                                                                                                                                >  1:414406/1‑150 (MQ=255)
                        gcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                               <  1:1035401/150‑1 (MQ=255)
                        gcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                               <  1:1818508/150‑1 (MQ=255)
                         cgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGc                                                                                                                              <  2:758065/150‑1 (MQ=255)
                          ggTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAACAACAGCTGGAACCGCAAGAGAGAGGTCACAATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGTAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCa                                                                                                                             <  2:1395006/150‑1 (MQ=255)
                              ttCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGaa                                                                                                                         <  1:1118518/150‑1 (MQ=255)
                               tCGAACTTTTTTCATCAAGCAAAAACACCTCGCAAGAACAGATGTAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                                                                        <  1:1567607/150‑1 (MQ=255)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                       >  2:198318/1‑150 (MQ=255)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                       <  1:795028/150‑1 (MQ=255)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                       <  1:228356/150‑1 (MQ=255)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGAATGAAGc                                                                                                                       >  1:1641011/1‑150 (MQ=255)
                                cGAAATTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                       <  1:634315/150‑1 (MQ=255)
                                      tttttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCa                                                                                                                 <  1:793195/150‑1 (MQ=255)
                                        tttCATCAAGCATAAACACCTCGCAAGAACAGATGTAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                                                                               <  1:888103/150‑1 (MQ=255)
                                           catcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                                                                            <  1:246734/150‑1 (MQ=255)
                                            atcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcga                                                                                                           >  2:1714790/1‑150 (MQ=255)
                                                aGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCCGCAcgcagcc                                                                                                       >  2:712992/1‑150 (MQ=255)
                                                    tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                   >  1:273722/1‑150 (MQ=255)
                                                    tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                   <  2:210930/150‑1 (MQ=255)
                                                    tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACAATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                   <  1:128969/150‑1 (MQ=255)
                                                     aaaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaagc                                                                                                  <  2:1690070/150‑1 (MQ=255)
                                                      aaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCa                                                                                                 >  1:1192798/1‑150 (MQ=255)
                                                      aaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCa                                                                                                 >  2:1344821/1‑150 (MQ=255)
                                                        cacCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCaaa                                                                                               >  1:1349134/1‑150 (MQ=255)
                                                         acCTCGCAAGAACAGATGGAACCGCAAGATAGAGGTCACTATACTACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAAc                                                                                              <  2:1169026/150‑1 (MQ=255)
                                                          ccTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACg                                                                                             >  1:1160178/1‑150 (MQ=255)
                                                           cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGCTTGCTCATTACAGCATGCCGCTACCCAGCCCGCCGCCGAAGCCAACgg                                                                                            >  1:986195/1‑150 (MQ=255)
                                                                  taaCAGATGGAACGGCAAGAGAGAGGTCACTATACAACGGAG‑GGGGGAAGGATGTAAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCTGGGc                                                                                     <  1:169971/149‑1 (MQ=255)
                                                                  gAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                                     <  2:763422/150‑1 (MQ=255)
                                                                   aaCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCa                                                                                    <  2:564566/150‑1 (MQ=255)
                                                                      aGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAaca                                                                                 >  2:41598/1‑150 (MQ=255)
                                                                       gATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAacac                                                                                >  2:506732/1‑150 (MQ=255)
                                                                         tGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                                              <  1:1175719/150‑1 (MQ=255)
                                                                          ggAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                                                                             >  2:952166/1‑150 (MQ=255)
                                                                          ggAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                                                                             >  2:1224996/1‑150 (MQ=255)
                                                                          ggAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                                                                             >  1:801222/1‑150 (MQ=255)
                                                                           gAACCGCAAGAGTGAGGTCACTCTACACCGGAC‑GGGTTAAGGATGTTAGTCGGCGGGAATGAAAGCGGGGCGTGAACGCCCCGAGATTGATCATTACATCATTAAGCTACCCAGCACTTAGCCGAAGCAAACGGCCAGGGCAACACCCag                                                                            <  1:560329/150‑1 (MQ=255)
                                                                             aCCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                                                                          <  2:949044/150‑1 (MQ=255)
                                                                              ccGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAATACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGt                                                                         >  1:1738389/1‑150 (MQ=255)
                                                                              ccGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATCTTAGTCGGCGGGAAAGAAAGCGGCGCGTGACCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCCCGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGt                                                                         >  1:360869/1‑150 (MQ=255)
                                                                               cGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTa                                                                        >  1:903137/1‑150 (MQ=255)
                                                                               cGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTa                                                                        <  2:1257603/150‑1 (MQ=255)
                                                                                gCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGACGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAc                                                                       >  1:975587/1‑150 (MQ=255)
                                                                                  aaGAGAGAGGTCACTATACAACGGAC‑GGGGGTAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTAccc                                                                     >  1:1450121/1‑150 (MQ=255)
                                                                                     agagagGTCACTATACAACGGAC‑GGGGGAAGGATGTAAGTCGGCGGGACAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACATCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCggg                                                                  <  2:1652387/150‑1 (MQ=255)
                                                                                            tCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCCATAGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaaga                                                           >  1:1518461/1‑150 (MQ=255)
                                                                                             cACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaagaa                                                          >  2:1578896/1‑150 (MQ=255)
                                                                                              aCTATCCAACGGAC‑GGGGGAAGGAGGTTAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAAc                                                         <  2:1434006/150‑1 (MQ=255)
                                                                                               cTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg                                                        >  1:449154/1‑150 (MQ=255)
                                                                                                tataCAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACgg                                                       <  2:425497/150‑1 (MQ=255)
                                                                                                 ataCAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGa                                                      <  1:1235459/150‑1 (MQ=255)
                                                                                                    cAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCACCACCCAGAGTACCCGGGATGAAGAACGGATgg                                                   >  2:123987/1‑150 (MQ=255)
                                                                                                    cAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATgg                                                   >  2:1361503/1‑150 (MQ=255)
                                                                                                     aaCGGAC‑GGGGGAAGGATGATAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGt                                                  <  1:143551/150‑1 (MQ=255)
                                                                                                      acggac‑ggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCACGGCAACACCCAGAGTACCCTGCATCAAGCACGGATGGtt                                                 >  2:597762/1‑150 (MQ=255)
                                                                                                       cggac‑ggGGGATGGATGTTTGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTTCAGCATGTATCTACCCAGCACGAAGGCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                                <  1:1564480/150‑1 (MQ=255)
                                                                                                       cggac‑ggGGGAAGGATGTTAGTCGGCGGGAAAGACAGCGGTGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGAGGGTTg                                                <  1:1692006/150‑1 (MQ=255)
                                                                                                       cggac‑ggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                                <  1:190314/150‑1 (MQ=255)
                                                                                                          ac‑ggGGGAAGTATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAg                                             <  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)
                                                                                                                                                      |                                                                                                                                                     
GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGAC‑GGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGTA  >  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: ‑