Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,564,640 G→C 5.3% 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 5.3% 339.7 / 8.7 95G40R (GGA→CGA) yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (42/48);  new base C (2/3);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

TGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAA  >  NC_000913/2564494‑2564786
                                                                                                                                                  |                                                                                                                                                  
tGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGa                                                                                                                                                  <  1:252375/149‑1 (MQ=255)
tGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGa                                                                                                                                                  <  3:118204/149‑1 (MQ=255)
  gATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGt                                                                                                                                                >  3:255591/1‑149 (MQ=255)
   aTTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGtt                                                                                                                                               <  2:555930/149‑1 (MQ=255)
       tATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTaaaa                                                                                                                                           >  2:307647/1‑149 (MQ=255)
       tATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGACCGCAAACGAGTTAACTCCCCGGCGAGTTaact                                                                                                                                           <  1:334029/149‑3 (MQ=255)
        aTGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTaaaaa                                                                                                                                          <  1:475883/149‑1 (MQ=255)
        aTGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTaaaaa                                                                                                                                          <  3:229130/149‑1 (MQ=255)
        aTGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCCCCGGCGAGTTaactc                                                                                                                                          <  3:733715/149‑4 (MQ=255)
         tGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAaac                                                                                                                                         >  1:64138/1‑149 (MQ=255)
          gagaTGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAaaca                                                                                                                                        >  3:271980/1‑149 (MQ=255)
                ggAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAgg                                                                                                                                  >  3:381950/1‑149 (MQ=255)
                ggAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAgg                                                                                                                                  <  2:290194/149‑1 (MQ=255)
                  agaTAAAAATCATGAGTTACGAAATCAAAACTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGga                                                                                                                                <  3:237264/149‑1 (MQ=255)
                   gacaaaagTCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGCCAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGgag                                                                                                                               >  4:487589/4‑149 (MQ=255)
                     tAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGAGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGagag                                                                                                                             <  2:130934/149‑1 (MQ=255)
                      aaaaaTCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAgg                                                                                                                            >  2:537437/1‑149 (MQ=255)
                      aaaaaTCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACAAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAgg                                                                                                                            >  4:151229/1‑149 (MQ=255)
                       aaaaTCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGGTGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGt                                                                                                                           >  2:661465/1‑149 (MQ=255)
                       aaaaTCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGt                                                                                                                           <  1:127103/149‑1 (MQ=255)
                          aTCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGttat                                                                                                                        <  1:412721/149‑1 (MQ=255)
                           tCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGttatt                                                                                                                       <  4:271980/149‑1 (MQ=255)
                              tGAGTTACGAAAACAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACt                                                                                                                    <  4:160156/149‑1 (MQ=255)
                                aGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTaa                                                                                                                  <  3:151229/149‑1 (MQ=255)
                                    aCGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGggc                                                                                                              <  3:621368/149‑1 (MQ=255)
                                       aaaTCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCCGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGggcggc                                                                                                           <  4:607076/149‑1 (MQ=255)
                                         aTCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGt                                                                                                         >  4:231561/1‑149 (MQ=255)
                                          tCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGtt                                                                                                        >  4:489556/1‑149 (MQ=255)
                                            aaaaTTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAg                                                                                                      >  2:376395/1‑149 (MQ=255)
                                             aaaTTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTTCTAAGGGCGGCGTTAGc                                                                                                     <  4:515020/149‑1 (MQ=255)
                                                tttGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCAt                                                                                                  <  4:272887/149‑1 (MQ=255)
                                                  tGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTg                                                                                                <  2:177191/149‑1 (MQ=255)
                                                    cGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGc                                                                                              <  2:66467/149‑1 (MQ=255)
                                                     gACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCt                                                                                             >  2:671358/1‑149 (MQ=255)
                                                         tCCTCAAAGGCGCTGCGTTGGAAGGACAATATAAGGGCGCTTAACGTGGTGCCAAATGCGAGGAGATCGAAAACGAGTTAACTCGCCGGGGACTTAAAAACAACAAGGGAGAGGTTATCACTAAGGGCGGCGTTAGCTATTGGCTGGaa                                                                                         >  4:468586/1‑149 (MQ=255)
                                                         tCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCt                                                                                             <  3:607156/145‑1 (MQ=255)
                                                         tCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCt                                                                                             >  4:607156/1‑145 (MQ=255)
                                                         tCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGaa                                                                                         <  2:494232/149‑1 (MQ=255)
                                                           cTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGt                                                                                                                                                >  4:634346/1‑92 (MQ=255)
                                                           cTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGt                                                                                                                                                <  3:634346/92‑1 (MQ=255)
                                                            tCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGt                                                                                      >  3:142492/1‑149 (MQ=255)
                                                                  gcgcTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACgg                                                                                <  2:545439/149‑1 (MQ=255)
                                                                    gcTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGa                                                                              <  4:276877/149‑1 (MQ=255)
                                                                    gcTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGa                                                                              <  4:475036/149‑1 (MQ=255)
                                                                           tGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGATTTAAAAACAACAAgg                                                                                                                                  <  4:634740/90‑1 (MQ=255)
                                                                           tGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGACTTAAAAACAACAAgg                                                                                                                                  >  3:634740/1‑90 (MQ=255)
                                                                            ggAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAAtt                                                                      >  1:622543/1‑149 (MQ=255)
                                                                            ggAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAAtt                                                                      >  4:653510/1‑149 (MQ=255)
                                                                             gAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAAttt                                                                     >  4:114842/1‑149 (MQ=255)
                                                                               aGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCg                                                                   >  2:594465/1‑149 (MQ=255)
                                                                                ggACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCGTTGGCTGGAAGGTAGACGGGAGCCAAATTTCGa                                                                  <  2:101563/149‑1 (MQ=255)
                                                                                ggACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGa                                                                  <  4:708987/149‑1 (MQ=255)
                                                                                 gACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGac                                                                 >  4:710790/1‑149 (MQ=255)
                                                                                 gACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGac                                                                 >  3:758306/1‑149 (MQ=255)
                                                                                       ataAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTg                                                           <  4:574971/149‑1 (MQ=255)
                                                                                              cgcTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAAc                                                    <  3:13158/149‑1 (MQ=255)
                                                                                               gcTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACt                                                   >  3:529140/1‑149 (MQ=255)
                                                                                                   aaCGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGc                                               <  1:211952/149‑1 (MQ=255)
                                                                                                         gggCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATg                                         >  4:681242/1‑149 (MQ=255)
                                                                                                          ggCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGt                                        >  1:639247/1‑149 (MQ=255)
                                                                                                             cAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCg                                     >  4:683693/1‑149 (MQ=255)
                                                                                                                aTGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTAATAAGGGCGGCGTTAGCCATGGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCg                                  >  3:466618/1‑149 (MQ=255)
                                                                                                                 tGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGt                                 <  4:758306/149‑1 (MQ=255)
                                                                                                                  gCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGtt                                <  4:658080/149‑1 (MQ=255)
                                                                                                                   cGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGttt                               <  1:69977/149‑1 (MQ=255)
                                                                                                                   cGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGttt                               <  3:164503/149‑1 (MQ=255)
                                                                                                                     aggagATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTAt                             <  1:273379/149‑1 (MQ=255)
                                                                                                                         gaTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGccc                         <  4:282548/149‑1 (MQ=255)
                                                                                                                          aTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAgg                                                                                                                                  >  3:377847/1‑43 (MQ=255)
                                                                                                                          aTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAgg                                                                                                                                  <  4:377847/43‑1 (MQ=255)
                                                                                                                          aTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTa                                                                                     <  1:246200/88‑1 (MQ=255)
                                                                                                                          aTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTa                                                                                     >  2:246200/1‑88 (MQ=255)
                                                                                                                          aTCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCt                        >  1:531604/1‑149 (MQ=255)
                                                                                                                           tCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACtgt                                                 >  1:27354/1‑123 (MQ=255)
                                                                                                                           tCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACtgt                                                 <  2:27354/123‑1 (MQ=255)
                                                                                                                           tCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTg                       >  1:169065/1‑149 (MQ=255)
                                                                                                                           tCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTg                       >  1:294309/1‑149 (MQ=255)
                                                                                                                           tCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTg                       >  4:485031/1‑149 (MQ=255)
                                                                                                                             gCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGAt                     <  3:683693/149‑1 (MQ=255)
                                                                                                                             gCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGAt                     >  4:543336/1‑149 (MQ=255)
                                                                                                                             gCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGAt                     <  4:592541/149‑1 (MQ=255)
                                                                                                                               aaaCGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGc                   >  3:50716/1‑149 (MQ=255)
                                                                                                                                  cGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGa                <  1:376395/149‑1 (MQ=255)
                                                                                                                                   tttttAACTCCCCGGCGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGAt               >  2:442340/4‑149 (MQ=255)
                                                                                                                                   tttttAACTCCCCGGCGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGAt               <  1:442340/146‑1 (MQ=255)
                                                                                                                                   gAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGAt               <  2:169065/149‑1 (MQ=255)
                                                                                                                                    ttttAACTCCCCGGCGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATg              >  2:766932/3‑149 (MQ=255)
                                                                                                                                     gTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGc             <  3:653510/149‑1 (MQ=255)
                                                                                                                                      ttAANTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCg            >  3:697029/1‑149 (MQ=255)
                                                                                                                                       tAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGt           <  4:255591/149‑1 (MQ=255)
                                                                                                                                        aaCTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTg          >  2:135486/1‑149 (MQ=255)
                                                                                                                                         aCTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTgg         >  3:394775/1‑149 (MQ=255)
                                                                                                                                         aCTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATACCCTGATGCCGATGCGTgg         >  2:401237/1‑149 (MQ=255)
                                                                                                                                          cTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGc        <  1:537437/149‑1 (MQ=255)
                                                                                                                                           tCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCg       <  4:110947/149‑1 (MQ=255)
                                                                                                                                                ggggAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTaaa  <  3:489556/149‑1 (MQ=255)
                                                                                                                                                ggggAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTaaa  <  1:66637/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCGGTAAA  >  NC_000913/2564494‑2564786

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