Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,452 G→T 6.2% intergenic (‑96/+20) 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,4520GT 6.2% 400.1 / 10.9 137intergenic (‑96/+20)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base G (54/74);  new base T (3/6);  total (57/80)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

Only 100 of 137 total aligned reads displayed.

GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTG  >  NC_000913/4365304‑4365600
                                                                                                                                                     |                                                                                                                                                    
 ccGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACCACGGACGGGGGAAGGATGGTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                    <  2:1539087/150‑1 (MQ=255)
 ccGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGTGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                                    <  2:692082/150‑1 (MQ=255)
   gTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGTGGCGTGAACg                                                                                                                                                  <  2:1403381/150‑1 (MQ=255)
    tACATAGCACCACGACAGACGCGGGTTTCTAACTTTTTTCATCAAGCCTAAACACCTCGCAAGAACCGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGc                                                                                                                                                 <  2:923293/150‑1 (MQ=255)
      cATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGccc                                                                                                                                               >  1:819992/1‑150 (MQ=255)
         aGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCAAAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgc                                                                                                                                            <  1:1125049/150‑1 (MQ=255)
            accacGACAGACGCTGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                                         <  2:1305804/150‑4 (MQ=255)
            accacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgcttt                                                                                                                                         <  1:1210969/150‑4 (MQ=255)
             ccacGACAGACGCGGTATTTGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                        <  2:1246953/150‑1 (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)
              cacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAATAGAGAGGTCACTATACAACTGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTg                                                                                                                                       <  2:1351563/150‑1 (MQ=255)
              cacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGAGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTg                                                                                                                                       >  1:808673/1‑150 (MQ=255)
               acGACAGACGCTGTATTCGAACTTTTTTCATGAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttct                                                                                                                                       <  1:1018560/150‑6 (MQ=255)
               acGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGa                                                                                                                                      <  1:872414/150‑1 (MQ=255)
                cGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGAGGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGAt                                                                                                                                     <  1:615486/150‑1 (MQ=255)
                cGACAGACGCGGTATTCGAACATTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATCCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttcttt                                                                                                                                     <  1:281535/150‑8 (MQ=255)
                 gACAGACGCGGTATTCGAACTTTTTTCAGCAAGCATAACCACCTCGCAAGAACAGATGGAACCGAAAGAGATAGGGCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATc                                                                                                                                    <  2:1416064/150‑1 (MQ=255)
                  aCAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGACCCGCAAGCGAGAGGCCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCa                                                                                                                                   <  1:1745756/150‑1 (MQ=255)
                  aCAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCa                                                                                                                                   <  2:415890/150‑1 (MQ=255)
                   cactcgcgGTATTAGAACCTTTTTCAACAACCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                                                                  <  1:1193143/146‑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)
                       cgcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACa                                                                                                                              >  1:414406/1‑150 (MQ=255)
                        gcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                             <  1:1035401/150‑1 (MQ=255)
                         cgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGc                                                                                                                            <  2:758065/150‑1 (MQ=255)
                             aTTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGa                                                                                                                        >  1:210930/1‑150 (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:228356/150‑1 (MQ=255)
                                cGAAATTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                     <  1:634315/150‑1 (MQ=255)
                                      tttttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCa                                                                                                               <  1:793195/150‑1 (MQ=255)
                                        tttCATCAAGCATAAACACCTCGCAAGAACAGATGTAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                                                                             <  1:888103/150‑1 (MQ=255)
                                            atcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcga                                                                                                         >  2:1714790/1‑150 (MQ=255)
                                             tcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAATCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaa                                                                                                        >  2:1800280/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)
                                                      aaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCa                                                                                               >  1:1192798/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)
                                                                  gAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                                   <  2:763422/150‑1 (MQ=255)
                                                                   aaCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCa                                                                                  <  2:564566/150‑1 (MQ=255)
                                                                     cAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAac                                                                                >  2:728785/1‑150 (MQ=255)
                                                                      aGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAaca                                                                               >  2:41598/1‑150 (MQ=255)
                                                                       gATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAacac                                                                              >  2:506732/1‑150 (MQ=255)
                                                                        taggAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAcc                                                                             <  1:1599342/148‑1 (MQ=255)
                                                                        aTGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGTAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAcc                                                                             <  1:1392190/150‑1 (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)
                                                                             aCCGCAAGAGTGAGGTCACTATACAACGGACGGGGGAATGATTTTAGTCGGCGGGAAAGAAAGCGTGGCGTGAACGCCCCGCGATTGATCATTACTGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                                                                        <  2:1451524/150‑1 (MQ=255)
                                                                             aCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                                                                        <  2:949044/150‑1 (MQ=255)
                                                                             aCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGGGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                                                                        <  1:1498469/150‑1 (MQ=255)
                                                                              ccGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGCGt                                                                       >  1:1685275/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)
                                                                                       atcgatCACTATACAACGGACGGGGGAAGGATCTTAGGCGGCGGGAAATAAAGCGGGGCGTGACCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCTAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAt                                                              <  1:456802/145‑1 (MQ=255)
                                                                                               cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg                                                      >  1:449154/1‑150 (MQ=255)
                                                                                                 ataCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGa                                                    >  1:400484/1‑150 (MQ=255)
                                                                                                 ataCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGa                                                    <  1:1235459/150‑1 (MQ=255)
                                                                                                    cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATgg                                                 >  2:1361503/1‑150 (MQ=255)
                                                                                                     aaCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGt                                                <  2:1691271/150‑1 (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)
                                                                                                               gggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGa                                       >  2:170303/1‑150 (MQ=255)
                                                                                                                 gAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGCACGGATGGTTGAAGACGAAt                                     >  1:99461/1‑150 (MQ=255)
                                                                                                                      aTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAcc                                <  2:1417078/150‑1 (MQ=255)
                                                                                                                         ttAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGAt                             >  1:1723256/1‑150 (MQ=255)
                                                                                                                            gTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGACGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACg                          <  1:802311/150‑1 (MQ=255)
                                                                                                                               ggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgt                       <  2:819992/150‑1 (MQ=255)
                                                                                                                                  gggAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgtagt                    <  2:1192798/150‑1 (MQ=255)
                                                                                                                                       agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg               >  2:662611/1‑150 (MQ=255)
                                                                                                                                       agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg               >  2:254164/1‑150 (MQ=255)
                                                                                                                                        gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt              <  1:606337/150‑1 (MQ=255)
                                                                                                                                        gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt              >  2:1564368/1‑150 (MQ=255)
                                                                                                                                        gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCAATACGGGTAGTACCCgt              >  2:800995/1‑150 (MQ=255)
                                                                                                                                         atcgCGGGGCGTTCACGCCCCGCGATTGATCATTTCAGCATGAAGCTACCCAGCCCGAAGCCGCAGCCAACGGCCAGGGCACCACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTTCCGATACGGGTAGTACCCgtg             >  1:775528/4‑150 (MQ=255)
                                                                                                                                         atcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGTTGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg             >  2:703653/4‑150 (MQ=255)
                                                                                                                                         aaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg             <  2:449154/150‑1 (MQ=255)
                                                                                                                                         aaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg             <  1:950209/150‑1 (MQ=255)
                                                                                                                                             cGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGATGACGAATTTACCGATACGGGTAGTACCCGTGTcat         >  2:162864/1‑150 (MQ=255)
                                                                                                                                                 gCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatccat     <  2:273722/150‑1 (MQ=255)
                                                                                                                                                  cGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatccatc    >  1:1760004/1‑150 (MQ=255)
                                                                                                                                                    tGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTg  <  2:1234969/150‑1 (MQ=255)
                                                                                                                                                    tGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTg  <  1:920895/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTG  >  NC_000913/4365304‑4365600

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: