Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,321,510 C→A 5.7% K407N (AAG→AAT)  PP_1154 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,321,5100CA 5.7% 225.7 / 3.0 70K407N (AAG→AAT) PP_1154sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (35/31);  new base A (2/2);  total (37/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.11e-01

CGCGCAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCG  >  NC_002947/1321362‑1321657
                                                                                                                                                     |                                                                                                                                                   
cgcgCAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGcc                                                                                                                                                     <  6:324001/149‑1 (MQ=255)
   gcAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTg                                                                                                                                                  >  3:137558/1‑149 (MQ=255)
     aacaaGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGc                                                                                                                                                 >  7:5436/1‑148 (MQ=255)
      acaaGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGcc                                                                                                                                                >  2:192238/1‑148 (MQ=255)
           gCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCAAGGCCATcggt                                                                                                                                          <  8:136115/149‑3 (MQ=255)
             cTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGCATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCgcgc                                                                                                                                         >  7:105047/1‑148 (MQ=255)
             cTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGCTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGAGTCGACTTCGTCGGTGATATCGCTGGCGGTACACCGG‑AAACCGATGGCCTTGCCATCgcgc                                                                                                                                         >  2:96616/1‑148 (MQ=255)
                   cAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGT‑AAACCGATGGCCTTGCCATCGCGCCAGATa                                                                                                                                   <  8:5436/148‑1 (MQ=255)
                      gTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAgcgc                                                                                                                               <  6:337325/149‑1 (MQ=255)
                      gTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAgcgc                                                                                                                               <  1:192238/149‑1 (MQ=255)
                         ccGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCgcgg                                                                                                                             <  6:48409/148‑1 (MQ=255)
                                 gCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGATGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGTTACCCCGG‑AAACCGATGGCCATGCCATCGCGCCAGATAGCGCGGGCGGATATc                                                                                                                    >  7:296347/1‑149 (MQ=255)
                                  cTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGGCAAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATc                                                                                                                    <  1:78976/149‑1 (MQ=255)
                                   ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCg                                                                                                                  >  7:216420/1‑149 (MQ=255)
                                   ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCg                                                                                                                  >  5:191810/1‑149 (MQ=255)
                                   ttgttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTACTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGCACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCg                                                                                                                  >  1:274057/1‑149 (MQ=255)
                                     gttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGCCTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AACCCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                 >  3:9594/1‑148 (MQ=255)
                                     gttgCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                 >  2:63225/1‑148 (MQ=255)
                                     gttgCGGTTGGCCAGGCCGGGCAGTGGGTCGTGCAGCGAAAGGGGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGCTGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                 <  8:266078/148‑1 (MQ=255)
                                        gCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAg                                                                                                              >  7:351977/1‑148 (MQ=255)
                                              ggCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCgg                                                                                                        <  7:332687/148‑1 (MQ=255)
                                                  aGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                 >  7:168370/1‑135 (MQ=255)
                                                  aGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCgg                                                                                                                 <  8:168370/135‑1 (MQ=255)
                                                   gTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTgg                                                                                                  <  4:137558/149‑1 (MQ=255)
                                                      cGGGCAGTGGGTCGTGCAGAGAAATGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGcc                                                                                                <  1:179987/148‑1 (MQ=255)
                                                         tCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTc                                                                                            <  6:191810/149‑1 (MQ=255)
                                                            gTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGTGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTgg                                                                                         >  2:99371/1‑149 (MQ=255)
                                                            gTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCCGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCAGTCTgg                                                                                         >  4:240613/1‑149 (MQ=255)
                                                                 tCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTcgcg                                                                                    >  2:79567/1‑149 (MQ=255)
                                                                   gTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCgg                                                                                   >  4:98411/1‑148 (MQ=255)
                                                                     gCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTaa                                                                                <  6:291416/149‑1 (MQ=255)
                                                                      cAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTaa                                                                                <  2:206306/148‑1 (MQ=255)
                                                                      cAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAg                                                                               >  7:76094/1‑149 (MQ=255)
                                                                      cAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGCGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTaa                                                                                <  1:372780/148‑1 (MQ=255)
                                                                       aGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGt                                                                              <  4:365595/149‑1 (MQ=255)
                                                                         cGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTAc                                                                            <  1:63225/149‑1 (MQ=255)
                                                                           aaaGGTGTTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAg                                                                          >  3:189091/1‑149 (MQ=255)
                                                                             aGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAgcg                                                                        >  1:159935/1‑149 (MQ=255)
                                                                                  cTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGt                                                                   >  5:110189/1‑149 (MQ=255)
                                                                                    gAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTgg                                                                                                  >  2:330592/1‑116 (MQ=255)
                                                                                    gAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTgg                                                                                                  <  1:330592/116‑1 (MQ=255)
                                                                                     aaTGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTgg                                                                >  1:270668/1‑149 (MQ=255)
                                                                                      aTGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGc                                                               <  8:351977/149‑1 (MQ=255)
                                                                                            gcgTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATtg                                                         >  7:181749/1‑149 (MQ=255)
                                                                                                 ggCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATtgctgc                                                     <  4:189091/148‑1 (MQ=255)
                                                                                                        aCTTCGTCGGTGATATCGCTGGCGGCACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctcggc                                             <  2:274057/149‑1 (MQ=255)
                                                                                                         cTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATtgctgc                                                     >  5:211523/1‑140 (MQ=255)
                                                                                                         cTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATtgctgc                                                     <  6:211523/140‑1 (MQ=255)
                                                                                                         cTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGcggctcggct                                            <  2:159935/149‑1 (MQ=255)
                                                                                                            cgtcgGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCt                                                                                                    >  5:126904/1‑90 (MQ=255)
                                                                                                            cgtcgGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCt                                                                                                    <  6:126904/90‑1 (MQ=255)
                                                                                                             gtcgGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCa                                         <  8:76094/148‑1 (MQ=255)
                                                                                                                ggTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGt                                       >  3:318507/1‑147 (MQ=255)
                                                                                                                 gTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAg                                    >  7:68201/1‑149 (MQ=255)
                                                                                                                     tatCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCg                                >  5:259267/1‑149 (MQ=255)
                                                                                                                       tCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTc                              >  5:61077/1‑149 (MQ=255)
                                                                                                                       tCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTc                              <  8:68201/149‑1 (MQ=255)
                                                                                                                          cTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAg                           >  6:42867/1‑149 (MQ=255)
                                                                                                                             gCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGcc                         <  4:9594/148‑1 (MQ=255)
                                                                                                                               ggTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGcc                         <  4:197958/146‑1 (MQ=255)
                                                                                                                               ggTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGcc                         >  3:197958/1‑146 (MQ=255)
                                                                                                                               ggTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAg                      >  8:16492/1‑149 (MQ=255)
                                                                                                                                           gcgcGATGGCAAGGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgcggc           >  6:77955/4‑149 (MQ=255)
                                                                                                                                           aaCCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgcggc           <  1:96616/149‑1 (MQ=255)
                                                                                                                                             gcGATGGCAAGGCCATCGCGCCAGATAGCGGAGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCa                            <  8:205591/129‑1 (MQ=255)
                                                                                                                                             gcGATGGCAAGGCCATCGCGCCAGATAGCGCAGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCa                            >  7:205591/2‑130 (MQ=255)
                                                                                                                                             ccGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGt         <  6:110189/149‑1 (MQ=255)
                                                                                                                                              cGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTg        <  3:98411/149‑1 (MQ=255)
                                                                                                                                              cGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTACGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTg        >  6:24689/1‑149 (MQ=255)
                                                                                                                                              cGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGGTGGCCAATTGCTGCCGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTg        <  1:79567/149‑1 (MQ=255)
                                                                                                                                                aTGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGt      >  2:86847/1‑149 (MQ=255)
                                                                                                                                                    ccTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgcggc           >  3:255039/1‑140 (MQ=255)
                                                                                                                                                    ccTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAgcggc           <  4:255039/140‑1 (MQ=255)
                                                                                                                                                     cTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCg  >  8:246615/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGCGCAACAAGGCCTGCTCCAGGTGCCGGGCGAGCTTGTTGCGGTTGGCCAGTCCGGTCAGTGGGTCGTGCAGCGAAAGGTGCTGAATGCGCGCGTGGGCGTCGACTTCGTCGGTGATATCGCTGGCGGTACCCCGG‑AAACCGATGGCCTTGCCATCGCGCCAGATAGCGCGGGCGGATATCCGGCAGTAGCGGTTCTGGCCATTCTGGTCGCGGTAAGTACAGCGCAGGTTGGCCAATTGCTGCGGCTCGGCTGCAGTCAGGGCGTCCAGCCAAGGTGACAGCGGCGTGGTATCG  >  NC_002947/1321362‑1321657

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