Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,109,246 C→T 7.3% A3T (GCA→ACA)  PP_0968 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,109,2460CT 7.3% 202.0 / 4.2 68A3T (GCA→ACA) PP_0968hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base T (3/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCTGGTGGAGTTGGCGGTAGAGCGCATGCGCG  >  NC_002947/1109099‑1109392
                                                                                                                                                   |                                                                                                                                                  
cGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGcc                                                                                                                                                   >  6:319380/1‑149 (MQ=255)
cGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGcc                                                                                                                                                   >  7:179975/1‑149 (MQ=255)
       cAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACt                                                                                                                                            >  7:298177/1‑149 (MQ=255)
        aaGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTg                                                                                                                                           >  3:30959/1‑149 (MQ=255)
            ttttAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAg                                                                                                                                       <  5:319380/149‑1 (MQ=255)
                 gggCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAAt                                                                                                                                   <  4:201063/148‑1 (MQ=255)
                   ccggcTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATAc                                                                                                                                 >  3:206445/2‑148 (MQ=255)
                      gcTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAgg                                                                                                                             >  6:153706/1‑149 (MQ=255)
                        tCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAggtg                                                                                                                           <  8:298177/149‑1 (MQ=255)
                           ggTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGtt                                                                                                                        <  2:153214/149‑1 (MQ=255)
                               ccgggccgggTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTAcc                                                                                                                    <  4:206445/149‑1 (MQ=255)
                                   gccgggTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCtt                                                                                                                <  4:18759/149‑1 (MQ=255)
                                    ccgggTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCttt                                                                                                               <  8:179975/149‑1 (MQ=255)
                                     cgggTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCTCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTg                                                                                                              >  4:4896/1‑149 (MQ=255)
                                           acaGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCttt                                                                                                               <  6:191468/142‑1 (MQ=255)
                                           acaGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCttt                                                                                                               >  5:191468/1‑142 (MQ=255)
                                            caGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGTGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATCCCAGGTGGTTTACCTCTTTGCCTGCaa                                                                                                       >  4:189967/1‑149 (MQ=255)
                                                  ggCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGt                                                                                                 <  4:30959/149‑1 (MQ=255)
                                                     tGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCtt                                                                                                                <  7:17458/131‑1 (MQ=255)
                                                     tGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCtt                                                                                                                >  8:17458/1‑131 (MQ=255)
                                                      gCAGGGCCAGGAACAACAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGc                                                                                             >  1:4536/1‑149 (MQ=255)
                                                              aGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCaa                                                                                     <  8:16374/149‑1 (MQ=255)
                                                                    aGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCTCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATtcgtc                                                                               <  5:328829/149‑1 (MQ=255)
                                                                     gCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATtcgtcg                                                                              <  5:153706/149‑1 (MQ=255)
                                                                     gCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATtcgtcg                                                                              >  5:87296/1‑149 (MQ=255)
                                                                     gCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGATGCCTCAATtcgtcg                                                                              >  8:203466/1‑149 (MQ=255)
                                                                      cAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATtcgtcg                                                                              >  2:193984/1‑148 (MQ=255)
                                                                       aaTAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACctct                                                                                                                 <  7:129595/112‑1 (MQ=255)
                                                                       aaTAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACctct                                                                                                                 >  8:129595/1‑112 (MQ=255)
                                                                           aCAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAggtggt                                                                                                                         >  7:263781/1‑100 (MQ=255)
                                                                           aCAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAggtggt                                                                                                                         <  8:263781/100‑1 (MQ=255)
                                                                             aGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGTCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTg                                                                      >  3:33167/1‑149 (MQ=255)
                                                                             aGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGTCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTg                                                                      >  5:79131/1‑149 (MQ=255)
                                                                             aGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCttt                                                                                                               >  6:216077/1‑108 (MQ=255)
                                                                             aGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCttt                                                                                                               <  5:216077/108‑1 (MQ=255)
                                                                             aGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCGGCCTCAATTCGTCGGAGGATTg                                                                      >  6:11911/1‑149 (MQ=255)
                                                                                        cAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCa                                                           <  1:193984/149‑1 (MQ=255)
                                                                                                tGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACa                                                   <  7:254000/149‑1 (MQ=255)
                                                                                                  gggTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTGGGATGATTgg                                                                     >  2:129683/1‑129 (MQ=255)
                                                                                                  gggTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTgg                                                                     <  1:129683/129‑1 (MQ=255)
                                                                                                     tCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCt                                              >  8:52041/1‑149 (MQ=255)
                                                                                                     tCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCt                                              <  2:152527/149‑1 (MQ=255)
                                                                                                       gCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCa                                                                                                   <  8:58990/94‑1 (MQ=255)
                                                                                                       gCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCa                                                                                                   >  7:58990/1‑94 (MQ=255)
                                                                                                        cATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGTCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGAc                                           >  2:18597/1‑149 (MQ=255)
                                                                                                        cATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAg                                                                <  4:29587/128‑1 (MQ=255)
                                                                                                        cATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAg                                                                >  3:29587/1‑128 (MQ=255)
                                                                                                         aTCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGATGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGAc                                           >  6:215080/1‑148 (MQ=255)
                                                                                                           ccTTGAAAGCACGCAACGAGCCTTGCAACAAAGACGCATCCCCACAACTTGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGAAACCTCAATTCCTCGGATGAATGGTACAGTGGCATGAATACAGAGATGAcc                                          >  8:89732/1‑147 (MQ=255)
                                                                                                            cTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCaa                                        <  6:31928/148‑1 (MQ=255)
                                                                                                             ttGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCtt                                                                                                                >  4:82616/1‑75 (MQ=255)
                                                                                                             ttGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCtt                                                                                                                <  3:82616/75‑1 (MQ=255)
                                                                                                                aaaGCACGCAACGAGCGTTGCAACAAAGTGGCATGTCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAgcgc                                    <  4:33167/148‑1 (MQ=255)
                                                                                                                aaaGCACGCAACGAGCGTTGCAACAAAGTGGCATGTCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAgcgc                                    <  6:79131/148‑1 (MQ=255)
                                                                                                                 aaGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTc                                  >  7:244527/1‑149 (MQ=255)
                                                                                                                     ttgCAACGCTCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCtg                              >  3:248231/3‑149 (MQ=255)
                                                                                                                     aCGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCtg                              >  7:110469/1‑149 (MQ=255)
                                                                                                                      cGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCtgg                             >  3:272095/1‑149 (MQ=255)
                                                                                                                                gTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCTGGTGGAGTTGGc                   <  7:52041/149‑1 (MQ=255)
                                                                                                                                 ttGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTc                                  <  8:61927/133‑1 (MQ=255)
                                                                                                                                 ttGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTc                                  >  7:61927/1‑133 (MQ=255)
                                                                                                                                   gCAACAAAGTGGCATGCTCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCTGGTGGAGTTGGCgg                 <  3:4896/148‑1 (MQ=255)
                                                                                                                                    caacaaAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCTGGTGGAGTTGGCGGTa               <  5:215080/149‑1 (MQ=255)
                                                                                                                                     aacaaAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATtcg                                                                                 >  1:321870/1‑82 (MQ=255)
                                                                                                                                     aacaaAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATtcg                                                                                 <  2:321870/82‑1 (MQ=255)
                                                                                                                                      acaaAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACagag                                                >  4:17101/1‑114 (MQ=255)
                                                                                                                                      acaaAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACagag                                                <  3:17101/114‑1 (MQ=255)
                                                                                                                                        aaaGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCTGGTGGAGTTGGCGGTAGAgc           >  1:224462/1‑149 (MQ=255)
                                                                                                                                                aTGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCTGGTGGAGTTGGCGGTAGAGCACATgcgc   <  8:110469/149‑1 (MQ=255)
                                                                                                                                                  gCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAAACCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAACCCCTCCCTGGTGGAGTTGGCGGTAGAGCGCATGcgcg  >  7:238626/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGATGAGCAAGGTTTTAGGGCGGCTCGGGTGCCGGGCCGGGTCACAGCCCGGCTGCAGGGCCAGGAACAGCAATAACAGCAAACTTGCCAGTCGCATGGGGTCGCATCCTTGAAAGCACGCAACGAGCGTTGCAACAAAGTGGCATGCCCACAACTGGCAGAAATACCAGGTGGTTTACCTCTTTGCCTGCAATCCAGTAAGCTGCCTCAATTCGTCGGATGATTGGTTGAGTGGCATGAATACAGAGCTGACCAAGCGCTCCCTGGTGGAGTTGGCGGTAGAGCGCATGCGCG  >  NC_002947/1109099‑1109392

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