Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,015,088 C→T 7.0% Q65* (CAA→TAA)  pqiC → DUF330 family putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,015,0880CT 7.0% 227.9 / 6.3 71Q65* (CAA→TAA) pqiCDUF330 family putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (38/28);  new base T (2/3);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

CAGCTCCGGCGAAATTAATAAAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGGGAAGCGCCCAG  >  NC_000913/1014943‑1015234
                                                                                                                                                 |                                                                                                                                                  
cAGCTCCGGCGAAATTAATAAAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCaaa                                                                                                                                                 <  4:221569/149‑1 (MQ=255)
     ccGGCGAAATTAATAAAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTGATAAACCaca                                                                                                                                            <  1:413668/149‑3 (MQ=255)
     ccGGCGAAATTAATAAAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGCGAATGGTGTGGTTTATCAAACCAGt                                                                                                                                            <  4:222666/149‑1 (MQ=255)
      cGGCGAAATTAATAAAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTg                                                                                                                                           >  1:60934/1‑149 (MQ=255)
        gCGAAATTAATAAAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGAt                                                                                                                                         <  4:444269/149‑1 (MQ=255)
                  taaAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGt                                                                                                                                  <  3:222073/146‑1 (MQ=255)
                  taaAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGt                                                                                                                                  >  4:222073/1‑146 (MQ=255)
                  taaAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTAtg                                                                                                                               >  2:315245/1‑149 (MQ=255)
                  taaAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTAtg                                                                                                                               >  3:243884/1‑149 (MQ=255)
                  taaAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTAtg                                                                                                                               >  4:159179/1‑149 (MQ=255)
                    aaaaCTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTAtgtg                                                                                                                             >  1:426564/1‑149 (MQ=255)
                        cTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTg                                                                                                                         >  3:297358/1‑149 (MQ=255)
                           ttACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCca                                                                                                                      >  2:251299/1‑149 (MQ=255)
                           ttACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCca                                                                                                                      <  4:266894/149‑1 (MQ=255)
                              ccAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCcaaca                                                                                                                   >  3:185614/1‑149 (MQ=255)
                                  ttACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCcaacaacaa                                                                                                               <  4:243884/149‑1 (MQ=255)
                                        gtggtgCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAGCAACAACTtgtg                                                                                                         <  2:60934/149‑1 (MQ=255)
                                               agagCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGc                                                                                                  <  3:415523/149‑1 (MQ=255)
                                                    ggTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGtt                                                                                             >  2:91715/1‑149 (MQ=255)
                                                    ggTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCAAGCCCGtt                                                                                             >  2:410042/1‑149 (MQ=255)
                                                      tACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTgg                                                                                           >  2:268903/1‑149 (MQ=255)
                                                       acacAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGa                                                                                          <  4:297358/149‑1 (MQ=255)
                                                        cacaAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGAt                                                                                         >  1:315427/1‑149 (MQ=255)
                                                         acaAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATc                                                                                        >  1:61528/1‑149 (MQ=255)
                                                           aaaGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCaa                                                                                      >  2:188945/1‑149 (MQ=255)
                                                                  cGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTgc                                                                               <  1:191655/149‑1 (MQ=255)
                                                                  cGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTgc                                                                               >  2:390569/1‑149 (MQ=255)
                                                                    ccagccagGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGAt                                                                                         <  1:240354/137‑1 (MQ=255)
                                                                    ccagccagGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGAt                                                                                         >  2:240354/1‑137 (MQ=255)
                                                                       gccagGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTg                                                                                            >  1:159478/1‑131 (MQ=255)
                                                                       gccagGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTg                                                                                            <  2:159478/131‑1 (MQ=255)
                                                                       gccagGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAac                                                                          >  3:414935/1‑149 (MQ=255)
                                                                       gccagGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAac                                                                          <  2:222736/149‑1 (MQ=255)
                                                                          agGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACAcc                                                                       <  2:62646/149‑1 (MQ=255)
                                                                          agGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACAcc                                                                       >  2:420541/1‑149 (MQ=255)
                                                                           gggCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACAccc                                                                      >  1:165620/1‑149 (MQ=255)
                                                                               aaTCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGt                                                                  <  3:159179/149‑1 (MQ=255)
                                                                                aTCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGtt                                                                 <  1:420541/149‑1 (MQ=255)
                                                                                       gTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAAcc                                                          <  1:293936/149‑1 (MQ=255)
                                                                                            gggTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGt                                                     >  1:385113/1‑149 (MQ=255)
                                                                                              gTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTAc                                                   >  4:194832/1‑149 (MQ=255)
                                                                                                  agCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCaa                                               >  4:32154/1‑149 (MQ=255)
                                                                                                       gTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGcc                                          >  1:332879/1‑149 (MQ=255)
                                                                                                        tCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGAAGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGccc                                         >  2:452380/1‑149 (MQ=255)
                                                                                                            tGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCt                                     >  3:312733/1‑149 (MQ=255)
                                                                                                              ttCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCt                                                         <  2:363075/127‑1 (MQ=255)
                                                                                                              ttCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCt                                                         >  1:363075/1‑127 (MQ=255)
                                                                                                                 cTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAAGCTGAg                                                      >  4:307476/1‑127 (MQ=255)
                                                                                                                 cTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGggtg                                <  2:165620/149‑1 (MQ=255)
                                                                                                                 cTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGggtg                                >  3:272879/1‑149 (MQ=255)
                                                                                                                 cTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGggtg                                >  3:165879/1‑149 (MQ=255)
                                                                                                                 cTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTCGCCAACCTGAg                                                      <  3:307476/127‑1 (MQ=255)
                                                                                                                  tGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGggtgg                               >  2:414341/1‑149 (MQ=255)
                                                                                                                      tATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGGTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGc                           <  1:268903/149‑1 (MQ=255)
                                                                                                                          tGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTcc                       >  1:133140/1‑149 (MQ=255)
                                                                                                                               gggAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTgc                                                                               >  3:185133/1‑88 (MQ=255)
                                                                                                                               gggAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTgc                                                                               <  4:185133/88‑1 (MQ=255)
                                                                                                                               gggAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGcc                  <  4:272879/149‑1 (MQ=255)
                                                                                                                                    tGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTgg             <  3:196503/149‑1 (MQ=255)
                                                                                                                                      gtgtGGTTTATTAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGtt                                                                                 >  1:401978/1‑79 (MQ=255)
                                                                                                                                      gtgtGGTTTATTAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGtt                                                                                 >  3:307258/1‑79 (MQ=255)
                                                                                                                                      gtgtGGTTTATTAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGtt                                                                                 <  2:401978/79‑1 (MQ=255)
                                                                                                                                      gtgtGGTTTATTAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGtt                                                                                 <  4:307258/79‑1 (MQ=255)
                                                                                                                                       tgtgGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGTGaa          <  1:410042/149‑1 (MQ=255)
                                                                                                                                       tgtgGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGGGaa          <  1:389952/149‑1 (MQ=255)
                                                                                                                                       tgtgGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGGGaa          >  4:22768/1‑149 (MQ=255)
                                                                                                                                           gTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGGGAAgcgc      >  1:158531/1‑149 (MQ=255)
                                                                                                                                            tttATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGGGAAGCGcc     >  2:59652/1‑149 (MQ=255)
                                                                                                                                            tttATCAAACCAGTGATGTGAAGTATGTGATTGCAAACAACATCTTGTGGGCCAGCCCGTTGGATCTCCAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGGGAAGCGcc     <  1:91715/149‑1 (MQ=255)
                                                                                                                                              tATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGGGAAGCGCCCa   <  3:22768/149‑1 (MQ=255)
                                                                                                                                               aTCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGGGAAGCGCCCAg  >  2:295617/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CAGCTCCGGCGAAATTAATAAAAACTATTACCAGTTACCTGTGGTGCAGAGCGGTACACAAAGTACCGCCAGCCAGGGCAATCGTCTGTTATGGGTAGAGCAGGTCACTGTTCCTGACTATCTGGCGGGGAATGGTGTGGTTTATCAAACCAGTGATGTGAAGTATGTGATTGCCAACAACAACTTGTGGGCCAGCCCGTTGGATCAACAGTTGCGCAACACCCTGGTTGCCAACCTGAGTACGCAACTGCCCGGCTGGGTGGTTGCCTCCCAGCCTCTGGGAAGCGCCCAG  >  NC_000913/1014943‑1015234

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