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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,4530AC 7.3% 400.0 / 15.3 142intergenic (‑97/+19)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (68/63);  new base C (6/5);  total (74/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 144 total aligned reads displayed.

GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGT  >  NC_000913/4365304‑4365601
                                                                                                                                                     |                                                                                                                                                    
   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)
           caccacGACAGACGCGGTATTGGAACTTTTTTCATCTACCATAAACACCTCGCAAGAACAGATGGAACCTCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctt                                                                                                                                           <  2:1451493/150‑3 (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:1755450/150‑1 (MQ=255)
             ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                         <  1:871254/150‑1 (MQ=255)
             ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                         >  2:1537418/1‑150 (MQ=255)
             ccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATAGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                         >  2:1579430/1‑150 (MQ=255)
              cacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACATATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTg                                                                                                                                        <  2:1635205/150‑1 (MQ=255)
                 gACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATc                                                                                                                                     <  2:330791/150‑1 (MQ=255)
                    aGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGAGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                                  <  1:41883/150‑1 (MQ=255)
                    aGACGCGGTATTCGAACTTTTATCATCAAGCATAAACACCTCGCAAGAACTGATGGAACCGCAATTGAGAGGTCACTCTACAACGGACGGGGGACGCATGTTAGACGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAtt                                                                                                                                  <  2:1368960/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)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACAGGGGAAGGATGTTAGTCGGCGGGAAGGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                      >  2:1665723/1‑150 (MQ=255)
                                 gAACTTTTATCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCt                                                                                                                     <  1:520910/150‑1 (MQ=255)
                                    cTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAAACGAAAGATAGAGGTCACTATAAAACGGACGGGGTAAGGATGTTACTCGGCGGGAAAGAAAGCGGGGCTTGAACGCCCCGCGATTCATCATAACAGCATGAAGCTAcc                                                                                                                  >  1:320864/1‑150 (MQ=255)
                                           catcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTTACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                                                                           >  2:743341/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                                                                                                        >  2:65838/1‑150 (MQ=255)
                                                    tatacacCTCGCAAGAACAGTTGGAACCGCAAGAGTGAGGTCACTATACAACGGACGGGGGATGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCGCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                  <  1:621530/147‑1 (MQ=255)
                                                    tAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                                  >  1:348851/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)
                                                          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)
                                                             cGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGcc                                                                                         >  2:1480827/1‑150 (MQ=255)
                                                             cGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGcc                                                                                         >  2:47072/1‑150 (MQ=255)
                                                               cAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAg                                                                                       <  2:1846518/150‑1 (MQ=255)
                                                                  gAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                                    <  1:1731110/150‑1 (MQ=255)
                                                                  caaCAGATGGAACCGCAAGAGAGAGGTCACTATACACCGGACGGGGGAAGGATGTTAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                                    <  1:1693074/149‑1 (MQ=255)
                                                                   aaCAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCa                                                                                   <  2:362317/150‑1 (MQ=255)
                                                                       gATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAacac                                                                               >  2:1734479/1‑150 (MQ=255)
                                                                         tGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                                             >  1:1843635/1‑150 (MQ=255)
                                                                           gAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                                                                           <  1:38761/150‑1 (MQ=255)
                                                                           gAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                                                                           >  2:976796/1‑150 (MQ=255)
                                                                            aaCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCaga                                                                          <  2:696031/150‑1 (MQ=255)
                                                                             aCCGCAAGAGAGAGGTCACTATACAACGGACGCGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                                                                         >  1:685530/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)
                                                                                       agagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAt                                                               <  1:1784892/150‑1 (MQ=255)
                                                                                       agagGTCACTATACAACCGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAt                                                               <  2:351641/150‑1 (MQ=255)
                                                                                         agGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATga                                                             <  2:53325/150‑1 (MQ=255)
                                                                                            tCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaaga                                                          >  2:8422/1‑150 (MQ=255)
                                                                                                   aCAACGGACTGGGGTTGGTTGTTATTCGGGGGGAATGAAAGGGGGGCTTTAACGCCACGAGATTGATCTTTTCAGCATGAAGTTACCCAGCACGAAGCAGAAGCAAACGGGCAGGGCAACACCCAGAGTACCCGGGATGATGAACGGATg                                                   <  1:1588007/150‑1 (MQ=255)
                                                                                                   aCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATg                                                   >  1:1326814/1‑150 (MQ=255)
                                                                                                    cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATgg                                                  >  1:1656142/1‑150 (MQ=255)
                                                                                                       cggacggGGGAAGGATGTTAGTCGGCGGGAAAGGAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                               <  2:74633/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)
                                                                                                             ggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACg                                         >  2:1056370/1‑150 (MQ=255)
                                                                                                                   ggATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTa                                   >  2:1656578/1‑150 (MQ=255)
                                                                                                                    gATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                                  >  2:982917/1‑150 (MQ=255)
                                                                                                                     aggTTAGGCGGCGCGAAAGCACGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTGGAAGACGAATTTAcc                                 <  1:588297/148‑1 (MQ=255)
                                                                                                                        ttAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGAt                              >  2:1627076/1‑150 (MQ=255)
                                                                                                                           gTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACg                           >  2:1744011/1‑150 (MQ=255)
                                                                                                                           gTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGAGGAAGAACGGATGGTTGACGGCGAATTTACAGATACg                           >  2:116829/1‑150 (MQ=255)
                                                                                                                                 gggAAAGAACGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCCGGGCAACACCCAGAGTACCCGGGATGAAGATCGGATGGTTGAAGACGAATTTACCGATACGGgtagt                     >  1:805631/1‑150 (MQ=255)
                                                                                                                                  ggAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTCCCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgtagta                    >  2:156343/1‑150 (MQ=255)
                                                                                                                                    aaagaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTAcc                  <  2:1596572/150‑1 (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)
                                                                                                                                       gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt               >  2:1143732/1‑150 (MQ=255)
                                                                                                                                        aaagCGGTTCTTGAACGAACCGCGATTGATCATTACAGCATGACGCTACCCATCACTAAGCCGAAGCAAACGGCCAGTGCAACACCCAGAGTACCCGGGATGAATAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg              <  2:483303/150‑1 (MQ=255)
                                                                                                                                        aaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg              <  1:152089/150‑1 (MQ=255)
                                                                                                                                         aagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAACCGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtgt             >  2:1472860/1‑150 (MQ=255)
                                                                                                                                          cgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTc            >  1:1655512/2‑150 (MQ=255)
                                                                                                                                          cgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTc            >  1:1476600/2‑150 (MQ=255)
                                                                                                                                          cgCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAAGGGCCAGGGCAACACCTAGTGTTCCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTc            <  1:955834/149‑1 (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:626871/1‑150 (MQ=255)
                                                                                                                                                    gAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGt  <  1:1589984/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGT  >  NC_000913/4365304‑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: ‑