Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,984,262 A→T 5.1% L140Q (CTG→CAG)  PP_3512 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,984,2620AT 5.1% 249.0 / 2.5 77L140Q (CTG→CAG) PP_3512membrane protein
Reads supporting (aligned to +/- strand):  ref base A (32/41);  new base T (2/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGCACCACACCACCAG  >  NC_002947/3984114‑3984407
                                                                                                                                                    |                                                                                                                                                 
tACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACcag                                                                                                                                                  >  1:1179640/1‑150 (MQ=255)
          ccTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCgg                                                                                                                                        >  2:970801/1‑150 (MQ=255)
          ccTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCgg                                                                                                                                        <  1:1576637/150‑1 (MQ=255)
           cTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCTGCTGGTGcagcg                                                                                                                                       <  2:201881/150‑5 (MQ=255)
               tAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCttt                                                                                                                                   >  2:1010190/1‑150 (MQ=255)
               tAGGAGCGGGGTTACCCGGGAATAGGCCACACCCGCCAGCCAATAGTGCCTGGGAGATCACCGTTGTCCGTGCAGCCGCATTAGGGGGCAAGCACCCGCCGGGTGACAGAATCATGCCCCACGCCGCTGCACCAGCAGTTGCCGGGCttt                                                                                                                                   <  2:285915/150‑1 (MQ=255)
                  gAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTc                                                                                                                                <  1:1663981/150‑1 (MQ=255)
                  gAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTc                                                                                                                                <  1:970801/150‑1 (MQ=255)
                  gAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTc                                                                                                                                >  2:380497/1‑150 (MQ=255)
                  aaGCGGGTTTATCCGCGAATAGGCCATACCCGCCAGCCAATAGTGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTATCGGGCAAGCCCGCGTCTTGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTc                                                                                                                                <  2:836699/149‑1 (MQ=255)
                    gCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGAATGATGACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGa                                                                                                                              >  1:355499/1‑150 (MQ=255)
                     cGGGTTTACCCGCGAATAGGCTACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGaa                                                                                                                             <  2:355499/150‑1 (MQ=255)
                      gggTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGaa                                                                                                                             <  2:1209757/149‑1 (MQ=255)
                      gggTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGaa                                                                                                                             >  1:1209757/1‑149 (MQ=255)
                       ggTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACg                                                                                                                           <  1:73534/150‑1 (MQ=255)
                        gTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGc                                                                                                                          >  2:214585/1‑150 (MQ=255)
                            aCCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAg                                                                                                                      <  1:1083504/150‑1 (MQ=255)
                               cgcgAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTg                                                                                                                   <  1:149965/150‑1 (MQ=255)
                                gcgAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGt                                                                                                                  >  1:203186/1‑150 (MQ=255)
                                   aaTAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGa                                                                                                               <  2:1753235/150‑1 (MQ=255)
                                       ggCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCAc                                                                                                           >  2:1912463/1‑150 (MQ=255)
                                        gCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCaag                                                                                                          >  2:1458310/1‑148 (MQ=255)
                                        gCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACg                                                                                                          >  2:541276/1‑150 (MQ=255)
                                        gCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACg                                                                                                          >  2:983929/1‑150 (MQ=255)
                                        gCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCAGCATTAGCGAGCAAACCCGCGCCAGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACg                                                                                                          >  1:1555492/1‑150 (MQ=255)
                                         ccACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACgg                                                                                                         <  1:1359581/150‑1 (MQ=255)
                                                cGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAggg                                                                                                  <  1:1620934/150‑1 (MQ=255)
                                                cGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAggg                                                                                                  >  2:1157315/1‑150 (MQ=255)
                                                 gccagccaATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGc                                                                                                 <  1:1421049/150‑1 (MQ=255)
                                                    agccaATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAgg                                                                                              <  2:659902/150‑1 (MQ=255)
                                                     gccaATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGga                                                                                             <  1:1010190/150‑1 (MQ=255)
                                                     gccaATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGga                                                                                             <  1:983929/150‑1 (MQ=255)
                                                      ccaATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACg                                                                                                          >  2:204090/1‑136 (MQ=255)
                                                      ccaATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACg                                                                                                          <  1:204090/136‑1 (MQ=255)
                                                       caATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGcatg                                                                                           >  1:1823145/1‑150 (MQ=255)
                                                           aatGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCTAACGCCCAGTTGTAGACCACCGTGTAGGGCAGGATGATCa                                                                                       <  1:451164/150‑1 (MQ=255)
                                                           aatGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCa                                                                                       <  1:1698019/150‑1 (MQ=255)
                                                               ccTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACgaa                                                                                   <  2:203186/150‑1 (MQ=255)
                                                                cTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCggg                                                                                                                                       >  1:274251/1‑97 (MQ=255)
                                                                cTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCggg                                                                                                                                       <  2:274251/97‑1 (MQ=255)
                                                                cTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAAc                                                                                                                            <  2:78961/108‑1 (MQ=255)
                                                                cTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAAc                                                                                                                            >  1:78961/1‑108 (MQ=255)
                                                                 tGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACgaaga                                                                                 >  2:1844426/1‑150 (MQ=255)
                                                                   ggAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGaagaag                                                                               <  2:1179640/150‑1 (MQ=255)
                                                                       aTCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCa                                                                           <  1:214585/150‑1 (MQ=255)
                                                                        tCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAg                                                                          <  2:1410034/150‑1 (MQ=255)
                                                                        tCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAg                                                                          <  1:541276/150‑1 (MQ=255)
                                                                                    gTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGaa                                                              <  2:772752/150‑1 (MQ=255)
                                                                                       cAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGc                                                           >  1:1867562/1‑150 (MQ=255)
                                                                                       cAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGc                                                           <  2:1555492/150‑1 (MQ=255)
                                                                                       cAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGc                                                           <  2:1292113/150‑1 (MQ=255)
                                                                                        aGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGccc                                                                            <  2:969998/132‑1 (MQ=255)
                                                                                        aGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGccc                                                                            >  1:969998/1‑132 (MQ=255)
                                                                                        aGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGcc                                                          >  2:432681/1‑150 (MQ=255)
                                                                                        aGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGcc                                                          >  1:1012281/1‑150 (MQ=255)
                                                                                         gCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGGGTACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCt                                                         <  1:630123/150‑1 (MQ=255)
                                                                                          ccGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTg                                                        <  1:380497/150‑1 (MQ=255)
                                                                                               ttAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTgcagca                                                   <  1:1912463/150‑1 (MQ=255)
                                                                                                          gCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTcc                                        >  2:136214/1‑150 (MQ=255)
                                                                                                          gCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTcc                                        >  2:799441/1‑150 (MQ=255)
                                                                                                          gCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGACGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTgcag                                                     >  2:594667/1‑137 (MQ=255)
                                                                                                          gCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGACGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTgcag                                                     <  1:594667/137‑1 (MQ=255)
                                                                                                                  cGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTATGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCa                                >  2:789524/1‑150 (MQ=255)
                                                                                                                   gggTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCaa                               <  1:76346/150‑1 (MQ=255)
                                                                                                                        aCAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGcc                                                          >  1:318474/1‑118 (MQ=255)
                                                                                                                        aCAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGcc                                                          <  2:318474/118‑1 (MQ=255)
                                                                                                                               cATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCa                                                                                       <  1:1925494/82‑1 (MQ=255)
                                                                                                                               cATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCa                                                                                       >  2:1925494/1‑82 (MQ=255)
                                                                                                                                   ccACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAAt                                                             >  1:1756889/1‑104 (MQ=255)
                                                                                                                                   ccACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAAt                                                             <  2:1756889/104‑1 (MQ=255)
                                                                                                                                    cacaCGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCggcaggca              >  2:886762/1‑150 (MQ=255)
                                                                                                                                    cacaCGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCggcaggca              >  2:1694582/1‑150 (MQ=255)
                                                                                                                                         gccgcTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGcaccaca         >  2:1561753/1‑150 (MQ=255)
                                                                                                                                          gcccggCACCTGCTGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGcaccacac        >  1:1344643/6‑150 (MQ=255)
                                                                                                                                          ccgcTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGcaccacac        <  1:432681/150‑1 (MQ=255)
                                                                                                                                            ccggCACCTGCTGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAAt                                                             <  1:883828/92‑1 (MQ=255)
                                                                                                                                            ccggCACCTGCTGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAAt                                                             >  2:883828/4‑95 (MQ=255)
                                                                                                                                                cACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGCACCACACCACCAg  <  1:1844426/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGCACCACACCACCAG  >  NC_002947/3984114‑3984407

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑