Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,564,635 G→C 6.6% R38P (CGC→CCC)  yffS → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,564,6350GC 6.6% 167.2 / 3.8 61R38P (CGC→CCC) yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base C (2/2);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGT  >  NC_000913/2564489‑2564783
                                                                                                                                                  |                                                                                                                                                    
cccAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGcc                                                                                                                                                    <  1:769646/149‑1 (MQ=255)
   aaTGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCATCAAAGGCGCTGCGATGGAAGGACAATATACGGGCGCTCAACGTGCGGCCAACTGCGAGTAGATCGCAAACGAGTTAACTCGCCggg                                                                                                                                                 >  2:463327/1‑149 (MQ=255)
    aTGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCACACGAGTTAACTCGCCgggg                                                                                                                                                >  1:914022/1‑149 (MQ=255)
    aTGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCgggg                                                                                                                                                >  1:139880/1‑149 (MQ=255)
     tcgATTCAATTAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTTCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGa                                                                                                                                               <  1:648124/147‑1 (MQ=255)
     tGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAAGGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGa                                                                                                                                               <  1:96617/149‑1 (MQ=255)
        aTTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGtt                                                                                                                                            >  1:996661/1‑149 (MQ=255)
         ttCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACCTCCTCAAAGGCGCAGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTa                                                                                                                                           <  1:701437/149‑1 (MQ=255)
           gtATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCCCCGGCGAGTTaac                                                                                                                                         <  2:599243/148‑2 (MQ=255)
             aTGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCCCCGGCGAGTTaactc                                                                                                                                       <  2:254788/149‑4 (MQ=255)
              tGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAaac                                                                                                                                      >  2:56868/1‑149 (MQ=255)
               gagaTGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAaaca                                                                                                                                     >  2:314503/1‑149 (MQ=255)
                    gggAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAg                                                                                                                                <  1:452486/149‑1 (MQ=255)
                      gagaTAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAggg                                                                                                                              >  1:240188/1‑149 (MQ=255)
                        gaTAAAAATCATGAGTTACGAAATCAAAATTAGCGACATCCTCAAAGGCGCTGCGATGGAAGGACATTATAATGGCGCTCAACGTCGGGCCAAATCCGAGGAGATAGCAAACGAGTTAACTCGCCGGGGTCTTAAAAACAACAAGGgag                                                                                                                            >  1:336002/1‑149 (MQ=255)
                         aTAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCACAGGCGCTGCGATGGAAGGACAATATACGGGCGCTCAACTTGGGGGCACATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTAAAAAACAACAAGGaaga                                                                                                                           >  1:561832/1‑149 (MQ=255)
                         aTAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGgaga                                                                                                                           <  1:915560/149‑1 (MQ=255)
                         aTAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGgaga                                                                                                                           <  2:298834/149‑1 (MQ=255)
                         aTAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGgaga                                                                                                                           <  2:551095/149‑1 (MQ=255)
                         aTAAAAATCATGAGTTACGAAATCAAAATTAGCGACATCCTCAAAGGAGCTGCGATGGAAGGACAATTTAAGGGCGCTCAACGTGGGGCCAAATGCGAGGATATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGgaga                                                                                                                           <  1:940205/149‑1 (MQ=255)
                          tAAAAATCATGAGTTACGAAATCAAAATTAGCGACATCCTCACAGTCGCTGCGAGGGAAGGACAATATACGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGagag                                                                                                                          <  1:671002/149‑1 (MQ=255)
                           aaaaaTCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAgg                                                                                                                         >  1:792288/1‑149 (MQ=255)
                                   tGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGTGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACt                                                                                                                 <  2:892375/149‑1 (MQ=255)
                                     cgTTACGAAATCAAAATTTCCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTaa                                                                                                               <  2:451896/148‑1 (MQ=255)
                                            aaaTCAAAATTTGCGCCAACCTCAAATGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGggcggc                                                                                                        <  2:1017507/149‑1 (MQ=255)
                                            aaaTCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAAGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGggcggc                                                                                                        <  2:1018964/149‑1 (MQ=255)
                                                   aaTTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGcc                                                                                                 >  1:325785/1‑149 (MQ=255)
                                                      ttGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCAtt                                                                                              >  1:406759/1‑149 (MQ=255)
                                                      ttGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGGTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCAtt                                                                                              <  2:596806/149‑1 (MQ=255)
                                                       tGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTg                                                                                             <  1:56868/149‑1 (MQ=255)
                                                        gCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTgg                                                                                            <  1:153429/149‑1 (MQ=255)
                                                           aCATCCTCAAAGCCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTg                                                                                         <  2:165060/149‑1 (MQ=255)
                                                            cATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCGCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTgg                                                                                        <  1:209647/149‑1 (MQ=255)
                                                             aTCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGa                                                                                       >  1:340904/1‑149 (MQ=255)
                                                              tCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGTCGTTAGCCATTGGCTGGaa                                                                                      <  2:992488/149‑1 (MQ=255)
                                                                cTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAgg                                                                                    >  1:941231/1‑149 (MQ=255)
                                                                 tCAAAGGCGCTGCGATGGAAGGCCAATATAAGGGCGCTCTACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTAAACTCGCCTGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGt                                                                                   <  2:195185/149‑1 (MQ=255)
                                                                             cGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCa                                                                       >  1:383440/1‑149 (MQ=255)
                                                                                     ggACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGa                                                               <  1:862225/149‑1 (MQ=255)
                                                                                     ggACAAAATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGa                                                               <  2:405940/149‑1 (MQ=255)
                                                                                      gACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGac                                                              >  2:873035/1‑149 (MQ=255)
                                                                                             taAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTgg                                                       <  1:873035/149‑1 (MQ=255)
                                                                                              aaGGGCGCCCAACGTGGGGCCAAATGCGAGGAGCTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGc                                                      <  1:271770/149‑1 (MQ=255)
                                                                                               cgggCGCTCAACGTGGGGCCAAATGCGAGGAGCTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACATGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCGGGAAGGTTGACGGGAGCCAAATTTCGACACTCTGGCt                                                     <  2:410430/148‑1 (MQ=255)
                                                                                                        aaCGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGATGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGc                                            >  1:389059/1‑149 (MQ=255)
                                                                                                        aaCGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGc                                            >  2:635985/1‑149 (MQ=255)
                                                                                                          cGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGa                                          <  2:505397/149‑1 (MQ=255)
                                                                                                           gTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGCGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGAc                                         <  2:732050/149‑1 (MQ=255)
                                                                                                                gCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGtt                                    >  1:881825/1‑149 (MQ=255)
                                                                                                                   aaaTGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGTGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCgg                                 >  1:506746/1‑149 (MQ=255)
                                                                                                                          aggtgATCGCAAACGAGTTAACTGGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTAt                          <  1:755199/149‑1 (MQ=255)
                                                                                                                             agctCGCAAACGAGTTAACTCGACTGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGcc                       <  2:792288/146‑1 (MQ=255)
                                                                                                                                      aCGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCg              >  1:349456/1‑149 (MQ=255)
                                                                                                                                        tttttAACTCCCCGGCGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGAt            >  2:814221/4‑149 (MQ=255)
                                                                                                                                        tttttAACTCCCCGGCGAGTTAAAAACAACAAGGGAGAGGTTATAACTAAGGGCGTCGTTAGCCATTGGCTCGAAGCTAGCCGGGAGCCAAATTTCGACACTCTGCCTGAACTGTGCGACATGTTCGGCGTTTATTCACTCATGCCGaa            >  2:834511/4‑148 (MQ=255)
                                                                                                                                        gAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGAt            >  2:389562/1‑149 (MQ=255)
                                                                                                                                        gAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGAt            <  2:383440/149‑1 (MQ=255)
                                                                                                                                        gAGGTAACTCGCCGGGGAGTTAAAAACTACAAGGGAGAGGGTATTACTCAGGGCGGCGGTAGCCATTGGCTGGAAGGGAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTAGGCCCTGATGCCGAt            <  2:996661/149‑1 (MQ=255)
                                                                                                                                           ttAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACACTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCg         >  2:17937/1‑149 (MQ=255)
                                                                                                                                              aatCGCCGGGGAGTTAAAAACATCAGGGTAGATGTTATTCCTAAGGGCCGAGTTAGCCATTGGCTGGATGGTAGACGTGAGCCAAAATTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTgg      <  1:993385/147‑1 (MQ=255)
                                                                                                                                              aCTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTgg      >  2:825939/1‑149 (MQ=255)
                                                                                                                                                  gCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGt  <  2:941231/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGT  >  NC_000913/2564489‑2564783

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: ‑