Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 807,908 G→A 5.3% Y42Y (TAC→TAT)  PP_0694 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947807,9080GA 5.3% 291.2 / 3.3 93Y42Y (TAC→TAT) PP_0694hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/47);  new base A (3/2);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

AGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGTTGGGTTCG  >  NC_002947/807761‑808056
                                                                                                                                                     |                                                                                                                                                    
aGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCgg                                                                                                                                                      >  6:4389/1‑148 (MQ=255)
 gCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTa                                                                                                                                                    >  2:310399/1‑149 (MQ=255)
   aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGTACAGG‑GTGCGGTACt                                                                                                                                                  <  7:349652/149‑1 (MQ=255)
   aaGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGCTCATGTCCGGCACGATGTAGTAGAAACCACGGGTATCGCACACAAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACt                                                                                                                                                  >  6:366928/1‑149 (MQ=255)
     ggTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTgg                                                                                                                                                >  3:319049/1‑149 (MQ=255)
        gCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTcc                                                                                                                                             >  1:151453/1‑149 (MQ=255)
         cTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                                            <  4:61754/149‑1 (MQ=255)
          tCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACAACAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTg                                                                                                                                           >  5:183270/1‑149 (MQ=255)
            gACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                                            <  7:252142/146‑1 (MQ=255)
            gACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCt                                                                                                                                            >  8:252142/1‑146 (MQ=255)
             aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTgg                                                                                                                                          >  3:32476/1‑147 (MQ=255)
             aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                                                 <  3:223967/140‑3 (MQ=25)
             aCCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATAccc                                                                                                                                                 >  4:223967/1‑138 (MQ=25)
                aGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGtt                                                                                                                                     <  4:32476/149‑1 (MQ=255)
                   aTGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                                                  <  2:116847/149‑1 (MQ=255)
                    tGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                                                  >  2:336574/1‑148 (MQ=255)
                    tGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACAAGTTCG‑AAGCCCAGGCCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGa                                                                                                                                  >  5:194476/1‑148 (MQ=255)
                            aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTggg                                                                                                                          <  4:319049/148‑1 (MQ=255)
                            aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                                                         >  1:355153/1‑149 (MQ=255)
                            aaCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                                                         <  5:170569/149‑1 (MQ=255)
                             aCAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGc                                                                                                                         <  1:171383/148‑1 (MQ=255)
                                gCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTAc                                                                                                                     <  6:265899/149‑1 (MQ=255)
                                 cGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTa                                                                                                                      >  3:142164/1‑147 (MQ=255)
                                       gcgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCGTCAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCa                                                                                                               <  1:310399/149‑1 (MQ=255)
                                        cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                                             <  8:76081/149‑1 (MQ=255)
                                        cgcTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAgg                                                                                                             <  2:151453/149‑1 (MQ=255)
                                           tGCGCGGTCTCGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCGGGGTc                                                                                                          >  3:253965/1‑149 (MQ=255)
                                             cgcgGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTcgc                                                                                                        <  2:84123/149‑1 (MQ=255)
                                              ggggTGTTGTTGACCTGCGCGGGGGTCATGTCCGGCACGAAGTATTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTcgcg                                                                                                       <  6:194476/147‑1 (MQ=255)
                                                ggTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCAGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGcc                                                                                                     >  2:143636/1‑149 (MQ=255)
                                                    ttgttgACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                                 >  7:10279/1‑149 (MQ=255)
                                                    gacatgACCGCCGCGCAGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                                 >  4:93474/5‑149 (MQ=255)
                                                    gacatgACCGCCGCGCAGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCt                                                                                                 >  2:221444/5‑149 (MQ=255)
                                                          aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTg                                                                                           <  1:336574/149‑1 (MQ=255)
                                                          aCCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTAG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGCAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTg                                                                                           >  6:289396/1‑149 (MQ=255)
                                                            cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                                          >  3:35609/1‑148 (MQ=255)
                                                            cTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTgg                                                                                          >  8:67872/1‑148 (MQ=255)
                                                               cgcgGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAg                                                                                       <  5:66872/148‑1 (MQ=255)
                                                                 cggcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCt                                                                                    >  8:270216/1‑149 (MQ=255)
                                                                   gcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                                  <  7:67872/149‑1 (MQ=255)
                                                                   gcggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                                  <  2:138920/149‑1 (MQ=255)
                                                                    cggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCGGGGTCGCGCCGGCTGACGTGGAAGCCTTTg                                                                                 <  4:253965/149‑1 (MQ=255)
                                                                    cggTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GGGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTg                                                                                 <  4:51058/149‑1 (MQ=255)
                                                                       tCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGaa                                                                               >  2:60983/1‑148 (MQ=255)
                                                                          tGTCCGGAACTAAGTAGTAGAAACCCCGGGTATCGCACACCCGTTCG‑AAGCCCCGGGCCTTGAACAGG‑GGGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                                           <  8:28121/149‑1 (MQ=255)
                                                                           gTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGc                                                                           >  4:97630/1‑148 (MQ=255)
                                                                            tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                                  <  5:87628/140‑1 (MQ=255)
                                                                            tCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCttt                                                                                  >  6:87628/1‑140 (MQ=255)
                                                                              cGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTcgc                                                                       <  1:60983/149‑1 (MQ=255)
                                                                               ggCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTTACGTGGAAGCCTTTGAACAGCTcgcg                                                                      <  6:183270/149‑1 (MQ=255)
                                                                                 cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                                     >  2:286983/1‑148 (MQ=255)
                                                                                 cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                                     >  7:284163/1‑148 (MQ=255)
                                                                                 cACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGa                                                                    <  5:4389/149‑1 (MQ=255)
                                                                                    gAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCgg                                                                     >  2:141098/1‑145 (MQ=255)
                                                                                         agtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCgg                                                             <  4:35609/148‑1 (MQ=255)
                                                                                          gtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcg                                                           >  7:220357/1‑149 (MQ=255)
                                                                                          gtagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgc                                                            <  7:270216/148‑1 (MQ=255)
                                                                                           tagAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                                          >  1:149557/1‑149 (MQ=255)
                                                                                             gAAACCACGGGTATCGCACACCGGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTGGCGGAAGATCGgcgc                                                          >  7:59424/1‑147 (MQ=255)
                                                                                             gAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGgcgc                                                          >  6:33613/1‑147 (MQ=255)
                                                                                                aCCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCt                                                     >  2:342658/1‑149 (MQ=255)
                                                                                                 ccACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCTGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTg                                                    <  1:286983/149‑1 (MQ=255)
                                                                                                  cACGGGCATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGGCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTgg                                                   <  8:158044/149‑1 (MQ=255)
                                                                                                     cagTACCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                                                 >  1:56034/3‑148 (MQ=25)
                                                                                                      agTACCGCAC‑CCTGTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGa                                                 >  3:95580/2‑147 (MQ=25)
                                                                                                        tATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGt                                             >  5:21947/1‑149 (MQ=255)
                                                                                                          tCGCACACCATAGCG‑TACCCCAGGGCCTTGCACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACACCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGtt                                            <  8:165315/148‑1 (MQ=255)
                                                                                                          tCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGTCAGtt                                            <  8:284163/148‑1 (MQ=255)
                                                                                                            gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTc                                           >  2:82327/1‑147 (MQ=255)
                                                                                                            gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTc                                           <  1:82327/147‑1 (MQ=255)
                                                                                                            gCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCgg                                         <  1:141098/149‑1 (MQ=255)
                                                                                                             cacaCCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa                                        <  3:97630/149‑1 (MQ=255)
                                                                                                                ccctgTTCA‑AGGCCCTGGGCTTCGAACTGGTGTGCGATACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGaaa                                      <  4:95580/145‑1 (MQ=37)
                                                                                                                acCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGTAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAg                                     <  8:220357/149‑1 (MQ=255)
                                                                                                                acCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGa                                        >  3:37701/1‑146 (MQ=255)
                                                                                                                 ccAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGa                                    <  2:355153/149‑1 (MQ=255)
                                                                                                                  cAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTCAACAGCTCGCGGAAGATCGGCGCGAGCTGGTACAGTTCGGAAAGa                                    >  4:49877/1‑148 (MQ=255)
                                                                                                                     ttCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCaa                                <  5:33613/149‑1 (MQ=255)
                                                                                                                        g‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCAGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGac                             >  3:128662/1‑148 (MQ=255)
                                                                                                                          aaGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATg                            >  8:50674/1‑149 (MQ=255)
                                                                                                                              ccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAg                               <  7:190377/142‑1 (MQ=255)
                                                                                                                              ccAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAg                               >  8:190377/1‑142 (MQ=255)
                                                                                                                                     cTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGTGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgggctggg                 >  2:28929/1‑149 (MQ=255)
                                                                                                                                     cTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgggctggg                 <  4:37701/149‑1 (MQ=255)
                                                                                                                                     cTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCAtgggctggg                 >  4:101239/1‑149 (MQ=255)
                                                                                                                                             gg‑gTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGtt         <  8:10279/149‑1 (MQ=255)
                                                                                                                                              g‑gTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTGAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACCGCTCGCGGATGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGTTg        <  3:49877/149‑1 (MQ=255)
                                                                                                                                              g‑gTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGTTg        <  2:149557/149‑1 (MQ=255)
                                                                                                                                                gTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGTTgg       <  4:142164/149‑1 (MQ=255)
                                                                                                                                                gTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGTTgg       <  6:21947/149‑1 (MQ=255)
                                                                                                                                                gTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGTTgg       <  3:93474/149‑1 (MQ=255)
                                                                                                                                                  gCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGtt         >  4:300298/1‑145 (MQ=255)
                                                                                                                                                  gCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGtt         <  3:300298/145‑1 (MQ=255)
                                                                                                                                                   cGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGTTGGGtt    <  8:59424/149‑1 (MQ=255)
                                                                                                                                                     gTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGTTGGGTTCg  >  4:368950/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCAAGGTGCTCGACCAGGATGAAGGTGAACAGCGACAGGCGCTGCGCGGTCTTGTTGACCTGCGCGGCGGTCATGTCCGGCACGAAGTAGTAGAAACCACGGGTATCGCACACCAGTTCG‑AAGCCCAGGGCCTTGAACAGG‑GTGCGGTACTGGTCCTGGAAGTTCGACAGCTGGGCGTACAGCTCAGGGTCGCGCCGGCTGACGTGGAAGCCTTTGAACAGCTCGCGGAAGATCGGCGCGAGCTGGGACAGTTCGGAAAGATCAAGATGCATGGGCTGGGCTCGCGTTGGGTTCG  >  NC_002947/807761‑808056

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