Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,058,488 C→T 7.3% L62L (TTG→TTA) ‡ PP_3576 ← transmembrane sensor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,058,4880CT 7.3% 199.4 / 4.3 68L62L (TTG→TTA) ‡PP_3576transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base C (32/31);  new base T (2/3);  total (34/34)
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.61e-01

GGCGGTCCAGGCATCCAGCGGCAGCTCACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATGCCAGCCAACGGGCTGCA  >  NC_002947/4058347‑4058632
                                                                                                                                             |                                                                                                                                                
ggcggTCCAGGCATCCAGCGGCAGCTCACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCCGGGCGTGTTTGAGCAAGGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCgg                                                                                                                                           <  5:233105/149‑1 (MQ=255)
ggcggTCCAGGCATCCAGCGGCAGCTCACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCg                                                                                                                                            <  2:239021/148‑1 (MQ=255)
     tCCAGGCATCCAGCGGCAGCTCACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAAt                                                                                                                                      >  2:16282/1‑149 (MQ=255)
       cAGGCATCCAGCGGCAGCTCACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGTTGCCCGGcgcatt                                                                                                                                     <  1:911/148‑6 (MQ=255)
        aGGCATCCAGCGGCAGCTCACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATcgc                                                                                                                                   <  1:310132/149‑1 (MQ=255)
            aTCCAGCGGCAGCTCACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTc                                                                                                                               >  2:235804/1‑149 (MQ=255)
                          cACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGTTGCCCGGc                                                                                                                                          >  8:286475/1‑124 (MQ=255)
                          cACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGTTGCCCGGc                                                                                                                                          <  7:286475/124‑1 (MQ=255)
                           aCGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGTTGCCCGGcgcatt                                                                                                                                     <  6:225943/128‑6 (MQ=255)
                           aCGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGTTGCCCGGcgcatt                                                                                                                                     >  5:225943/1‑123 (MQ=255)
                             ggCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCttt                                                                                                              >  2:11005/1‑149 (MQ=255)
                             ggCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCttt                                                                                                              >  6:328608/1‑149 (MQ=255)
                                           gCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTc                                                                                                <  1:235804/149‑1 (MQ=255)
                                           gCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTc                                                                                                >  6:175217/1‑149 (MQ=255)
                                               ccGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGc                                                                                            <  1:16282/149‑1 (MQ=255)
                                                     acaacGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGa                                                                                      <  6:282748/149‑1 (MQ=255)
                                                       aacGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGa                                                                                      >  3:345651/1‑147 (MQ=255)
                                                       aacGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGa                                                                                      <  4:345651/147‑1 (MQ=255)
                                                       aacGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGAct                                                                                    >  6:164492/1‑149 (MQ=255)
                                                                acacCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGgcgcg                                                                           >  6:278474/1‑149 (MQ=255)
                                                                 cacCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGgcgcgc                                                                          <  5:175217/149‑1 (MQ=255)
                                                                  acCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGc                                                                                                           >  6:64839/1‑115 (MQ=255)
                                                                  acCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGc                                                                                                           <  5:64839/115‑1 (MQ=255)
                                                                  acCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGTTGgcgcgc                                                                          >  7:109233/1‑148 (MQ=255)
                                                                    ccAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCAt                                                                       >  3:330086/1‑149 (MQ=255)
                                                                      aGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGgcgc                                                                            <  7:254165/142‑1 (MQ=255)
                                                                      aGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGgcgc                                                                            >  8:254165/1‑142 (MQ=255)
                                                                         gCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGAAGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTgcgc                                                                  >  8:54908/1‑149 (MQ=255)
                                                                             gTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGACAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGa                                                              >  3:53391/1‑149 (MQ=255)
                                                                              tttGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGa                                                                                      <  2:228918/124‑1 (MQ=255)
                                                                                  aGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTgcgc                                                                  >  3:158124/1‑140 (MQ=255)
                                                                                  aGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTgcgc                                                                  >  5:142103/1‑140 (MQ=255)
                                                                                  aGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTgcgc                                                                  <  4:158124/140‑1 (MQ=255)
                                                                                  aGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTgcgc                                                                  <  6:142103/140‑1 (MQ=255)
                                                                                   gCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGTTGGCGCGCCATTGCGCCAGACGCTGc                                                        <  8:109233/149‑1 (MQ=255)
                                                                                        gCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTgc                                                                    <  6:237433/132‑1 (MQ=255)
                                                                                           cgacgaCCGGCATCCGCACGAACCAGGGCGGCCATCGCAAATGCGCCGGGCAAACCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTCCTGGCGCGCCATTGCGCCAGACGCTCCAGGACac                                                 >  2:321908/1‑147 (MQ=255)
                                                                                            gtcgacCGGCATCCGGACGATCCAGGGCGGCCATCTCCAATGCGCGTGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTgc                                                                    <  4:2830/126‑1 (MQ=255)
                                                                                            gacgacCGGCATCCGGACGATCCAGGGCGGCCATCCCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCACGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTgc                                                                    >  3:2830/1‑128 (MQ=255)
                                                                                             acgacCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGc                                                                                                      <  1:301172/93‑1 (MQ=255)
                                                                                             acgacCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGc                                                                                                      >  2:301172/1‑93 (MQ=255)
                                                                                                  cGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCtt                                         >  3:142706/1‑149 (MQ=255)
                                                                                                  cGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCtt                                         >  6:70305/1‑149 (MQ=255)
                                                                                                    tcATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGCGCAACCCGCCAAATCGCGCACGCACCACTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTCCCCCAGACGCTGCAGGTCAGATTCGCTAGc                                       >  8:332449/2‑149 (MQ=255)
                                                                                                       tCCGGACGATCCAGGGCGGCCATCGCCAATGAGCCGGGCAACCCGGCAAATCGCCCTCGCAGCAGTGATGCCTTATCCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTc                                    >  7:315345/1‑149 (MQ=255)
                                                                                                        ccGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCg                                   >  6:322938/1‑149 (MQ=255)
                                                                                                         cGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGc                                  <  5:278474/149‑1 (MQ=255)
                                                                                                         cGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGc                                  <  1:321908/149‑1 (MQ=255)
                                                                                                           gACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGc                                  >  3:99780/1‑147 (MQ=255)
                                                                                                            aCGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGc                                  >  8:238793/1‑146 (MQ=255)
                                                                                                                tCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGAtt                                              >  3:63202/1‑130 (MQ=255)
                                                                                                                tCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGAtt                                              <  4:63202/130‑1 (MQ=255)
                                                                                                                  cAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCa                         <  4:330086/149‑1 (MQ=255)
                                                                                                                  cAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCa                         <  2:40598/149‑1 (MQ=255)
                                                                                                                  cAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCa                         <  4:99780/149‑1 (MQ=255)
                                                                                                                   aGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAg                        <  4:142706/149‑1 (MQ=255)
                                                                                                                   aGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAg                        >  4:258183/1‑149 (MQ=255)
                                                                                                                       cggcCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGAGGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAAt                    <  7:332449/149‑1 (MQ=255)
                                                                                                                        ggcCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGGAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATg                   >  1:245135/1‑149 (MQ=255)
                                                                                                                           cATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATgcca                <  2:153375/149‑1 (MQ=255)
                                                                                                                           cATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATgc                  >  1:153375/1‑147 (MQ=255)
                                                                                                                            aTCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATgccag               >  7:108133/1‑149 (MQ=255)
                                                                                                                                 cAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATGCCAGCCAAc          <  7:54908/149‑1 (MQ=255)
                                                                                                                                  aaTGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATGCCAGCCAACg         >  7:71603/1‑149 (MQ=255)
                                                                                                                                   aTGCGCCGGACAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTTCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATGCCAGCCATCgg        <  4:53391/149‑1 (MQ=255)
                                                                                                                                      cgcCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATGCCAGCCAACGGGc      <  5:70305/148‑1 (MQ=255)
                                                                                                                                         cGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATGCCAGCCAACGGGCTgca  <  7:238793/149‑1 (MQ=255)
                                                                                                                                         cGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATGCCAGCCAACGGGCTgca  <  8:108133/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
GGCGGTCCAGGCATCCAGCGGCAGCTCACGGCTCACCAGCCAGGCCACCGGCAACAACGCCGTCACACCCAGGGCTTGTTTGAGCAATGCCCGACGACCGGCATCCGGACGATCCAGGGCGGCCATCGCCAATGCGCCGGGCAACCCGGCAAATCGCGCTCGCAGCAGTGATGCCTTTTGCCAGGCGTCTTCGTGCAGGTGACTGCTGGCGCGCCATTGCGCCAGACGCTGCAGGTCAGATTCGCTTGCGTCGCCAGAATCCAGTAATGCCAGCCAACGGGCTGCA  >  NC_002947/4058347‑4058632

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