Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,453 A→C 5.9% 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 5.9% 362.0 / 9.0 126intergenic (‑97/+19)cutA/dcuAdivalent‑cation tolerance protein, copper sensitivity/C4‑dicarboxylate antiporter
Reads supporting (aligned to +/- strand):  ref base A (61/57);  new base C (2/6);  total (63/63)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 128 total aligned reads displayed.

GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGTA  >  NC_000913/4365304‑4365602
                                                                                                                                                     |                                                                                                                                                     
gCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGCTGGAACCGCAAGAGAGAGGTCCCTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGa                                                                                                                                                       >  1:6616/1‑150 (MQ=255)
  cggaCATAGCACCACGACAGCCGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGTGAGAGGTCACTATACATCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAAc                                                                                                                                                     <  2:228330/147‑1 (MQ=255)
  cGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAAc                                                                                                                                                     >  1:897017/1‑150 (MQ=255)
  cGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAAc                                                                                                                                                     >  1:1897917/1‑150 (MQ=255)
  cGTACAAATCACCACTACAGACGCGGTATTCGACCTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGACGGAACCGAAAGATAGAGGTCACTATACACCGGACGGGGGAAGCATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAAc                                                                                                                                                     <  2:1025153/150‑1 (MQ=255)
         accaccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAATAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgc                                                                                                                                              <  2:1562107/148‑1 (MQ=255)
         aGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgc                                                                                                                                              >  1:16896/1‑150 (MQ=255)
          gCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgcg                                                                                                                                             >  2:1176309/1‑150 (MQ=255)
          gCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCcgcg                                                                                                                                             <  1:1218391/150‑1 (MQ=255)
           caccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAAGAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGCACGCCCCGCGa                                                                                                                                            >  2:100925/1‑150 (MQ=255)
            accacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGATAGAAAGCGGGGCGTGAACGCCCCGCGAt                                                                                                                                           >  1:1039669/1‑150 (MQ=255)
             ccacGACAGCCGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                          <  2:1656636/150‑1 (MQ=255)
              cacGCCCGACGCTGGATTCGAAATTTTTTCATCAAGACTAACCACATCGCAAGACAAGAGGTACCCGCAAGAGAGAGGTCACTAAACACCGGCCGGGGGAAGGATGTTACTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTg                                                                                                                                         <  2:935877/150‑1 (MQ=255)
               acGACATACGCGGTATTCGAACTTATATCATCAAGCAAAAACACCTCGCAAGAACAGATGGAACCCCAAGAGAGAGGGCACTAAACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttctt                                                                                                                                        <  1:525500/150‑7 (MQ=255)
               acGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcgctttctt                                                                                                                                        <  1:1256077/150‑7 (MQ=255)
                   cAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGGCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAt                                                                                                                                    <  2:1644014/150‑1 (MQ=255)
                   cAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGCAACCGCAAGAGAGAGGTCACTATAAAACGGACGGGGGAAGTAAGTTAGTCTGCGGGAAAGAAAGAGGGGCGTGAACGCCCCGACATTGATCAt                                                                                                                                    >  2:1205125/1‑150 (MQ=255)
                      ccgcgGTATTCGAACTTTATTCATCACGCATAACCACCTCGCAAGAACAGATGGCCCCGCAAGAGAGCGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTAc                                                                                                                                 <  1:193554/149‑1 (MQ=255)
                      aCGCGGTATTCGAACTTTTTTCATCACGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTAc                                                                                                                                 <  1:373833/150‑1 (MQ=255)
                        gcgGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                               <  1:475905/150‑1 (MQ=255)
                        gcgGTATTCGAACTTTTATCATCAAGCATAAACCCCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                               <  1:1335502/150‑1 (MQ=255)
                          ggTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCa                                                                                                                             >  2:1440446/1‑150 (MQ=255)
                             aTTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGa                                                                                                                          >  1:575770/1‑150 (MQ=255)
                             aTTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGa                                                                                                                          <  2:248377/150‑1 (MQ=255)
                               tCGAACTTTTTTCATCACGCATACCCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACAATACACCGGACGGGGGAAGGATGTTAGTAGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                                                                        <  2:408916/150‑1 (MQ=255)
                               tCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAg                                                                                                                        <  1:1072209/150‑1 (MQ=255)
                                cGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                       >  2:1795294/1‑150 (MQ=255)
                                   aattttttCATCAAGCATAAACCCCTCGCCAGAACAGATGGAACCGCAAGAGAGAGGGCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTAc                                                                                                                    <  1:347164/148‑1 (MQ=255)
                                      tttttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCa                                                                                                                 <  2:1327645/150‑1 (MQ=255)
                                        tttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGc                                                                                                               >  1:1513440/1‑150 (MQ=255)
                                        tttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAAAGGACGGGGGAAGGATGTTAGTCGGCGCGTAAGAAATCGTGGCGTGAACGCCCAGCGATTGATCATTAGAGCATGAAGCTACCCAGc                                                                                                               >  2:294871/1‑150 (MQ=255)
                                          tcatcaAGCATAAACACCTCGCAAGAACAGCTGGAACCGCAAGAGAGAGCTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAc                                                                                                             <  1:54817/150‑1 (MQ=255)
                                          tcatcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAc                                                                                                             <  2:1111493/150‑1 (MQ=255)
                                           catcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcg                                                                                                            >  2:1338848/1‑150 (MQ=255)
                                            atcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATTCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACACAGCAaga                                                                                                           >  2:909117/1‑150 (MQ=255)
                                            atcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcga                                                                                                           >  2:1290560/1‑150 (MQ=255)
                                               aaccATAAACACCTCGCAAGATCAGATGGAAGCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagc                                                                                                        <  1:500765/147‑1 (MQ=255)
                                               aaGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCATCAcgaagc                                                                                                        <  1:606646/150‑1 (MQ=255)
                                                aGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagcc                                                                                                       >  2:1210774/1‑150 (MQ=255)
                                                aGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAA‑GATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccg                                                                                                      >  1:903704/1‑150 (MQ=255)
                                                     aaaCACCTCGCAAGAACAGAGGGAACCGCAAGAGAGAGGTCACTATCCAACGGACGGGGGAAGGAGGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaagc                                                                                                  <  1:1392070/150‑1 (MQ=255)
                                                        cacCTCGCACGAACAGATGTAACCGCAAGAGATAGTTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCaaa                                                                                               <  1:1733552/150‑1 (MQ=255)
                                                           cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACgg                                                                                            >  1:958013/1‑150 (MQ=255)
                                                             cGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGcc                                                                                          >  2:1086812/1‑150 (MQ=255)
                                                             cGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGcc                                                                                          >  1:604274/1‑150 (MQ=255)
                                                                aaGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgg                                                                                       >  2:187226/1‑150 (MQ=255)
                                                                  gAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGCCGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGc                                                                                     <  1:830658/150‑1 (MQ=255)
                                                                     cAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAac                                                                                  >  1:549227/1‑150 (MQ=255)
                                                                         tGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                                              >  2:484744/1‑150 (MQ=255)
                                                                         tGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACAccc                                                                              >  2:1833241/1‑150 (MQ=255)
                                                                           gAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                                                                            <  1:1192322/150‑1 (MQ=255)
                                                                           cccccGCACCAGAGCGGGCACTATACAACGGACGGGGGAAGGCTGTTAGTCGGCGGGACAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCag                                                                            <  2:1618693/147‑1 (MQ=255)
                                                                             aCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                                                                          <  2:1543850/150‑1 (MQ=255)
                                                                             aCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                                                                          <  1:1176309/150‑1 (MQ=255)
                                                                              ccGCAAGAGAGAGGTCACTATACACCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCAATACCGCATGCAGCTACCCAGCACGCAGCCGAAGCACACGGCCAGGGCAACACCCAGAGt                                                                         >  1:1656355/1‑150 (MQ=255)
                                                                                   agagagagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCg                                                                    <  1:659617/150‑1 (MQ=255)
                                                                                          ggacACTATACAACGGACGGGGTAAGCATGTATGTCGGCGGGAAAGAAAGCGGGGCGTGACCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaa                                                             <  2:1841144/147‑1 (MQ=255)
                                                                                          ggTCAATATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGATCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaa                                                             <  2:1044475/150‑1 (MQ=255)
                                                                                             cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaagaa                                                          >  1:1078046/1‑150 (MQ=255)
                                                                                             cACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaagaa                                                          >  2:1190574/1‑150 (MQ=255)
                                                                                              aCTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGCCCAGGGAAACACCCAGAGTACACGGGATTAAGAAc                                                         >  1:53410/1‑150 (MQ=255)
                                                                                               cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg                                                        <  1:1882639/150‑1 (MQ=255)
                                                                                               cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg                                                        >  2:1186399/1‑150 (MQ=255)
                                                                                               cTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg                                                        <  1:1440446/150‑1 (MQ=255)
                                                                                                    cAACGGACGGGGGATGTTTGTTTTGCGGCGGTAAAGATAGCGGGGCGTGAACGCCACGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATgg                                                   <  1:1404034/150‑1 (MQ=255)
                                                                                                    cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACCGCATGAAGATACCAAGCCCGATGCCGCAGCAACCGGCCAGGGCAACACCCAGAGTACCCGGGATGACGCCCGGATgg                                                   >  2:1192336/1‑150 (MQ=255)
                                                                                                     aaCGGCCGGGGGAAGGATGTTAGGCGGCGGGAAAGAAAGCGGGGCGGGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGt                                                  <  2:441929/150‑1 (MQ=255)
                                                                                                     aaCGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGt                                                  >  1:1126730/1‑150 (MQ=255)
                                                                                                       cggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                                <  2:363352/150‑1 (MQ=255)
                                                                                                       cggacggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                                >  2:836196/1‑150 (MQ=255)
                                                                                                             ggggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACg                                          >  2:511745/1‑150 (MQ=255)
                                                                                                               ggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCTATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGaa                                        <  1:994314/150‑1 (MQ=255)
                                                                                                               ggAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGaa                                        <  1:1086812/150‑1 (MQ=255)
                                                                                                                gAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCTTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAAt                                       <  1:562819/150‑1 (MQ=255)
                                                                                                                gAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAAt                                       >  2:582372/1‑150 (MQ=255)
                                                                                                                 aaGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAAtt                                      >  1:1039466/1‑150 (MQ=255)
                                                                                                                   gtagcttAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCACGCGATTGATCATTAAATCATGAAGCTACCCAGCACGAAGCCGAAGCATACGTCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGTATGGTTGAAGACGAATTTa                                    <  1:294871/145‑1 (MQ=255)
                                                                                                                   ggATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTa                                    >  2:1132872/1‑150 (MQ=255)
                                                                                                                    gATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                                   >  1:1030051/1‑150 (MQ=255)
                                                                                                                    gATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                                   >  1:1614814/1‑150 (MQ=255)
                                                                                                                     aTGTTAGTCGGCGGGAAAGAAAGCGGGGCGGGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAcc                                  <  2:1736211/150‑1 (MQ=255)
                                                                                                                       gTTAGTCGGCGGGAAAGAAAGCGGGGCGTGACCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGa                                <  2:733609/150‑1 (MQ=255)
                                                                                                                             cggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCCGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACggg                          <  1:1115996/150‑1 (MQ=255)
                                                                                                                             cggcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACggg                          >  2:918323/1‑150 (MQ=255)
                                                                                                                               gcggGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgta                        <  2:134560/150‑1 (MQ=255)
                                                                                                                                     aagaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTAccc                  >  1:1796294/1‑150 (MQ=255)
                                                                                                                                       gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt                >  2:1544216/1‑150 (MQ=255)
                                                                                                                                       gaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt                >  1:1325453/1‑150 (MQ=255)
                                                                                                                                        aaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg               <  1:1098239/150‑1 (MQ=255)
                                                                                                                                         aagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGACGACGAATTTACCGATACGGGTAGTACCCgtgt              >  2:1763840/1‑150 (MQ=255)
                                                                                                                                         aagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtgt              <  1:936358/150‑1 (MQ=255)
                                                                                                                                           gCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGCTACGGGTAGTACCCGTGTca            >  2:1850312/1‑150 (MQ=255)
                                                                                                                                            cGGGGCTTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcat           <  1:1107165/150‑1 (MQ=255)
                                                                                                                                            cGGGGCGTGAACGCCCCGGGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGGTGAAGACGAATTTACCGATACGGGTACTACCCGTATcgg           <  1:1609960/150‑3 (MQ=255)
                                                                                                                                               ggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatcca        >  1:526475/1‑150 (MQ=255)
                                                                                                                                               ggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGGcatcca        >  1:1039942/1‑150 (MQ=255)
                                                                                                                                                gCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTcatccat       >  2:1030154/1‑150 (MQ=255)
                                                                                                                                                   tGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTg    <  1:1833241/150‑1 (MQ=255)
                                                                                                                                                   tGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTg    >  1:1468948/1‑150 (MQ=255)
                                                                                                                                                     aaCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGTa  <  1:1190574/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGTACATAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCCATCTGTA  >  NC_000913/4365304‑4365602

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑