Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,992,827 T→C 5.9% intergenic (‑10/+47) pgsA ← / ← uvrC phosphatidylglycerophosphate synthetase/excinuclease UvrABC, endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,992,8270TC 5.9% 207.4 / 3.0 68intergenic (‑10/+47)pgsA/uvrCphosphatidylglycerophosphate synthetase/excinuclease UvrABC, endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base C (2/2);  total (35/33)
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.32e-01

TCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTTACGTAAACCTTGCA  >  NC_000913/1992685‑1992971
                                                                                                                                              |                                                                                                                                                
tcaCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTAt                                                                                                                                            >  3:186135/1‑149 (MQ=255)
 caCCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATc                                                                                                                                           >  7:250411/1‑149 (MQ=255)
   ccgccgCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCtg                                                                                                                                         <  7:158397/149‑1 (MQ=255)
       cgcgACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCtgttgt                                                                                                                                     <  7:237206/149‑1 (MQ=255)
           aCGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGt                                                                                                                                 <  4:186135/149‑1 (MQ=255)
             gctgAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTaa                                                                                                                               <  8:388925/146‑1 (MQ=255)
             gCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTaa                                                                                                                               <  1:160198/149‑1 (MQ=255)
               aGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGt                                                                                                                                 >  6:28499/1‑145 (MQ=255)
                  aaaTGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTa                                                                                                                          >  8:251604/1‑149 (MQ=255)
                    aTGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTAcc                                                                                                                        <  6:165983/149‑1 (MQ=255)
                    aTGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTAcc                                                                                                                        <  8:250411/149‑1 (MQ=255)
                         cgaggcggcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAATCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTAt                                                                                                                   <  6:169756/146‑1 (MQ=255)
                                gcAAACGGCGACCAGGTGACAGGCAGATAAAAGCCCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTa                                                                                                            >  8:87056/1‑149 (MQ=255)
                                gcAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTa                                                                                                            <  5:28499/149‑1 (MQ=255)
                                 cAAACGGCGACCATGTGACAGGCAGATAAAAGACCAATACAAAGATTGGGATATGGATGACACGGCACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATGTGTTGTCAGTAAGATTACCCCTTTGTTGCTAc                                                                                                           <  5:79303/149‑1 (MQ=255)
                                  aaaCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACTTAGGGATATTAAATTGCATAATGACGGGTAACTATCtgt                                                                                                                                        <  1:289799/119‑1 (MQ=255)
                                  aaaCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTa                                                                                                            >  8:215106/1‑147 (MQ=255)
                                  aaaCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACATAGGGATATTAAATTGCATAATGACGGGTAACTATCtgt                                                                                                                                        >  2:289799/1‑119 (MQ=255)
                                    aCGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACaga                                                                                                        <  1:270296/149‑1 (MQ=255)
                                           ccAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCa                                                                                                 <  6:350330/149‑1 (MQ=255)
                                               gTGACAGGCAGATAAAAGCCCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGt                                                                                             >  4:297437/1‑149 (MQ=255)
                                                    aGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCaa                                                                                        >  7:375259/1‑149 (MQ=255)
                                                    aGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCaa                                                                                        >  8:243431/1‑149 (MQ=255)
                                                        aGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATATGTTGTCATTAAGATTACCCCTATGTTGCTACAGATACATCTATGTTTCAACGAc                                                                                    >  3:270493/1‑149 (MQ=255)
                                                          aTAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGAc                                                                                    >  1:134873/1‑147 (MQ=255)
                                                           tAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGTATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGAc                                                                                    >  6:445/1‑146 (MQ=255)
                                                            aaaaGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGATGTCAGTAAGATTACCCCTATGTTGCTACAgaga                                                                                                      >  3:212039/1‑127 (MQ=255)
                                                            aaaaGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGATGTCAGTAAGATTACCCCTATGTTGCTACAgaga                                                                                                      <  4:212039/127‑1 (MQ=255)
                                                                gACCAATACAAAGAATGGGATAAGGATGACACTGAACAGTTTAAGCAACGTAGCGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGAt                                                                            <  5:19689/149‑1 (MQ=255)
                                                                   cAATACAAAGAATGGGATTAGGATGACAAGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAAGTATCTTTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACTACCAGAAGATGtt                                                                         >  2:360314/1‑149 (MQ=255)
                                                                    aaTACAAAGAATGGGATAATGATGACACGGAACTGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGTGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCttt                                                                        <  1:331332/149‑1 (MQ=255)
                                                                    aaTACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCttt                                                                        <  8:375259/149‑1 (MQ=255)
                                                                     aTACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGTATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCtttt                                                                       <  5:445/149‑1 (MQ=255)
                                                                     aTACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCtttt                                                                       <  7:251604/149‑1 (MQ=255)
                                                                     aTACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTAACCCTAAGTTGCTACAGAGACATCAATGTTTCAACGACAAGAAGATCttt                                                                        >  3:312227/1‑148 (MQ=255)
                                                                         aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCCACGTAGGGATATTAAATTGCATAATGACGGGCAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGc                                                                   >  4:26392/1‑149 (MQ=255)
                                                                         aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTa                                                                                                            <  7:345054/108‑1 (MQ=255)
                                                                         aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTa                                                                                                            >  8:345054/1‑108 (MQ=255)
                                                                         aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCtt                                                                         >  7:271178/1‑143 (MQ=255)
                                                                         aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCtt                                                                         <  8:271178/143‑1 (MQ=255)
                                                                         aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGc                                                                   >  2:398148/1‑149 (MQ=255)
                                                                         aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGc                                                                   >  2:323081/1‑149 (MQ=255)
                                                                         aaaGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGCAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGc                                                                   >  4:21369/1‑149 (MQ=255)
                                                                               tGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACaga                                                                                                        <  5:203214/106‑1 (MQ=255)
                                                                               tGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACaga                                                                                                        >  6:203214/1‑106 (MQ=255)
                                                                                        gATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAAt                                                    <  2:134873/149‑1 (MQ=255)
                                                                                            acacGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATAcc                                                 <  8:85508/148‑1 (MQ=255)
                                                                                             cacGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGATGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGCACAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCg                                               >  4:91382/1‑149 (MQ=255)
                                                                                                      gtgtAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTg                                      >  5:82797/1‑149 (MQ=255)
                                                                                                           aGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAAtt                                 >  1:56442/1‑149 (MQ=255)
                                                                                                           aGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAAt                                  <  8:132808/148‑1 (MQ=255)
                                                                                                                  tAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGa                          <  7:243431/149‑1 (MQ=255)
                                                                                                                   aGGGATATTAAATTGCATAATGCCGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGCCCAGAAGATCTTTTCTGCCAGACCTTGCGAACTACCCGGCACTTATGCAAATTCCTCGAc                         >  8:125803/1‑149 (MQ=255)
                                                                                                                    gggATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGAc                         <  1:398148/148‑1 (MQ=255)
                                                                                                                      gATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCt                      >  6:76460/1‑149 (MQ=255)
                                                                                                                          ttAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCg                  >  4:328999/1‑149 (MQ=255)
                                                                                                                            aaaTTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGt                 >  8:129855/1‑148 (MQ=255)
                                                                                                                             aaTTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTTa               >  1:152391/1‑149 (MQ=255)
                                                                                                                              aTTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTTAc              <  6:82797/149‑1 (MQ=255)
                                                                                                                              aTTGCATAATGACGGGCAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTTAc              <  3:21369/149‑1 (MQ=255)
                                                                                                                              aTTGCATAATGACGGGCAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTGCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTTAc              <  3:26392/149‑1 (MQ=255)
                                                                                                                                 gCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATAACCGGCACTTTTGCAATTTCCTCGACGCTGGCGTTACGt            >  1:348669/1‑148 (MQ=255)
                                                                                                                                       tGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGCACTTGCGAAATACCCGGCCCTTTTGCAATTTCCTCGACGCTGGCGTTACGTAAACCtt     >  6:245745/1‑149 (MQ=255)
                                                                                                                                         aCGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTTACGTAAACCTTg    <  7:129855/148‑1 (MQ=255)
                                                                                                                                          cGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTTACGTAAACCTTGc   >  5:233379/1‑148 (MQ=255)
                                                                                                                                          cGGGTAACTAGATGTTGTAAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCTGCACTTTTGCAATTTCCTCGACGCTGGCGTTACGTAAACCTTGCa  <  4:270493/149‑1 (MQ=255)
                                                                                                                                             gTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGc                   <  3:290058/129‑1 (MQ=255)
                                                                                                                                             gTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGc                   >  4:290058/1‑129 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TCACCGCCGCGACGCAGAAAATGAGCGCGGCGGCAAACGGCGACCAGGTGACAGGCAGATAAAAGACCAATACAAAGAATGGGATAAGGATGACACGGAACAGTGTAAGCAACGTAGGGATATTAAATTGCATAATGACGGGTAACTATCTGTTGTCAGTAAGATTACCCCTATGTTGCTACAGAGACATCAATGTTTCAACGACCAGAAGATCTTTTCTGCCAGACCTTGCGAAATACCCGGCACTTTTGCAATTTCCTCGACGCTGGCGTTACGTAAACCTTGCA  >  NC_000913/1992685‑1992971

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