Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,564,635 G→C 5.0% 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 5.0% 350.6 / 8.6 99R38P (CGC→CCC) yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (46/48);  new base C (2/3);  total (48/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 = 5.11e-01

TCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCG  >  NC_000913/2564488‑2564781
                                                                                                                                                   |                                                                                                                                                  
tCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGc                                                                                                                                                   >  2:66637/1‑149 (MQ=255)
tCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGc                                                                                                                                                   >  4:164503/1‑149 (MQ=255)
tCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGc                                                                                                                                                   >  2:69977/1‑149 (MQ=255)
   cAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCgg                                                                                                                                                >  3:475036/1‑149 (MQ=255)
      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                                                                                                                                     <  3:229130/149‑1 (MQ=255)
              aTGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTaaaaa                                                                                                                                     <  1:475883/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                                                                                                                             <  2:290194/149‑1 (MQ=255)
                      ggAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAgg                                                                                                                             >  3:381950/1‑149 (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                                                                                        >  4:607156/1‑145 (MQ=255)
                                                               tCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCt                                                                                        <  3:607156/145‑1 (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                                                                 >  4:653510/1‑149 (MQ=255)
                                                                                  ggAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAAtt                                                                 >  1:622543/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                          <  3:164503/149‑1 (MQ=255)
                                                                                                                         cGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGttt                          <  1:69977/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                >  4:543336/1‑149 (MQ=255)
                                                                                                                                   gCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGAt                <  4:592541/149‑1 (MQ=255)
                                                                                                                                   gCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGAt                <  3:683693/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)
                                                                                                                                                   |                                                                                                                                                  
TCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGCAAACGAGTTAACTCGCCGGGGAGTTAAAAACAACAAGGGAGAGGTTATTACTAAGGGCGGCGTTAGCCATTGGCTGGAAGGTAGACGGGAGCCAAATTTCGACACTCTGGCTGAACTGTGCGACATGTTCGGCGTTTATGCCCTGATGCCGATGCGTGGCG  >  NC_000913/2564488‑2564781

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: