Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,997 T→A 14.2% K246M (AAG→ATG) ‡ PP_4563 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,9970TA14.2% 175.8 / 16.6 70K246M (AAG→ATG) ‡PP_4563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/26);  new base A (3/7);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCA  >  NC_002947/5182851‑5183144
                                                                                                                                                  |                                                                                                                                                     
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                                                                                                                                            >  2:366241/1‑143 (MQ=255)
               aCTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                            <  1:366241/143‑1 (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                                                                                            >  3:95368/1‑149 (MQ=255)
                                                        tcGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTc                                                                                            >  7:257690/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                                                                                                                               >  6:3378/1‑76 (MQ=255)
                                                                                              gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑ctc                                                                                                                               <  5:3378/76‑1 (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                                 <  3:124079/144‑1 (MQ=255)
                                                                                                                       aaGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGc                                 >  4:124079/1‑144 (MQ=255)
                                                                                                                         gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCt                            >  1:283353/1‑147 (MQ=255)
                                                                                                                          ttGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAg                         <  7:187869/149‑1 (MQ=255)
                                                                                                                              ccTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAg                               >  2:237450/1‑139 (MQ=255)
                                                                                                                              ccTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAg                               <  1:237450/139‑1 (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)
                                                                                                                                                 cTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCa  <  1:21347/149‑1 (MQ=255)
                                                                                                                                                 cTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGAGCGCCAT‑GATCTGCTCTTGCAGGCTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGATCAGTGAAATCTTGCCGTCATGGCGCa  >  7:336004/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTGGCGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCA‑CTCTTGTCGCACACCCCCAACACCCCGCCTTCGGCTTCCAGGCGCGCCAT‑GATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCA  >  NC_002947/5182851‑5183144

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