Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,564,640 G→C 6.9% G40R (GGA→CGA)  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,6400GC 6.9% 155.3 / 4.0 58G40R (GGA→CGA) yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (21/33);  new base C (2/2);  total (23/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.94e-01

AATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAATG  >  NC_000913/2564492‑2564788
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                                <  2:551095/149‑1 (MQ=255)
                      aTAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGgaga                                                                                                                                <  1:915560/149‑1 (MQ=255)
                      aTAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGgaga                                                                                                                                <  2:298834/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:383440/149‑1 (MQ=255)
                                                                                                                                     gAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGAt                 >  2:389562/1‑149 (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)
                                                                                                                                                   gggACTTAACAACAACACGTGAGCGGTTATTACAAAGGGCGTCCTTAGCCATTGGCGGGAAGGTAGACGGCAGCCAAATTTCGACACTCTGGCTGACCTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAAt   <  2:1027543/149‑1 (MQ=255)
                                                                                                                                                    ggAGTTAAACACAACAAGGGAGAGGTAATTACTACGGGCCGCGTTAGCCATTGGCTGGAAGGTAGACGGGTGCCAAATTTCGACTCTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGAAAATg  <  2:880784/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAATG  >  NC_000913/2564492‑2564788

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