Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,365,453 A→C 8.2% intergenic (‑97/+19) cutA ← / ← dcuA copper binding protein CutA/C4‑dicarboxylate transporter DcuA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,365,4530AC 8.2% 114.7 / 3.0 62intergenic (‑97/+19)cutA/dcuAcopper binding protein CutA/C4‑dicarboxylate transporter DcuA
Reads supporting (aligned to +/- strand):  ref base A (18/35);  new base C (3/3);  total (24/38)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

TAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCC  >  NC_000913/4365312‑4365595
                                                                                                                                             |                                                                                                                                              
tAGCACCACGACAGACGCGGTCTTCGACCTCTTTTCATCAAGCATCAGCACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                               <  2:84013/143‑1 (MQ=255)
tAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGaa                                                                                                                                               <  1:873995/143‑1 (MQ=255)
 agaaccacGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACGATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAAc                                                                                                                                              <  2:400193/140‑1 (MQ=255)
  gCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCGAGAACAGATGGAACAGCCAGAGAGAGGTCACTATACAACGGACGGGGGAAGGAGGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGGACg                                                                                                                                             >  2:277030/1‑143 (MQ=255)
   caccacGACAGACGCGGTATTCGAACTTTTTTCATGAAGCATAAACACCTCGAACGAACATGTGGAACCGCAAGAAATAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGGGGGAAAGACAGCGGGGCGGGAACGc                                                                                                                                            >  2:225767/1‑143 (MQ=255)
       acGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGGCACTATCCGACGGGCGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTTCACGCCCcg                                                                                                                                        <  2:1046857/143‑1 (MQ=255)
            aGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGAtt                                                                                                                                   <  1:821811/143‑1 (MQ=255)
             gACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGGGGGGCGGGGGCGCCCCGCGAGTg                                                                                                                                  >  2:1036477/1‑143 (MQ=255)
              aCGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAACGCGGGGCGTGAACGCCCcgcgcgtga                                                                                                                                 >  2:633295/1‑138 (MQ=255)
                    tATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTa                                                                                                                           <  1:240716/143‑1 (MQ=255)
                       tCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAg                                                                                                                        >  1:1032818/1‑143 (MQ=255)
                           aCTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCCAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGGCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATg                                                                                                                    <  1:999499/143‑1 (MQ=255)
                           aCTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATg                                                                                                                    <  1:528166/143‑1 (MQ=255)
                            cTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACACCGGACGGGGGAAGGATGTTGGTGCAACACATATAATTCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATgc                                                                                                                   >  1:730037/1‑142 (MQ=255)
                             ttttttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACAGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGaa                                                                                                                  >  2:949771/1‑143 (MQ=255)
                               ttttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGc                                                                                                                >  2:820294/1‑143 (MQ=255)
                               ttttCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGAGGTTAGTCGGCGGGAAAGAAAGCGGGGCGGGAACGCCCCGCGATGGATCAGTACACCGTGCAGc                                                                                                                >  2:861678/1‑143 (MQ=255)
                                   catcaAGCAAAAGCCCCTCGTCAGCCGATATCGACCCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGTGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTAcc                                                                                                            <  1:1009843/143‑1 (MQ=255)
                                     tcaAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGACACTATACAACGGACGGGGGAAGGGGGTTAGTCGGCGGGAAAGAAAGCGGGGCGGGAAGGCCGCGCGAGGGATCAGGAGAGCAGGAGGCGAccgc                                                                                                          >  2:503107/1‑141 (MQ=255)
                                                   cTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCAcgaagccgaag                                                                                            >  1:742097/1‑143 (MQ=255)
                                                        aaGAACAGATGGAAGCGCAAGAGAGAGGTCACTATGCAACGGACGGGGTAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAAc                                                                                       <  1:721552/143‑1 (MQ=255)
                                                            agaGATGTACCGGCAAGAGAGCGGTGTCTATACAACGGCCGGTGTGGTGATGTTGGTCGGCGGGAAAGAACGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGcc                                                                                   <  2:368395/141‑1 (MQ=255)
                                                               taTGGACCCGCAAGGGAGAGGTCACTAGACAGCGGGCGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAgg                                                                                <  2:403642/142‑1 (MQ=255)
                                                                         cAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAGGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTAACCAGCACGAAGGCGAACCCCACGGCCCGGGCCACACCCa                                                                      >  1:182701/1‑143 (MQ=255)
                                                                         cAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCa                                                                      <  1:622811/143‑1 (MQ=255)
                                                                            gagagagGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCagag                                                                   >  2:392213/1‑143 (MQ=255)
                                                                                     cagtataGAAGGGACGGGGGGGGGATGTTGGTCGGCGGGAAAGAAGGCGGGGCGTGAGCGCCCCGCGATGGATCATAACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGa                                                          <  2:705265/140‑1 (MQ=255)
                                                                                      gcTACACGACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCATAGTACCCGGGAt                                                         <  2:941715/142‑1 (MQ=255)
                                                                                      aCTATACAACGGACGGGGGAAGGCTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGGACGCCCCGCGAGTGAGCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGCGGACCCGGGac                                                         >  2:92873/1‑142 (MQ=255)
                                                                                         ataCAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaa                                                      <  1:616597/143‑1 (MQ=255)
                                                                                            cAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATgaagaa                                                   <  1:820294/143‑1 (MQ=255)
                                                                                             aaggggcggGGGAGGGAGGTGAGTGGGGGGGAAGGGAAGGGGGGCGGGAACGGCCCGCGAGTGATCATTAGAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAAc                                                  <  1:308664/137‑1 (MQ=255)
                                                                                             aaCGGACGGGGGAAGGCTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCAGCGATTGATGATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCGAACGGGCAGGGCAACACCCAGAGAACCAGGGATGAAGGAc                                                  >  1:99341/1‑143 (MQ=255)
                                                                                              acggagggGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACg                                                 <  2:720827/143‑1 (MQ=255)
                                                                                                      gagAAGGATGTTAGTCGGCGTGAAAGAATGCGGGGCGGGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTg                                         <  2:274341/141‑1 (MQ=255)
                                                                                                          aGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGa                                     >  1:177218/1‑143 (MQ=255)
                                                                                                          aGGATGTCAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGa                                     <  2:283399/143‑1 (MQ=255)
                                                                                                           ggagggtAGGCGGGGGGAAGGGAGGCGGGGCGTGAGCGCCCCGCGATTGATCGTTACAGCACGAAGCTGCCCAGCGCGAGGCCGAAGCAAACGGGCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGAc                                    <  2:485229/137‑1 (MQ=255)
                                                                                                             aTGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGa                                  >  2:827044/1‑143 (MQ=255)
                                                                                                               gTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAAt                                <  2:661106/143‑1 (MQ=255)
                                                                                                                ttAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAAtt                               <  1:921389/143‑1 (MQ=255)
                                                                                                                 tAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAAttt                              <  1:864603/143‑1 (MQ=255)
                                                                                                                 tAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGGACCCGGGATGAAGAACGGAGGGTTGACGACGAAttt                              >  2:1071975/1‑143 (MQ=255)
                                                                                                                   gtgggcggGAAAGAACGCGGGCCGTCAACGCGCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                            <  2:769306/140‑1 (MQ=255)
                                                                                                                   gTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAc                            <  1:392213/143‑1 (MQ=255)
                                                                                                                    tgggcggGACAGAAAGCGGGGCGTGAACGCCCCGCGATTGACCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAcc                           <  1:936429/141‑1 (MQ=255)
                                                                                                                    tcagcaggAAGGAGAGCGGGGAGTGAACGCCCCGGGATTGAACATGACAGCATGAAGCTACCCACCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTAcc                           <  1:387297/137‑1 (MQ=255)
                                                                                                                            aaagaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTCACCGCTACggg                   >  2:621895/1‑143 (MQ=255)
                                                                                                                             atgaaagCGGGGCATGACCGCCCCGAGATTTATCATTACAGCATGAAGCTACCTAGGACGAAGCCGACTCAAACGGCCAGGGTTACACCCGGAGTACCTGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgt                  <  2:85682/141‑1 (MQ=255)
                                                                                                                             acgaaagCGGGGCGTGAACGCGCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGCCGAATTTACCGATACGGgt                  <  2:590682/141‑1 (MQ=255)
                                                                                                                              agaaagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgta                 <  2:1032818/143‑1 (MQ=255)
                                                                                                                                atcgCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGCAGAACGGATGGTTGAAGACGAAGTTACCGACACGGgtcgc               >  2:365836/4‑140 (MQ=255)
                                                                                                                                 aagGGGGGCGTGAACGCCCCGCGAATAATCACTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgtagta              <  1:565596/143‑1 (MQ=255)
                                                                                                                                 aagCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGgtagta              <  1:1017964/143‑1 (MQ=255)
                                                                                                                                   gCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGCAGCACACGGCCAGGGCAACACCCAGAGTCAGCGGGAACAACAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTAcc            >  1:448622/1‑143 (MQ=255)
                                                                                                                                   gCGGGGCGTTCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTAcc            >  1:879735/1‑143 (MQ=255)
                                                                                                                                   gCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGCAGACCAATTTACCGATCCCGGGAGGAcc            >  2:504561/1‑143 (MQ=255)
                                                                                                                                      gggCGTGCACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt         <  1:166818/143‑1 (MQ=255)
                                                                                                                                      gggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt         <  1:827044/143‑1 (MQ=255)
                                                                                                                                      gggCGGGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgt         <  1:494456/143‑1 (MQ=255)
                                                                                                                                       ggCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAAAACCCAGAGAACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtg        >  1:850160/1‑143 (MQ=255)
                                                                                                                                        gCGTGAACGCCCCGCGATTGGTCATTACCGCATGAAGCTACCCCGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTACTACCCgtat       <  1:393540/143‑3 (MQ=255)
                                                                                                                                        gCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCgtgt       <  1:344833/143‑1 (MQ=255)
                                                                                                                                          gTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTca     <  1:762549/143‑1 (MQ=255)
                                                                                                                                             aaCGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGCcaccc  >  2:78870/1‑143 (MQ=255)
                                                                                                                                             |                                                                                                                                              
TAGCACCACGACAGACGCGGTATTCGAACTTTTTTCATCAAGCATAAACACCTCGCAAGAACAGATGGAACCGCAAGAGAGAGGTCACTATACAACGGACGGGGGAAGGATGTTAGTCGGCGGGAAAGAAAGCGGGGCGTGAACGCCCCGCGATTGATCATTACAGCATGAAGCTACCCAGCACGAAGCCGAAGCAAACGGCCAGGGCAACACCCAGAGTACCCGGGATGAAGAACGGATGGTTGAAGACGAATTTACCGATACGGGTAGTACCCGTGTCATCC  >  NC_000913/4365312‑4365595

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: