Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,770,438 A→T 6.0% L145H (CTT→CAT)  PP_5060 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,770,4380AT 6.0% 201.9 / 2.7 67L145H (CTT→CAT) PP_5060hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base T (2/2);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

ACGCAGCAGCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGG  >  NC_002947/5770290‑5770584
                                                                                                                                                     |                                                                                                                                                     
aCGCAGCAGCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCaaa                                                                                                                                                       <  7:49457/149‑1 (MQ=255)
  gcagcaGCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAag                                                                                                                                                      >  5:135967/1‑148 (MQ=255)
        gcTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCAtt                                                                                                                                               <  4:289885/149‑1 (MQ=255)
            gCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTgg                                                                                                                                             <  6:108863/147‑1 (MQ=255)
             cGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGgcag                                                                                                                                          <  6:135967/149‑1 (MQ=255)
              gttgttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcact                                                                                                                                         >  6:31496/1‑147 (MQ=255)
              gttgttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGgcagc                                                                                                                                         <  7:282071/149‑1 (MQ=255)
               ttgttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGgcagca                                                                                                                                        <  1:168952/149‑1 (MQ=255)
                  ttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcacttcca                                                                                                                                     <  5:31496/149‑7 (MQ=255)
                  ttCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAAcg                                                                                                                                     <  6:256278/149‑1 (MQ=255)
                      tCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGtt                                                                                                                                 <  3:146839/149‑1 (MQ=255)
                         tCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAAcgc                                                                                                                                    >  8:43211/1‑143 (MQ=255)
                         tCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAAcgc                                                                                                                                    <  7:43211/143‑1 (MQ=255)
                          ccAGCGCTTCATCCTTCAGGCCATCGAGCGAGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAagcagc                                                                                                                                                  >  1:179363/1‑129 (MQ=255)
                          ccAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAagcagc                                                                                                                                                  <  2:179363/128‑1 (MQ=255)
                           cAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcacttcc                                                                                                                                      <  8:237214/139‑6 (MQ=255)
                           cAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTACAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGAAGCATTGGCAGAAACGAGTTCAGGc                                                                                                                            >  3:2844/1‑149 (MQ=255)
                           cAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTAGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAATGCTGCTTTGgcacttcc                                                                                                                                      >  7:237214/1‑134 (MQ=255)
                            aGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGc                                                                                                                            >  8:35310/1‑148 (MQ=255)
                              cgcTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCAt                                                                                                                         <  4:141880/149‑1 (MQ=255)
                                       cTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCAt                                                                                                                         >  1:215281/1‑140 (MQ=255)
                                       cTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCAt                                                                                                                         <  2:215281/140‑1 (MQ=255)
                                        ttCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGgcg                                                                                                              >  7:113989/1‑149 (MQ=255)
                                        ttCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGgcg                                                                                                              >  1:16130/1‑149 (MQ=255)
                                           aGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCAATCG‑GTCAGgcgcgc                                                                                                           >  2:165585/1‑149 (MQ=255)
                                            ggCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCACCGCGTTCCGGCCATCG‑GTCAGGcgcccg                                                                                                          >  8:120153/1‑149 (MQ=255)
                                                aTCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCt                                                                                                      <  4:211562/149‑1 (MQ=255)
                                                   gAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCgg                                                                                                   >  2:70554/1‑148 (MQ=255)
                                                   gAGC‑AGGGTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGGTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCAGCG‑GTCAGGCGCGCGGGCTCg                                                                                                    <  7:120153/148‑1 (MQ=255)
                                                          gTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCGGTCGCCGTTCCAATAGTCGATCTGCGCTGGCTGGATGTGCCAAAGCAGCATTGGCAGCAGCGCGTTCAGGCCTTCG‑GTCGGGCGCGCGGGCTCGATCatat                                                                                             >  8:266705/1‑149 (MQ=255)
                                                          gTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAgg                                                                                                                             >  8:43487/1‑118 (MQ=255)
                                                          gTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAt                                                                                             >  2:81926/1‑149 (MQ=255)
                                                              gcgcgcAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGcc                                                                                       <  8:113989/149‑1 (MQ=255)
                                                                gcgcAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGtt                                                                                                                                 >  2:86048/1‑108 (MQ=255)
                                                                gcgcAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGtt                                                                                                                                 <  1:86048/108‑1 (MQ=255)
                                                                    aacaacTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAagc                                                                                                                                                     <  1:116940/84‑1 (MQ=255)
                                                                    aacaacTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAagc                                                                                                                                                     >  2:116940/1‑84 (MQ=255)
                                                                    aacaacTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGc                                                                                 >  6:141840/1‑149 (MQ=255)
                                                                         ctgctgCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGttt                                                                            <  1:81926/149‑1 (MQ=255)
                                                                              gcAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAg                                                                       >  2:249496/1‑149 (MQ=255)
                                                                               cAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGaa                                                                        <  5:141840/147‑1 (MQ=255)
                                                                                  ccACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGATTCGAAGctc                                                                    >  7:81562/1‑148 (MQ=255)
                                                                                      aGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCAAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGctcgactc                                                               >  8:211703/1‑149 (MQ=255)
                                                                                           ccGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATgg                                                          >  1:152707/1‑149 (MQ=255)
                                                                                             ggTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCACCCGTTTCGAAGCTCGACTCGATGGCg                                                        <  1:165585/149‑1 (MQ=255)
                                                                                                  aTCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGaaaaa                                                   <  7:158997/149‑1 (MQ=255)
                                                                                                     gTGCCGTCGCCGTTCCAATAGTCGCTCTGCGCTTGCTGGAAGTGCCATAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAgc                                                 <  4:28591/148‑1 (MQ=255)
                                                                                                      tGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAgcgc                                               >  4:266164/1‑149 (MQ=255)
                                                                                                           tCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCgggg                                          >  6:228048/1‑149 (MQ=255)
                                                                                                              ccGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCgg                                       >  4:240964/1‑149 (MQ=255)
                                                                                                              ccGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCgg                                       <  2:16130/149‑1 (MQ=255)
                                                                                                              ccGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCgg                                       <  3:240964/149‑1 (MQ=255)
                                                                                                               cGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCggg                                      <  3:128720/149‑1 (MQ=255)
                                                                                                                 ttCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGtt                                    <  5:228048/149‑1 (MQ=255)
                                                                                                                   ccAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGtt                                    >  3:187974/1‑147 (MQ=255)
                                                                                                                   ccAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGtg                                    >  3:194736/1‑146 (MQ=255)
                                                                                                                     aaTAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAgg                                >  6:224377/1‑149 (MQ=255)
                                                                                                                      aTAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGc                               <  4:187974/149‑1 (MQ=255)
                                                                                                                             aTCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCGAATCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAACGCGCCGGGGCGGGTTCAGGCTAAGGc                         >  1:238670/1‑149 (MQ=255)
                                                                                                                              tCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCgg                       >  7:28898/1‑149 (MQ=255)
                                                                                                                                  ggcTTTCTGGAAGTGGCAAATCATCATTTGGATCAAAGCGTTCATGCCATTG‑GTGATGGGCGGGGGCTCGATCAAATTGCTCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTa                     <  1:249496/148‑1 (MQ=255)
                                                                                                                                         tGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCg                                                                          >  4:224481/1‑87 (MQ=255)
                                                                                                                                          ggAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTcgcccgcc           >  4:83185/1‑149 (MQ=255)
                                                                                                                                           gAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCccgccg          >  1:123511/1‑149 (MQ=255)
                                                                                                                                            aaGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATc                                                                                                 >  2:156071/1‑63 (MQ=255)
                                                                                                                                            aaGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATc                                                                                                 <  1:156071/63‑1 (MQ=255)
                                                                                                                                             aGTGCAAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCgg         <  4:2844/148‑1 (MQ=255)
                                                                                                                                               ggCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGc      <  3:1396/148‑1 (MQ=255)
                                                                                                                                                 ccAACGCAGCATTCGCAGCAACGCGTTCAGGCCATCG‑TTCAGGCGCGCGGGCACGATCAAAT‑‑CGCCCTCAGCCCGTTCGAAGCTCGACTCGATGGCGAAAAAGCCCCGTGGCGGGTTCAGGCTAAGGCGGTAGACGCCCGCCGGGGc      >  3:243899/1‑147 (MQ=255)
                                                                                                                                                 ccAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGc      >  4:219192/1‑147 (MQ=255)
                                                                                                                                                  cAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAg   <  7:35310/149‑1 (MQ=255)
                                                                                                                                                  cAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAg   >  2:187738/1‑149 (MQ=255)
                                                                                                                                                  cAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAg   <  2:152707/149‑1 (MQ=255)
                                                                                                                                                   aaaGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAgg  >  2:47060/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGCAGCAGCTTGCGTTGTTCGTCATCCAGCGCTTCATCCTTCAGGCCATCGAGC‑AGGTTGGCGCGCAACAACTGCTGCAACCACAGGTCCCGGTCGATCGTGCCGTCGCCGTTCCAATAGTCGATCTGCGCTTGCTGGAAGTGCCAAAGCAGCATTGGCAGCAACGCGTTCAGGCCATCG‑GTCAGGCGCGCGGGCTCGATCAAAT‑‑CGCCCTCAGCCGTTTCGAAGCTCGACTCGATGGCGAAAAAGCGCCGGGGCGGGTTCAGGCTAAGGCGGTAGTCGCCCGCCGGGGCGAGG  >  NC_002947/5770290‑5770584

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