Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 468,464 A→G 6.6% C831C (TGT→TGC PP_0386 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947468,4640AG 6.6% 265.1 / 6.3 94C831C (TGT→TGCPP_0386sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (37/49);  new base G (3/3);  total (42/52)
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.59e-01

TCAACGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGTGAT  >  NC_002947/468318‑468610
                                                                                                                                                   |                                                                                                                                                   
tCAACGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGca                                                                                                                                                     >  5:274690/1‑147 (MQ=255)
  aaCGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAgg                                                                                                                                                 >  6:306246/1‑149 (MQ=255)
   aCGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGt                                                                                                                                                <  4:141509/149‑1 (MQ=255)
   aCGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGt                                                                                                                                                <  8:262293/149‑1 (MQ=255)
     gTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGGCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCTCTGTGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                              <  6:124164/149‑1 (MQ=255)
     gTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                              <  5:123467/149‑1 (MQ=255)
     gTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                              <  6:309139/149‑1 (MQ=255)
       cgcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGt                                                                                                                                            <  7:27910/149‑1 (MQ=255)
        gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAAGGACCCCGCCCGGGTGGCCAAGGCTGACGTTGATGTGCTTGAGGCGGTCGAGGTCCAGGAACAGCACTCCGCCCTGGCGGTTGGACACCTGCGCACAGGTCAg                                                                                                                                             >  1:73299/1‑148 (MQ=255)
        gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGTGCGCAGGTgtcc                                                                                                                                            <  4:322505/148‑4 (MQ=255)
        gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGt                                                                                                                                            <  1:339722/148‑1 (MQ=255)
        gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGt                                                                                                                                            <  5:76857/148‑1 (MQ=255)
        gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTa                                                                                                                                           >  1:71054/1‑149 (MQ=255)
         cgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                              <  6:61840/145‑1 (MQ=255)
         cgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                              >  5:61840/1‑145 (MQ=255)
             ccTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                              >  5:29440/1‑141 (MQ=255)
             ccTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                              <  6:29440/141‑1 (MQ=255)
             ccTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGccc                                                                                                                                      <  7:278519/149‑1 (MQ=255)
              cTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCt                                                                                                                                     >  5:183399/1‑149 (MQ=255)
                gCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTg                                                                                                                                    >  5:85009/1‑148 (MQ=255)
                actctTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCCCAGGTGAGTACGCCCTg                                                                                                                                    >  5:224544/2‑148 (MQ=255)
                actctTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTg                                                                                                                                    <  6:224544/147‑1 (MQ=255)
                  tctTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAg                                                                                                                                 >  8:35962/1‑149 (MQ=255)
                    tttCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGcc                                                                                                                               >  6:230838/1‑149 (MQ=255)
                      tCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGc                                                                                                                                <  2:97017/146‑1 (MQ=255)
                           cgcTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                        <  6:183399/149‑1 (MQ=255)
                           cgcTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                        <  5:230838/149‑1 (MQ=255)
                                   gATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAgg                                                                                                                <  5:306246/149‑1 (MQ=255)
                                        ccTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACaa                                                                                                                  >  2:35470/1‑142 (MQ=255)
                                        ccTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACaa                                                                                                                  <  1:35470/142‑1 (MQ=255)
                                        ccTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGa                                                                                                           >  6:109425/1‑149 (MQ=255)
                                        ccTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGAGGTTGATGTCCTTGAAGCGGTCGAGGTCCAGGAACAGCACTTCGCCCTGGCGGTTGGACAACTGCGAAAACGTGACTACGAGCTGGAGCCGGTGTTCGACCAATGTGCGa                                                                                                           >  8:155306/1‑149 (MQ=255)
                                         cTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGAt                                                                                                          >  8:339904/1‑149 (MQ=255)
                                          ttGAGCAGCAGGTCGCAGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGAtt                                                                                                         >  2:259404/1‑149 (MQ=255)
                                             agcagcAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTgg                                                                                                       <  4:201441/148‑1 (MQ=255)
                                               cagcagGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTggg                                                                                                      >  1:217041/1‑147 (MQ=255)
                                               cagcagGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTggg                                                                                                      >  5:129245/1‑147 (MQ=255)
                                               cagcagGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCa                                                                                                    <  5:109425/149‑1 (MQ=255)
                                                  cagGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGc                                                                                                 <  1:362673/149‑1 (MQ=255)
                                                           cGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCTCTGGTGTGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGc                                                                                        >  1:312677/1‑149 (MQ=255)
                                                            gACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCg                                                                                       >  5:90979/1‑149 (MQ=255)
                                                                ggggggCCACGGCTGTCGTTGCTGTGCTTGTCGCGGTCGAGGTCCCGGAACAGCTCTGCGCCCTGGCGGTTGGACTCCTTCGCACAGGGGAGTACGCCCGGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTc                                                                                   <  8:265794/149‑1 (MQ=255)
                                                                  gTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCgt                                                                                 >  4:245877/1‑149 (MQ=255)
                                                                   tGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCgtg                                                                                <  3:274041/149‑1 (MQ=255)
                                                                      ccAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGcc                                                                                                                               <  6:265047/100‑1 (MQ=255)
                                                                      ccAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGcc                                                                                                                               >  5:265047/1‑100 (MQ=255)
                                                                       cAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGGCGtggtgg                                                                            <  6:274690/149‑1 (MQ=255)
                                                                         aGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGc                                                                          <  3:348806/149‑1 (MQ=255)
                                                                         aGGCTGTCGTTGATGTGCTGGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGc                                                                          <  2:73299/149‑1 (MQ=255)
                                                                           gCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCtt                                                                        <  2:217041/149‑1 (MQ=255)
                                                                             tGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGa                                                                      <  3:116247/149‑1 (MQ=255)
                                                                              gTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTCCGCCCTGGAGCCGGGTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGGTTGAt                                                                     >  8:266763/1‑149 (MQ=255)
                                                                               tCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGa                                                                      >  3:172081/1‑147 (MQ=255)
                                                                               tCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATa                                                                    <  7:35962/149‑1 (MQ=255)
                                                                                 gTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                        <  8:221462/96‑1 (MQ=255)
                                                                                 gTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                        >  5:9592/1‑96 (MQ=255)
                                                                                 gTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                        >  7:221462/1‑96 (MQ=255)
                                                                                 gTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                        <  6:9592/96‑1 (MQ=255)
                                                                                  ttGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTc                                                                 >  3:261516/1‑149 (MQ=255)
                                                                                     aTGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAAGCGGGTCGTGGTGGGCTTGATAGTCGa                                                               <  1:11802/149‑1 (MQ=255)
                                                                                     aTGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAg                                                              <  5:150972/149‑1 (MQ=255)
                                                                                       gtgCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGc                                                                                        >  3:26844/1‑121 (MQ=255)
                                                                                       gtgCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGc                                                                                        <  4:26844/121‑1 (MQ=255)
                                                                                          cTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTg                                                          <  1:259404/148‑1 (MQ=255)
                                                                                          cTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTg                                                          <  6:85009/148‑1 (MQ=255)
                                                                                            tGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCg                                                                                  <  8:132998/122‑1 (MQ=255)
                                                                                            tGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCg                                                                                  >  7:132998/1‑122 (MQ=255)
                                                                                            tGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGcc                                                       >  4:118716/1‑149 (MQ=255)
                                                                                            tGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGcc                                                       >  8:210356/1‑149 (MQ=255)
                                                                                              aaGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTg                                                     <  6:90979/149‑1 (MQ=255)
                                                                                                 cGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGc                                                  <  2:71054/149‑1 (MQ=255)
                                                                                                   gTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGt                                                >  7:92284/1‑149 (MQ=255)
                                                                                                    tCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGGACAGGTGAGTACGCCCTGGAGCCGGTTTTCGATCAAGCTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTAGAGTTTGGCCTGGGCGtc                                               >  6:107489/1‑148 (MQ=255)
                                                                                                      gAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCt                                             <  7:339904/149‑1 (MQ=255)
                                                                                                        ggTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGGGTGCtg                                            >  2:342049/1‑147 (MQ=255)
                                                                                                          tCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGa                                                                      <  5:292058/120‑1 (MQ=255)
                                                                                                          tCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGa                                                                      >  6:292058/1‑120 (MQ=255)
                                                                                                             aGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCt                                      <  7:210356/149‑1 (MQ=255)
                                                                                                                   aGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGAt                                <  2:312677/149‑1 (MQ=255)
                                                                                                                     cACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGTGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGt                              >  8:49854/1‑149 (MQ=255)
                                                                                                                       cTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGTGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCt                            <  7:49854/149‑1 (MQ=255)
                                                                                                                              cTGGTGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCaaaaa                     >  2:156893/1‑149 (MQ=255)
                                                                                                                                ggCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAAcc                   >  5:254084/1‑149 (MQ=255)
                                                                                                                                     ttGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCAc              <  4:172081/149‑1 (MQ=255)
                                                                                                                                       actcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTTCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACg             >  5:260039/4‑148 (MQ=255)
                                                                                                                                        gACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGaa           <  4:261516/149‑1 (MQ=255)
                                                                                                                                         tcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGc                                                                                        <  6:306830/70‑1 (MQ=255)
                                                                                                                                         tcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGc                                                                                        >  5:306830/2‑71 (MQ=255)
                                                                                                                                         tcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCaaaa                      <  5:68116/136‑1 (MQ=255)
                                                                                                                                         tcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCaaaa                      >  6:68116/2‑137 (MQ=255)
                                                                                                                                         acacCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGgag                                  <  4:26051/125‑1 (MQ=255)
                                                                                                                                             cTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGga                                   >  8:33602/1‑120 (MQ=255)
                                                                                                                                             cTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGga                                   <  7:33602/120‑1 (MQ=255)
                                                                                                                                              tGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGGGGGt     >  4:141522/1‑149 (MQ=255)
                                                                                                                                                 gcACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGTGAt  <  8:92284/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCAACGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGTGAT  >  NC_002947/468318‑468610

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: