Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,992 A→T 14.9% F248I (TTC→ATC)  PP_4563 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,9920AT14.9% 156.9 / 18.9 67F248I (TTC→ATC) PP_4563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/24);  new base T (3/7);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGG  >  NC_002947/5182851‑5183140
                                                                                                                                             |                                                                                                                                                      
gtgtCGCGCAGCTGGACTTCCGCAACTTCGTACAGGAGTAAGTCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTg                                                                                                                                                 <  3:201420/149‑1 (MQ=255)
    cGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                             <  2:276557/149‑1 (MQ=255)
         agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                        <  7:246038/149‑1 (MQ=255)
          ggggtaCTTCCGCAACTTCTTACAGCCGTAAGGCCGTAAGGC‑TTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                        <  1:132889/142‑1 (MQ=255)
           cTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgt                                                                                                                                      <  2:370074/149‑2 (MQ=255)
              gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtca                                                                                                                                    <  3:207730/148‑4 (MQ=255)
              gACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAACgtc                                                                                                                                     <  4:64707/147‑3 (MQ=255)
               aCTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                        <  1:366241/143‑1 (MQ=255)
               aCTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                        >  2:366241/1‑143 (MQ=255)
               aCTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAggg                                                                                                                                  <  6:362455/149‑1 (MQ=255)
                    cGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑ct                                                                                                                            >  2:67899/1‑149 (MQ=255)
                            cGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTc                                                                                                                      <  3:317918/147‑1 (MQ=255)
                                   gCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACAcccc                                                                                                             >  1:332797/1‑149 (MQ=255)
                                         ggCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAacac                                                                                                       >  5:206240/1‑149 (MQ=255)
                                           ccGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACAccc                                                                                                     >  8:146582/1‑149 (MQ=255)
                                                aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCt                                                                                                >  2:66381/1‑149 (MQ=255)
                                                  gCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTc                                                                                              <  8:242853/149‑1 (MQ=255)
                                                    tttCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCgg                                                                                            >  1:73111/1‑149 (MQ=255)
                                                      tctcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAAACCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTCCTGATCAACGAAGGGTCA‑ATCTTGTCGAACACCCCCAACACCACGCATTctcat                                                                                          >  7:79413/1‑145 (MQ=255)
                                                        tcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTc                                                                                        >  7:257690/1‑149 (MQ=255)
                                                        tcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTc                                                                                        >  3:95368/1‑149 (MQ=255)
                                                           aTGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAg                                                                                     <  7:146582/149‑1 (MQ=255)
                                                                   gTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCa                                                                             >  2:81930/1‑149 (MQ=255)
                                                                    tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTACCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAt                                                                            >  1:286638/1‑149 (MQ=255)
                                                                      aTGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑Ga                                                                         <  2:195354/149‑1 (MQ=255)
                                                                               gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCtt                                                                >  6:228066/1‑149 (MQ=255)
                                                                               gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCtt                                                                <  6:206240/149‑1 (MQ=255)
                                                                                 tCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTg                                                               >  4:271161/1‑148 (MQ=255)
                                                                                   aGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCa                                                             <  2:286638/148‑1 (MQ=255)
                                                                                      cccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGgt                                                         <  8:257690/149‑1 (MQ=255)
                                                                                       ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGgtg                                                        >  6:204704/1‑149 (MQ=255)
                                                                                           cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTcc                                                     <  5:228066/148‑1 (MQ=255)
                                                                                           cTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTcc                                                     <  1:66381/148‑1 (MQ=255)
                                                                                            ttGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGc                                                   >  2:185554/1‑149 (MQ=255)
                                                                                              gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑ctc                                                                                                                           <  5:3378/76‑1 (MQ=255)
                                                                                              gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑ctc                                                                                                                           >  6:3378/1‑76 (MQ=255)
                                                                                               aacaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAg                                                >  4:186195/1‑149 (MQ=255)
                                                                                                acaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCTCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAgg                                               <  1:67899/149‑1 (MQ=255)
                                                                                                acaacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAgg                                               <  4:95368/149‑1 (MQ=255)
                                                                                                 caacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAggg                                              <  2:73111/149‑1 (MQ=255)
                                                                                                  aacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGggc                                             >  5:291708/1‑149 (MQ=255)
                                                                                                       gCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGgcgc                                        >  8:187869/1‑149 (MQ=255)
                                                                                                       gCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGgcgc                                        >  6:151458/1‑149 (MQ=255)
                                                                                                           aTGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAc                                    >  4:258085/1‑149 (MQ=255)
                                                                                                            tGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAcc                                   <  5:151458/149‑1 (MQ=255)
                                                                                                                cttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTc                                                                                              <  3:90768/87‑1 (MQ=255)
                                                                                                                cttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTc                                                                                              >  4:90768/1‑87 (MQ=255)
                                                                                                                 ttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTg                              >  2:60499/1‑149 (MQ=255)
                                                                                                                      gAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCATCGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGc                             <  1:81930/146‑1 (MQ=255)
                                                                                                                      gAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGc                         >  3:7412/1‑149 (MQ=255)
                                                                                                                       aaGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGc                             >  4:124079/1‑144 (MQ=255)
                                                                                                                       aaGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGc                             <  3:124079/144‑1 (MQ=255)
                                                                                                                         gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCt                        >  1:283353/1‑147 (MQ=255)
                                                                                                                          ttGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAg                     <  7:187869/149‑1 (MQ=255)
                                                                                                                              ccTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAg                           <  1:237450/139‑1 (MQ=255)
                                                                                                                              ccTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAg                           >  2:237450/1‑139 (MQ=255)
                                                                                                                               cTTGCTGACGTTGAACAGCTTGCTGATCAACGCTGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGaaa                >  8:250783/1‑149 (MQ=255)
                                                                                                                                 tGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATc              >  2:21347/1‑149 (MQ=255)
                                                                                                                                   cctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCAGCTCTTGTCGCACCCCCCCAACCCCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCt             >  1:59912/5‑149 (MQ=255)
                                                                                                                                   cTGACGTTGTACAGATTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCTGGCTCGCGAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCTGCGCCTACCGGCTGCAGGCTCAGTGAAATCtt            <  8:79413/149‑1 (MQ=255)
                                                                                                                                    tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGTCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg           >  6:44524/1‑149 (MQ=255)
                                                                                                                                    tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg           >  6:374047/1‑149 (MQ=255)
                                                                                                                                    ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg           >  2:361452/4‑149 (MQ=255)
                                                                                                                                      aCGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg           <  2:332797/147‑1 (MQ=255)
                                                                                                                                       cGTTGATCAGCAAGCTGTTCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGcc         <  1:361452/148‑1 (MQ=255)
                                                                                                                                           gAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCt    >  2:364939/1‑149 (MQ=255)
                                                                                                                                           gAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCt    >  1:182824/1‑149 (MQ=255)
                                                                                                                                             aCAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTgg  <  2:283353/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGG  >  NC_002947/5182851‑5183140

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