Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,452 G→T 6.0% 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.0% 392.4 / 10.9 142intergenic (‑96/+20)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base G (70/62);  new base T (3/6);  total (73/69)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 145 total aligned reads displayed.

TGCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGT  >  NC_000913/4365303‑4365601
                                                                                                                                                     |                                                                                                                                                     
 gCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGa                                                                                                                                                      <  1:1825436/150‑1 (MQ=255)
    gTACATAGCACCACTACAGACGCGGTATTCGAACTTTTATCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACg                                                                                                                                                   <  2:1629452/150‑1 (MQ=255)
       cATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGccc                                                                                                                                                >  1:834390/1‑150 (MQ=255)
        aTAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGcccc                                                                                                                                               <  1:256555/150‑1 (MQ=255)
         tAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACAATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcg                                                                                                                                              <  2:1050588/150‑1 (MQ=255)
          accaccacGACAGACGCTGTATTCGAACTTTTTTCATCAAGTATAAACACGTCGCAAGAACAGATGGACCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAACGCGGGGCGTGAACGCCCcgc                                                                                                                                             <  1:923385/148‑1 (MQ=255)
           gCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAGCACCGCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgcg                                                                                                                                            >  2:156159/1‑150 (MQ=255)
           gCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgcg                                                                                                                                            >  1:1776793/1‑150 (MQ=255)
            caccacGACAGACGCGGTATTCGAACTTTTTTCATCAACCATAAACACCTAGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGa                                                                                                                                           >  2:1588007/1‑150 (MQ=255)
             accacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGCGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGACCGCCCCGCGag                                                                                                                                          >  2:1836084/1‑149 (MQ=255)
              ccacGACTGACGCGGTATTCCAACTTTTTTCATCAATCAAAAACCCCTCCCACCAACAGAGGGAACCGCAAGATATAGCCCACTATACAACGGACGGGGGAAGTATGTTAGTCGGCGGGAAACAAAGCGGGGCGTCAACCCCCCGCGAtt                                                                                                                                         <  2:1868861/150‑1 (MQ=255)
              ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCTTGAACGCCCCGCGAtt                                                                                                                                         >  2:894519/1‑150 (MQ=255)
              ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttc                                                                                                                                         <  2:613963/150‑5 (MQ=255)
              ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                         >  2:1537418/1‑150 (MQ=255)
              ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATAGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                         >  2:1579430/1‑150 (MQ=255)
               caggtcAGACGCGGTATTCGAACTTTTTTCATCACGCATAATCACGTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTg                                                                                                                                        <  2:722568/145‑1 (MQ=255)
                  gACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATc                                                                                                                                     <  2:330791/150‑1 (MQ=255)
                        cgcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACa                                                                                                                               >  2:1439892/1‑150 (MQ=255)
                         gcgGTATTCGAACTTTTATCATCAAGCATAAACCCCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                              <  1:1493529/150‑1 (MQ=255)
                              aTTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGAGGGAAAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCATAACGGCATGa                                                                                                                         >  1:580835/1‑150 (MQ=255)
                                tCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                                                                       <  1:1070183/150‑1 (MQ=255)
                                 cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                      >  1:1406253/1‑150 (MQ=255)
                                 cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                      >  1:589411/1‑150 (MQ=255)
                                 cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACAGGGGAAGGATGTTAGTCGGCGGGAAGGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                      >  2:1665723/1‑150 (MQ=255)
                                  gAACTTTTATCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCt                                                                                                                     <  1:520910/150‑1 (MQ=255)
                                     cTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAAACGAAAGATAGAGGTCACTATAAAACGGACGGGGTAAGGATGTTACTCGGCGGGAAAGAAAGCGGGGCTTGAACGCCCCGCGATTCATCATAACAGCATGAAGCTAcc                                                                                                                  >  1:320864/1‑150 (MQ=255)
                                      ttttttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAAAGCCCCGCGATTGATCATTACAGCATGAAGCTAccc                                                                                                                 >  1:974339/1‑150 (MQ=255)
                                             atcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGCAGCTACCCAGCAcga                                                                                                          >  1:639295/1‑150 (MQ=255)
                                             atcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcga                                                                                                          >  1:754085/1‑150 (MQ=255)
                                               caAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTTCAGCATGAAGCTACCCAGCAcgcag                                                                                                        >  2:1176701/1‑150 (MQ=255)
                                               caAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaag                                                                                                        >  1:440392/1‑150 (MQ=255)
                                               caAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaag                                                                                                        >  2:65838/1‑150 (MQ=255)
                                               caAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTAAACCCCGGACGGGGGAAGGATGTACGTCGGCGGGAAAGAACGCGGGGGCTGCACGCCCCGCGATTGCTCATTCGAGCCTGAAGCTTCCCAGCAcgaag                                                                                                        >  2:1703895/1‑150 (MQ=255)
                                                     tatacacCTCGCAAGAACAGTTGGAACCGCAAGAGTGAGGTCACTATACAACGGACGGGGGATGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCGCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                  <  1:621530/147‑1 (MQ=255)
                                                     tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                  >  1:628826/1‑150 (MQ=255)
                                                     tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGACGTCACTATACAACCGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                  <  1:838606/150‑1 (MQ=255)
                                                       aaCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACCGCATGAAGCTACCCAGCACGAAGCCGAAGCa                                                                                                >  2:320810/1‑150 (MQ=255)
                                                        acacCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgccgcac                                                                                               >  1:1711627/1‑144 (MQ=255)
                                                          acCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAAc                                                                                             <  2:754085/150‑1 (MQ=255)
                                                           ccTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGTATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCACCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACg                                                                                            <  2:1749801/150‑1 (MQ=255)
                                                           ccTCGCAAGAACAGATGCAACCGCAAGAGAGAGGTCAATAAACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAACGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACg                                                                                            <  2:1667096/150‑1 (MQ=255)
                                                            cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGTGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACgg                                                                                           <  1:105538/150‑1 (MQ=255)
                                                            cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACgg                                                                                           <  1:101003/150‑1 (MQ=255)
                                                            cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACgg                                                                                           >  1:474940/1‑150 (MQ=255)
                                                             tCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGc                                                                                          >  2:862858/1‑150 (MQ=255)
                                                                   caaCAGATGGAACCGCAAGAGAGAGGTCACTATACACCGGACGGGGGAAGGATGTTAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                                    <  1:1693074/149‑1 (MQ=255)
                                                                    aaCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCa                                                                                   <  2:362317/150‑1 (MQ=255)
                                                                          tGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                                             >  1:1843635/1‑150 (MQ=255)
                                                                            gAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                                                                           >  2:976796/1‑150 (MQ=255)
                                                                            gAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                                                                           <  1:38761/150‑1 (MQ=255)
                                                                             aaCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCaga                                                                          <  2:696031/150‑1 (MQ=255)
                                                                             aaCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCaga                                                                          >  2:1459993/1‑150 (MQ=255)
                                                                               ccGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCAGGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGt                                                                        >  1:349214/1‑150 (MQ=255)
                                                                                cGCAAGAGAGAGGTCACTATACAACGGACGGGAGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACCCCCAGAGTa                                                                       >  2:654107/1‑150 (MQ=255)
                                                                                       gagagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGa                                                                >  1:1606120/1‑150 (MQ=255)
                                                                                       gagagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGCGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGa                                                                <  1:249541/150‑1 (MQ=255)
                                                                                        agagGTCACTATACAACCGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAt                                                               <  2:351641/150‑1 (MQ=255)
                                                                                         gcgggcACTATACAACGGACGGGGGAAGGATGTTAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATg                                                              <  2:1381010/145‑1 (MQ=255)
                                                                                          agGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATga                                                             <  2:53325/150‑1 (MQ=255)
                                                                                             tCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaaga                                                          >  2:8422/1‑150 (MQ=255)
                                                                                                cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg                                                       <  1:818490/150‑1 (MQ=255)
                                                                                                   aaCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGACCCGCGATTGATCATTACAGCATGAAGCTACACAGCACCAAGCCGAAGCAAACGGCCAGGGCAACAACCAGAGTACCCGGGATGAAGAACGGAt                                                    >  1:1744152/2‑150 (MQ=255)
                                                                                                    aCAACGGACTGGGGTTGGTTGTTATTCGGGGGGAATGAAAGGGGGGCTTTAACGCCACGAGATTGATCTTTTCAGCATGAAGTTACCCAGCACGAAGCAGAAGCAAACGGGCAGGGCAACACCCAGAGTACCCGGGATGATGAACGGATg                                                   <  1:1588007/150‑1 (MQ=255)
                                                                                                    aCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATg                                                   >  1:1326814/1‑150 (MQ=255)
                                                                                                     cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATgg                                                  >  1:1656142/1‑150 (MQ=255)
                                                                                                      aaCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCTTTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGt                                                 <  2:1839816/150‑1 (MQ=255)
                                                                                                        cggacggGGGAAGGATGTTAGTCGGCGGGAAAGGAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                               <  2:74633/150‑1 (MQ=255)
                                                                                                         ggacggGGGAAGTATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGTTTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGa                                              <  1:1539689/150‑1 (MQ=255)
                                                                                                           acggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAg                                            <  2:1099205/150‑1 (MQ=255)
                                                                                                           acggGGGAAGGATGTTAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAg                                            <  1:698431/150‑1 (MQ=255)
                                                                                                            cggcggAAGGATGGGAGGCGGCGGGAACGAAAGCGGGCCGTGAACGCCCGGCGATTGAGCATTACAGCCTGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTCCCCGGGATGAAGAACGGATGGTTGAAGa                                           <  2:1143262/150‑1 (MQ=255)
                                                                                                             gggggAAGTATGTTAGTCGGCGGTAAAGAACGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGTAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGTATGAAGAACGGATGGTTGAAGAc                                          <  1:1457404/150‑1 (MQ=255)
                                                                                                             gggggAAGGATGTTAGTTGGCGGGTAAGACAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGTTTGAAGAc                                          <  1:1449946/150‑1 (MQ=255)
                                                                                                              ggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGACGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAGGAAGAACGGATGGTTGAAGACg                                         >  2:1154549/1‑150 (MQ=255)
                                                                                                                    ggATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTa                                   >  2:1656578/1‑150 (MQ=255)
                                                                                                                     gATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                                  >  2:982917/1‑150 (MQ=255)
                                                                                                                            gTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACg                           >  2:1744011/1‑150 (MQ=255)
                                                                                                                            gTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAGGAAGAACGGATGGTTGACGGCGAATTTACAGATACg                           >  2:116829/1‑150 (MQ=255)
                                                                                                                                   ggAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTCCCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgtagta                    >  2:156343/1‑150 (MQ=255)
                                                                                                                                       cgaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg                <  2:1500904/149‑1 (MQ=255)
                                                                                                                                       caatcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAACCGGCCAGGGCAACCCCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGCCGAATTTACCGATACGGGTAGTACCCg                >  1:917387/6‑150 (MQ=255)
                                                                                                                                       agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg                <  2:188218/150‑1 (MQ=255)
                                                                                                                                       agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCg                <  2:1615812/150‑1 (MQ=255)
                                                                                                                                        gaaagCGGGGCGTGCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCCACCCCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTTCCGATACGGGTAGTACCCgg               >  2:1657802/1‑149 (MQ=255)
                                                                                                                                        gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCACACGGCCAGGGCCACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt               >  1:337850/1‑150 (MQ=255)
                                                                                                                                         aaagCGGTTCTTGAACGAACCGCGATTGATCATTACAGCATGACGCTACCCATCACTAAGCCGAAGCAAACGGCCAGTGCAACACCCAGAGTACCCGGGATGAATAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg              <  2:483303/150‑1 (MQ=255)
                                                                                                                                          aagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAACCGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtgt             >  2:1472860/1‑150 (MQ=255)
                                                                                                                                          aagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAACCGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtgt             >  2:433471/1‑150 (MQ=255)
                                                                                                                                           cgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTc            >  1:1655512/2‑150 (MQ=255)
                                                                                                                                           cgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTc            >  1:1476600/2‑150 (MQ=255)
                                                                                                                                             cGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcat          >  2:1589984/1‑150 (MQ=255)
                                                                                                                                             cGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcat          <  2:1631114/150‑1 (MQ=255)
                                                                                                                                              ggggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatc         >  2:1624337/1‑150 (MQ=255)
                                                                                                                                               gggCGTGAACGCCCCGCGATTGATCATTACAGGATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatcc        <  2:685530/150‑1 (MQ=255)
                                                                                                                                               gggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatcc        >  2:594967/1‑150 (MQ=255)
                                                                                                                                                ggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatcca       >  1:543233/1‑150 (MQ=255)
                                                                                                                                                  cGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatccatc     <  1:1616885/150‑1 (MQ=255)
                                                                                                                                                   gTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCt    <  2:834390/150‑1 (MQ=255)
                                                                                                                                                    tGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTg   <  1:1836084/150‑1 (MQ=255)
                                                                                                                                                     gAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGt  <  1:1589984/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGT  >  NC_000913/4365303‑4365601

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: