Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,041,167 A→C 7.0% intergenic (+49/‑208) yedZ → / → zinT inner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,041,1670AC 7.0% 189.3 / 5.4 57intergenic (+49/‑208)yedZ/zinTinner membrane heme subunit for periplasmic YedYZ reductase/zinc and cadmium binding protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (23/30);  new base C (2/2);  total (25/32)
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.44e-01

GCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATAT  >  NC_000913/2040902‑2041333
                                                                                                                                                                                                                                                                         |                                                                                                                                                                      
gcaaaaaCTGGGCAATCAATTGCAACAGTNNNTTCATTTGTGCTATCTTTTCGCTATTCTGGCCCCAATACATTATCTGTGGGCCGTGAAGATTATTTCCCCGCAGCCTCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTGTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                               <  3:151863/272‑1 (MQ=255)
gCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                               <  2:378698/275‑1 (MQ=255)
gCGAAAACTCGGCAAGCACGGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                               <  2:319407/275‑1 (MQ=255)
       cTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGc                                                                                                                                                          <  1:157490/273‑1 (MQ=255)
                   gggCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCa                                                                                                                                            <  2:75971/274‑1 (MQ=255)
                                           caTCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatc                                                                                                                    <  2:137355/274‑1 (MQ=255)
                                              cTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTTCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGTTACCATCAGGTTATGTTTCGTTTGTCATCAAc                                                                                                                 >  4:242027/1‑275 (MQ=255)
                                                                                    cGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGa                                                                                                                                                                        <  1:145613/182‑1 (MQ=255)
                                                                                    cGTGAAGATTATTTCACCGCAGCCGCTCATCCACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCANATTTCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACCCAGATCGCATCAGGTTATNGCTCGGTTGTTATCAACCCCAAATANCCAAAAAAATTTTTTTTTTTTtaataaat                                                                           >  3:87698/1‑273 (MQ=255)
                                                                                                 tCACCGCAGCCGCTCATCCACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGa                                                               <  1:87698/274‑1 (MQ=255)
                                                                                                                                             aGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGcc                                                                                                                                                                   >  4:158045/1‑130 (MQ=255)
                                                                                                                                                 ttACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                               <  2:158046/130‑1 (MQ=255)
                                                                                                                                                     ggTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTTTCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGt          <  2:242028/275‑1 (MQ=255)
                                                                                                                                                           agaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAgcgc                                                                                                                                                                    >  8:391289/1‑115 (MQ=255)
                                                                                                                                                           agaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAgcgc                                                                                                                                                                    <  5:742971/115‑1 (MQ=255)
                                                                                                                                                           agaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAgcgc                                                                                                                                                                    >  6:742971/1‑115 (MQ=255)
                                                                                                                                                           agaagTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAgcgc                                                                                                                                                                    <  7:391289/115‑1 (MQ=255)
                                                                                                                                                             aagTTTCGTTCCCTTTTTAACCGGTTACGCAAGAAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCAGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAatat  <  1:1480/275‑1 (MQ=255)
                                                                                                                                                                 tGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATac                                                                                                                                                      >  4:545491/1‑123 (MQ=255)
                                                                                                                                                                     ttCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAg                                                                                                                                                  <  2:545492/123‑1 (MQ=255)
                                                                                                                                                                              tAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                               >  6:683984/1‑101 (MQ=255)
                                                                                                                                                                              tAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTAt                                                                                                                                                               <  5:683984/101‑1 (MQ=255)
                                                                                                                                                                                   ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGa                                                                                                                                                                        >  7:1080136/1‑87 (MQ=255)
                                                                                                                                                                                   ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGa                                                                                                                                                                        <  6:240872/87‑1 (MQ=255)
                                                                                                                                                                                   ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGa                                                                                                                                                                        >  5:240872/1‑87 (MQ=255)
                                                                                                                                                                                   ggTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGa                                                                                                                                                                        <  8:1080136/87‑1 (MQ=255)
                                                                                                                                                                                     ttACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCg                                                                                                                             >  5:994508/1‑128 (MQ=255)
                                                                                                                                                                                     ttACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCg                                                                                                                             <  6:994508/128‑1 (MQ=255)
                                                                                                                                                                                     ttACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTg                                                                                                                         >  8:1059863/1‑132 (MQ=255)
                                                                                                                                                                                     ttACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTg                                                                                                                         <  7:1059863/132‑1 (MQ=255)
                                                                                                                                                                                      tACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCa                                                                                                        <  7:694754/148‑1 (MQ=255)
                                                                                                                                                                                      tACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCa                                                                                                        >  8:694754/1‑148 (MQ=255)
                                                                                                                                                                                            aGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAg                                                                                                                                        >  8:146192/1‑110 (MQ=255)
                                                                                                                                                                                            aGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAg                                                                                                                                        <  7:146192/110‑1 (MQ=255)
                                                                                                                                                                                            aGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGccacat                                                                                                                                                               >  6:584881/1‑83 (MQ=37)
                                                                                                                                                                                            aGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGccacat                                                                                                                                                               <  5:584881/87‑5 (MQ=37)
                                                                                                                                                                                                 gTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGt                                                                                                                            >  7:940499/1‑117 (MQ=255)
                                                                                                                                                                                                 gTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGt                                                                                                                            <  8:940499/117‑1 (MQ=255)
                                                                                                                                                                                                        acaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGc                                                                                                                                                          <  7:1034716/80‑1 (MQ=25)
                                                                                                                                                                                                        acaaATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATAGGGCGCTCGCGCCACATCCGGc                                                                                                                                                          >  8:1034716/1‑80 (MQ=25)
                                                                                                                                                                                                           aaTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAgcgc                                                                                                                                                                    >  7:935932/1‑67 (MQ=255)
                                                                                                                                                                                                           aaTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAgcgc                                                                                                                                                                    <  8:935932/67‑1 (MQ=255)
                                                                                                                                                                                                           aaTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAAt                                                                                                          >  8:869515/1‑125 (MQ=255)
                                                                                                                                                                                                           aaTTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAAt                                                                                                          <  7:869515/125‑1 (MQ=255)
                                                                                                                                                                                                              tATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTa                                                                                                                                    <  6:104967/96‑1 (MQ=255)
                                                                                                                                                                                                              tATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTa                                                                                                                                    >  5:104967/1‑96 (MQ=255)
                                                                                                                                                                                                                tCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTAt                                                                                                                                   >  6:934228/1‑95 (MQ=255)
                                                                                                                                                                                                                tCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTAt                                                                                                                                   <  5:934228/95‑1 (MQ=255)
                                                                                                                                                                                                                tCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtca                                                                                                                      <  7:501401/108‑1 (MQ=255)
                                                                                                                                                                                                                tCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtca                                                                                                                      >  8:501401/1‑108 (MQ=255)
                                                                                                                                                                                                                          ttCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatca                                                                                                                   >  6:672825/1‑101 (MQ=255)
                                                                                                                                                                                                                          ttCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGtcatca                                                                                                                   <  5:672825/101‑1 (MQ=255)
                                                                                                                                                                                                                           tCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAg                                                                   >  5:256852/1‑148 (MQ=255)
                                                                                                                                                                                                                              tatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGt                                                                                                                                      >  7:337859/1‑78 (MQ=255)
                                                                                                                                                                                                                              tatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGt                                                                                                                                      <  8:337859/78‑1 (MQ=255)
                                                                                                                                                                                                                              tatTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAAcc                                                                                                                >  2:555056/1‑100 (MQ=255)
                                                                                                                                                                                                                                  gCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCaa                                                                                                            <  4:555055/100‑1 (MQ=255)
                                                                                                                                                                                                                                    cGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGccc                                                                                                                                                                  <  8:589956/44‑1 (MQ=255)
                                                                                                                                                                                                                                    cGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGccc                                                                                                                                                                  >  7:589956/1‑44 (MQ=255)
                                                                                                                                                                                                                                       cAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGaa                                                       <  6:256852/148‑1 (MQ=255)
                                                                                                                                                                                                                                                     tCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGttt                                                                                                                          <  8:1007915/67‑1 (MQ=255)
                                                                                                                                                                                                                                                     tCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGttt                                                                                                                          >  7:1007915/1‑67 (MQ=255)
                                                                                                                                                                                                                                                                         |                                                                                                                                                                      
GCGAAAACTCGGCAAGCACTGGCAACAGTTGCATAATTTCGTCTATCTTGTCGCGATTCTGGCCCCAATACATTATCTGTGGTCCGTGAAGATTATTTCACCGCAGCCGCTCATCTACGCTGGGCTGGCTGTACTGCTTTTAGCCTTACGGTATAAGAAGTTGCGTTCCCTTTTTAACCGGTTACGCAAGCAGGTTCACAACAAATTATCTGTATAACTTCATATTGCCGACAAAGTACGCGTTATCCCTGCCGGATGTGGCGCGAGCGCCCTATCCGGCATACACAGATCGCATCAGGTTATGTTTCGTTTGTCATCAACCTCAAATCACCAAAATAATTTATTTTATTTCAATAAGATAGATAAGTAGAACTGAGAAAGCCATGCTCTCGTTTCCTAAGAGTTGTTGCATTTTGCTATATGTTACAATAT  >  NC_000913/2040902‑2041333

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