Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,002,528 A→T 6.2% I105F (ATT→TTT)  PP_4409 → site‑specific recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,002,5280AT 6.2% 249.2 / 5.3 80I105F (ATT→TTT) PP_4409site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base T (2/3);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

ACATAATGCCAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCC  >  NC_002947/5002386‑5002673
                                                                                                                                               |                                                                                                                                                  
aCATAATGCCAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGc                                                                                                                                               >  3:108023/1‑148 (MQ=255)
  aTAATGCCAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGaa                                                                                                                                            <  1:16265/149‑1 (MQ=255)
    aaTGCCAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGc                                                                                                                                          <  8:142875/149‑1 (MQ=255)
      tGCCAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCATTTCTACAGTTGTTCGCAAAATTTGCGAAGCt                                                                                                                                         >  3:63397/1‑148 (MQ=255)
       gCCAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAATTTTGCGaacaa                                                                                                                                         >  5:329566/1‑144 (MQ=255)
       gCCAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAATTTTGCGaacaa                                                                                                                                         <  6:329566/147‑4 (MQ=255)
        ccAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAATTTTGCGaacaactg                                                                                                                                      <  8:161805/149‑7 (MQ=255)
        ccAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTaaa                                                                                                                                      <  6:135602/149‑1 (MQ=255)
           aCTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCACGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCg                                                                                                                                   >  6:10170/1‑149 (MQ=255)
             ttCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTc                                                                                                                                 >  1:137264/1‑149 (MQ=255)
                    aCTTCTTCGAAGAAGGTATCAATCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATg                                                                                                                          <  8:212929/149‑1 (MQ=255)
                     cttcttCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATgg                                                                                                                         >  1:78924/1‑149 (MQ=255)
                      ttcttcGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGt                                                                                                                        <  2:333849/149‑1 (MQ=255)
                        cttcGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTcc                                                                                                                                >  7:14477/1‑139 (MQ=255)
                        cttcGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTcc                                                                                                                                <  8:14477/139‑1 (MQ=255)
                         ttcGAAGAAGGTCTCATGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAg                                                                                                                     >  3:158224/1‑149 (MQ=255)
                         ttcGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAg                                                                                                                     >  7:117491/1‑149 (MQ=255)
                            gaagaaGGTCTCAAGCATATTGAGCCGAGAAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAACCAc                                                                                                                  >  8:26139/1‑149 (MQ=255)
                              agaagGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACgg                                                                                                                >  8:139519/1‑149 (MQ=255)
                              agaagGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACgg                                                                                                                >  7:177484/1‑149 (MQ=255)
                                     ctcAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCAc                                                                                                                  >  5:373549/1‑140 (MQ=255)
                                     ctcAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCAc                                                                                                                  <  6:373549/140‑1 (MQ=255)
                                        aaGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATg                                                                                                                          <  3:298370/129‑1 (MQ=255)
                                             tatTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGc                                                                                                                                               <  5:380215/103‑1 (MQ=255)
                                             tatTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGc                                                                                                                                               >  6:380215/1‑103 (MQ=255)
                                             tatTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAat                                                                                                 <  4:108023/149‑1 (MQ=255)
                                               ttGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGAAAGCACGGCATCGTCAAAAAAatat                                                                                               >  8:228055/1‑149 (MQ=255)
                                                  aGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCa                                                                                             <  8:117491/148‑1 (MQ=255)
                                                  aGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCa                                                                                             <  4:304077/148‑1 (MQ=255)
                                                      gAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCCTCGTCAAAAAAATATCACGCt                                                                                         >  6:197204/1‑148 (MQ=255)
                                                        gCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCaa                                                                                      <  6:29832/149‑1 (MQ=255)
                                                             tCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTc                                                                                  >  1:308277/1‑147 (MQ=255)
                                                             tCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTc                                                                                  >  4:278311/1‑147 (MQ=255)
                                                             tCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCa                                                                                <  7:26139/149‑1 (MQ=255)
                                                                       acacCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCt                                                                       >  5:46709/1‑148 (MQ=255)
                                                                         acCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTg                                                                      >  6:344832/1‑147 (MQ=255)
                                                                           ccTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAg                                                                               >  2:79594/1‑136 (MQ=255)
                                                                           ccTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAg                                                                               <  1:79594/136‑1 (MQ=255)
                                                                           ccTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAaga                                                                  <  7:139519/149‑1 (MQ=255)
                                                                             tGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGc                                                                >  6:317517/1‑149 (MQ=255)
                                                                                aTCGAGAGTATCACTCATACTTGACAGTTGGCAGCACGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCt                                                              <  5:10170/148‑1 (MQ=255)
                                                                                   gagaGTATCACTCATACTTGACAGTTGGCAGCGACGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTgg                                                           >  1:43591/1‑149 (MQ=255)
                                                                                     gagTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTaa                                                        <  2:78924/149‑1 (MQ=255)
                                                                                            aCTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTc                                                 >  2:100700/1‑149 (MQ=255)
                                                                                             cTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTcc                                                <  5:344832/149‑1 (MQ=255)
                                                                                             cTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTc                                                 >  4:284700/1‑148 (MQ=255)
                                                                                              tCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCGCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTc                                                 >  5:375693/1‑147 (MQ=255)
                                                                                              tCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCCAGCCTTTCCTGAAGAGCCTTGGTaa                                                        <  1:192199/141‑1 (MQ=255)
                                                                                              tCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTaa                                                        >  2:192199/1‑140 (MQ=255)
                                                                                              tCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTc                                                 >  4:163584/1‑147 (MQ=255)
                                                                                              tCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCAAGCATTTCCTGAAGAGCCTTGGTAATTACTTcc                                                >  1:156675/1‑148 (MQ=255)
                                                                                                  aCTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGaaa                                           >  5:243644/1‑149 (MQ=255)
                                                                                                   cTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAAt                                          <  8:177484/149‑1 (MQ=255)
                                                                                                    ttGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAAcc                                                                                                                                    >  5:371492/1‑59 (MQ=255)
                                                                                                    ttGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAAcc                                                                                                                                    <  6:371492/59‑1 (MQ=255)
                                                                                                    ttGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATa                                         >  2:293011/1‑149 (MQ=255)
                                                                                                    ttGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATa                                         >  6:342781/1‑149 (MQ=255)
                                                                                                    ttGACAGTTGGCAGCGAGAGCAATTCAAGAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCACCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGAACAAGTCCAGCCTTTCCTGAAGAGCCTTGGTATTTACTTCCAGAAAt                                          >  1:209563/1‑148 (MQ=255)
                                                                                                      gACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATAc                                        >  2:285487/1‑148 (MQ=255)
                                                                                                       aCAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACa                                       >  3:184812/1‑148 (MQ=255)
                                                                                                           ttGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGac                                  >  2:293488/1‑149 (MQ=255)
                                                                                                            tGGGAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGac                                  <  5:280043/148‑1 (MQ=255)
                                                                                                               cAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGactact                              <  4:184812/149‑1 (MQ=255)
                                                                                                                 gcgcGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTACCTGAAGAGACTTGGTAATTACTTCCAGAAATACAAAGactact                              >  5:190997/1‑147 (MQ=255)
                                                                                                                         cAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACg                     >  5:84943/1‑149 (MQ=255)
                                                                                                                          aaTTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATTTCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACgg                    >  5:185440/1‑149 (MQ=255)
                                                                                                                               aaCAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTg                <  5:342781/148‑1 (MQ=255)
                                                                                                                                    tttttCTCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCAGGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAAtt          <  2:156675/146‑1 (MQ=255)
                                                                                                                                    ttgttCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAAtt          <  2:43591/149‑1 (MQ=255)
                                                                                                                                    ttgttCGCAAAATTTGCGAAGCTAAACCGTCCCCAATCGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAAtt          <  4:63397/149‑1 (MQ=255)
                                                                                                                                     tgttCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAAtt          <  4:158224/148‑1 (MQ=255)
                                                                                                                                     tgttCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAAtt          <  1:293011/148‑1 (MQ=255)
                                                                                                                                     tgttCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAAtt          <  6:243644/148‑1 (MQ=255)
                                                                                                                                      cttCGCAAATTTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTg         >  4:244701/2‑148 (MQ=255)
                                                                                                                                          gCAAATTTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGaaaa     <  3:244701/148‑1 (MQ=255)
                                                                                                                                          gCAAAATTTGCGAATCTAAACCGTCCCCAATTGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAg    <  7:228055/149‑1 (MQ=255)
                                                                                                                                          gCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGaaaa     <  2:308277/148‑1 (MQ=255)
                                                                                                                                           cAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTg                <  7:188570/136‑1 (MQ=255)
                                                                                                                                           cAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTg                >  8:188570/1‑136 (MQ=255)
                                                                                                                                           cAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGc   <  2:137264/149‑1 (MQ=255)
                                                                                                                                            aaaaTTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGcc  >  1:221620/1‑149 (MQ=255)
                                                                                                                                            aaaaTTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGc   <  1:285487/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ACATAATGCCAACTTCATCGACTTCTTCGAAGAAGGTCTCAAGCATATTGAGCCGAGCAAATCCGAAGTAAACACCCTGTATCGAGAGTATCACTCATACTTGACAGTTGGCAGCGCGAGCAATTCAACAGTTGTTCGCAAAATTTGCGAAGCTAAACCGTCCCCAATGGTAAGCACGGCATCGTCAAAAAAATATCACGCTCAAGTCCAGCCTTTCCTGAAGAGCCTTGGTAATTACTTCCAGAAATACAAAGACTACTGCGAGAACGGGCTGGTAATTGAAAAGCC  >  NC_002947/5002386‑5002673

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