Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 201,385 T→C 7.3% N2297N (AAT→AAC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947201,3850TC 7.3% 195.7 / 4.2 68N2297N (AAT→AAC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (38/25);  new base C (2/3);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAA  >  NC_002947/201239‑201531
                                                                                                                                                  |                                                                                                                                                   
aacctgGACAAGAACACGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGTCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAAt                                                                                                                                                     <  7:143219/147‑1 (MQ=255)
  cctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGc                                                                                                                                                 >  5:106236/1‑149 (MQ=255)
   ctgGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGcc                                                                                                                                                <  6:216636/149‑1 (MQ=255)
        cAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                                           >  7:202846/1‑149 (MQ=255)
        cAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                                           <  7:340918/149‑1 (MQ=255)
         aaGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                                           >  4:93368/1‑148 (MQ=255)
         aaGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                                           <  3:93368/148‑1 (MQ=255)
           gACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGcc                                                                                                                                         >  2:172719/1‑148 (MQ=255)
              cccGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGc                                                                                                                                              >  1:328624/1‑140 (MQ=255)
              cccGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGc                                                                                                                                              <  2:328624/140‑1 (MQ=255)
                cGGGCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGAC                                                                                                                                   <  8:276955‑M1/149‑1 (MQ=255)
                   tCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCT                                                                                                                                 >  2:69318‑M1/1‑148 (MQ=255)
                      gCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATc                                                                                                                             <  3:295038/149‑1 (MQ=255)
                        aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                                           <  7:288733/133‑1 (MQ=255)
                        aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTgg                                                                                                                                           >  8:288733/1‑133 (MQ=255)
                        aCTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGt                                                                                                                           >  6:118323/1‑149 (MQ=255)
                         cTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGCTCGt                                                                                                                           >  6:254275/1‑148 (MQ=255)
                            ccGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGAcc                                                                                                                       >  2:81838/1‑149 (MQ=255)
                                          ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGc                                                                                                         >  5:264085/1‑149 (MQ=255)
                                          ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGc                                                                                                         >  8:207591/1‑149 (MQ=255)
                                            aGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCt                                                                                                        >  2:1205/1‑148 (MQ=255)
                                                aCTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAgg                                                                                                    >  4:165357/1‑148 (MQ=255)
                                                  tCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCa                                                                                                 >  7:132124/1‑149 (MQ=255)
                                                    atctcaacgAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtg                                                                                               <  7:296745‑M1/145‑34 (MQ=255)
                                                     gggcAACGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtgg                                                                                              <  1:326165‑M1/148‑35 (MQ=255)
                                                      ggCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCAt                                                                                            <  4:243247/149‑1 (MQ=255)
                                                         aaCGAAGGCGACCTGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtg                                                                                               >  8:296745‑M1/1‑111 (MQ=255)
                                                         aaCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATc                                                                                           >  5:263275/1‑147 (MQ=255)
                                                         aaCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCaa                                                                                         <  6:106236/149‑1 (MQ=255)
                                                         aaCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCaa                                                                                         <  5:118323/149‑1 (MQ=255)
                                                            gAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGc                                                                                      >  1:345611/1‑149 (MQ=255)
                                                                    cATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGaaa                                                                              <  1:69318/149‑1 (MQ=255)
                                                                         tCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACACCGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgc                                                                         >  6:238995/1‑149 (MQ=255)
                                                                         tCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑c                                                                                               >  3:233601/1‑127 (MQ=255)
                                                                         tCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑c                                                                                               <  4:233601/127‑1 (MQ=255)
                                                                         tCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgc                                                                         <  1:1205/149‑1 (MQ=255)
                                                                          aaaGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTTAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgcg                                                                        >  2:223307/2‑149 (MQ=255)
                                                                           aaGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgcgc                                                                       >  8:39934/1‑149 (MQ=255)
                                                                            aGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgcgcg                                                                      >  2:232605/1‑149 (MQ=255)
                                                                             ggTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAgcgcgc                                                                     >  7:226537/1‑149 (MQ=255)
                                                                              gTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACAGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGcc                                                                    >  4:318379/1‑149 (MQ=255)
                                                                               tCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACACCGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCg                                                                   >  6:348969/1‑149 (MQ=255)
                                                                                cACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCACGTCA‑CCATCAAGGCCCGTTAAACCAGCGCGCCg                                                                   >  2:90534/1‑148 (MQ=255)
                                                                                 acaaTCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACAGACCTTCAATGTGCACGTCAATCAGCCGCTGGACCACGACCTGATCTTGACCCTGAGCAACAACGATCAGGTAA‑CCATCAAGGCCGGTGAAACCAGCGCGCAGTa                                                                 >  6:270736/4‑149 (MQ=14)
                                                                                    aTCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtgg                                                                                              <  8:319322‑M1/118‑35 (MQ=255)
                                                                                    aTCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACCTTCACCGTGCACGTCAACCAGCCGCTGGCCCACGACCTGggcgttgttgctcagggtcacgaccaggtcgtgg                                                                                              >  7:319322‑M1/1‑84 (MQ=255)
                                                                                       aCGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTTAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACAccc                                                            <  1:223307/148‑1 (MQ=255)
                                                                                       aCGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCa                                                           <  1:172719/149‑1 (MQ=255)
                                                                                              aCCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGAcgcg                                                    <  8:132124/149‑1 (MQ=255)
                                                                                               ccAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGAcgcg                                                    >  4:296137/1‑148 (MQ=255)
                                                                                                  gACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCACCCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCaa                                                 <  1:81838/149‑1 (MQ=255)
                                                                                                  gACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAg                                                >  6:209362/1‑149 (MQ=255)
                                                                                                   aCCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAgg                                               <  5:209362/149‑1 (MQ=255)
                                                                                                   aCCTCGGTAGACGAGAACGTCAAACCGACCTTCACTGTTCACTTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCTCGCCGTACACCCACGACGCGCAAgg                                               <  5:270736/149‑1 (MQ=255)
                                                                                                     cTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGc                                              >  3:321652/1‑148 (MQ=255)
                                                                                                        ggTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACCCCCACGACGCGCAAGGCGATg                                          >  8:221258/1‑149 (MQ=255)
                                                                                                            gCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGt                                      <  2:345611/149‑1 (MQ=39)
                                                                                                             ccGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACACTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTc                                     >  7:258567/1‑149 (MQ=255)
                                                                                                              cGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCt                                    <  8:226537/149‑1 (MQ=39)
                                                                                                                 gaACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATc                                 <  1:46161/149‑1 (MQ=40)
                                                                                                                     gTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCCTCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGa                             >  3:98365/1‑149 (MQ=255)
                                                                                                                     gTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACTACGCGCAACGCCATCACGTCTATCAGGa                             >  4:54486/1‑149 (MQ=255)
                                                                                                                          aCCGACCTTCACTCTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTg                        <  6:311777/149‑1 (MQ=255)
                                                                                                                           ccGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTgg                       >  8:141999/1‑149 (MQ=255)
                                                                                                                              aCCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCa                    <  1:232605/149‑1 (MQ=255)
                                                                                                                                  tCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGAt                 >  5:117784/1‑148 (MQ=255)
                                                                                                                                  tCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGAt                 >  3:78609/1‑148 (MQ=255)
                                                                                                                                          cACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACACTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTgg         <  6:117784/148‑1 (MQ=255)
                                                                                                                                                aaTCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCaa  >  5:237108/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AACCTGGACAAGACCCCGGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAATCAGCCGCTGGCCCACGACCTGATCGTGACCCTGAGCAACAACGCTCAGGTCA‑CCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTATCAGGACGCTGGCCAGATCAGCTTGGGCATCAA  >  NC_002947/201239‑201531

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