Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 3,316,755 G→A R258C (CGT→TGT)  nusA ← transcription termination/antitermination L factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi3,316,7550GA91.5% 340.3 / 16.6 117R258C (CGT→TGT) nusAtranscription termination/antitermination L factor
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base A (61/46);  total (64/53)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

Only 100 of 121 total aligned reads displayed.

CATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGA  >  NC_000913_3_pae_pgi/3316612‑3316901
                                                                                                                                               |                                                                                                                                                  
cATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGcc                                                                                                                                              <  1:697172/150‑1 (MQ=255)
       tgtTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATa                                                                                                                                       >  1:1106399/1‑150 (MQ=255)
           tATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCTGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTa                                                                                                                                   <  2:584405/150‑1 (MQ=255)
              cTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGc                                                                                                                                <  1:344641/150‑1 (MQ=255)
              cTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGc                                                                                                                                >  1:278148/1‑150 (MQ=255)
              cTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGc                                                                                                                                >  3:649395/1‑150 (MQ=255)
                tCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCaa                                                                                                                              >  4:172398/1‑150 (MQ=255)
                  aTCCACCACGATAGAAGCAACGTCTGTCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGc                                                                                                                            >  1:1373461/1‑150 (MQ=255)
                  aTCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGAATACCTACGCAAGc                                                                                                                            >  2:67722/1‑150 (MQ=255)
                    ccaccaCGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCAc                                                                                                                          <  2:899998/150‑1 (MQ=255)
                    ccaccaCGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCAc                                                                                                                          <  3:390388/150‑1 (MQ=255)
                       ccacGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTa                                                                                                                       >  1:916817/1‑150 (MQ=255)
                        cacGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTAc                                                                                                                      <  3:373686/150‑1 (MQ=255)
                        cacGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTAc                                                                                                                      >  3:129034/1‑150 (MQ=255)
                          cGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCg                                                                                                                    <  2:1426753/150‑1 (MQ=255)
                             tAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGa                                                                                                                  >  1:77711/1‑149 (MQ=255)
                             tAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGAt                                                                                                                 <  4:805541/150‑1 (MQ=255)
                               gAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCAGGCATACCTACGCAAGCACCTACCGGATCg                                                                                                               >  2:1066120/1‑150 (MQ=255)
                                aaGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGa                                                                                                              >  4:747579/1‑150 (MQ=255)
                                aaGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGa                                                                                                              >  2:1078426/1‑150 (MQ=255)
                                aaGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCCACGCAAGCACCCACCGGCTCGa                                                                                                              >  4:177191/1‑150 (MQ=255)
                                 aGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGAt                                                                                                             <  1:1348772/150‑1 (MQ=255)
                                   cAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCAt                                                                                                                                        >  4:967683/1‑121 (MQ=255)
                                            ccGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGt                                                                                                  <  3:14102/150‑1 (MQ=255)
                                             cGGTGCCAGTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGtt                                                                                                 <  2:1261478/150‑1 (MQ=255)
                                              ggTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTg                                                                                                >  3:7594/1‑150 (MQ=255)
                                               gTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCg                                                                                                                    <  2:1695559/129‑1 (MQ=255)
                                               gTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCg                                                                                                                    >  1:1695559/1‑129 (MQ=255)
                                                   cATTGCGTTAATCACGGACTGCGCCGGGTTATCATCCCCCGCCACCATATCGACGCCTTCGTCCCCCATTTCAGTAGACACCGCCTGAACACACGCGCCACTCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGtttt                                                                                           <  1:67722/150‑1 (MQ=255)
                                                   cATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGtttt                                                                                           >  2:1447630/1‑150 (MQ=255)
                                                     ttGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACAcacgcggtgt                                                                                                                                            >  1:1431490/1‑95 (MQ=255)
                                                     ttGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCa                                                                                         >  2:537096/1‑150 (MQ=255)
                                                        cGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGAGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcg                                                                                      >  1:1441430/1‑150 (MQ=255)
                                                        cGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcg                                                                                      >  2:45624/1‑150 (MQ=255)
                                                         gTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACcgc                                                                                     >  4:746041/1‑150 (MQ=255)
                                                           tAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGa                                                                                   <  2:916817/150‑1 (MQ=255)
                                                             aTCACGAACTGCGCCCGGTTATCATCCCACAGAACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCGGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGAtt                                                                                 <  4:1024021/150‑1 (MQ=255)
                                                               cACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGCCACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTACCGTTGGCTTTCCCCGCGAtttc                                                                               >  2:975325/1‑149 (MQ=255)
                                                                 cGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCg                                                                             <  4:7594/150‑1 (MQ=255)
                                                                   aaCTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCa                                                                           <  4:1012767/150‑1 (MQ=255)
                                                                      tGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACga                                                                        >  2:762991/1‑150 (MQ=255)
                                                                       gcgcCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACgag                                                                       >  1:347160/1‑150 (MQ=255)
                                                                         gcCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTg                                                                                                <  4:894702/123‑1 (MQ=255)
                                                                         gcCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTg                                                                                                >  3:894702/1‑123 (MQ=255)
                                                                              gTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCgg                                                                <  1:1566694/150‑1 (MQ=255)
                                                                               ttATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAAGCCGGa                                                               >  2:1588736/1‑150 (MQ=255)
                                                                                tATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGAt                                                              <  1:416264/150‑1 (MQ=255)
                                                                                 atcatcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGAt                                                                                                             <  4:759653/102‑1 (MQ=255)
                                                                                 atcatcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGAt                                                                                                             >  3:759653/1‑102 (MQ=255)
                                                                                  tcatcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcg                                                            <  1:1022924/150‑1 (MQ=255)
                                                                                    atcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgcg                                                          >  1:1568173/1‑150 (MQ=255)
                                                                                    atcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgcg                                                          >  1:1568170/1‑150 (MQ=255)
                                                                                    atcCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATcgcg                                                          >  2:1551569/1‑150 (MQ=255)
                                                                                       ccACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGc                                                       >  3:790127/1‑150 (MQ=255)
                                                                                        cacaGGACGATATCGATACGCGCGCCACCCATGTCAGTAGACACCGCCTGAACACACGCGCCACGCATCCCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGcc                                                      >  2:667801/1‑150 (MQ=255)
                                                                                         acaGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCg                                                     >  4:749631/1‑150 (MQ=255)
                                                                                             gACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCgctgc                                                 >  2:518749/1‑150 (MQ=255)
                                                                                                  tatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTCCGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTaa                                            <  3:739489/150‑1 (MQ=255)
                                                                                                  tatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTaa                                            <  3:689601/150‑1 (MQ=255)
                                                                                                  tatCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTaa                                            <  1:518749/150‑1 (MQ=255)
                                                                                                   atCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTATCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAAt                                           >  1:708920/1‑150 (MQ=255)
                                                                                                    tCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAAtt                                          >  4:629964/1‑150 (MQ=255)
                                                                                                      gATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTc                                        <  1:764399/150‑1 (MQ=255)
                                                                                                      gATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTc                                        <  1:1551569/150‑1 (MQ=255)
                                                                                                      gATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGAACGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCCGCTTTAAttgt                                        >  4:427051/1‑148 (MQ=255)
                                                                                                        tACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCtttcatttcat                                      <  2:1366406/150‑2 (MQ=255)
                                                                                                         aCGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAAt                                     <  1:266000/150‑1 (MQ=255)
                                                                                                           gCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCGCGCGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCa                                   >  1:865109/1‑150 (MQ=255)
                                                                                                           gCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCa                                   >  4:597127/1‑150 (MQ=255)
                                                                                                             tCGCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCAct                                 >  3:161204/1‑150 (MQ=255)
                                                                                                               gCCACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActtc                               >  4:55682/1‑150 (MQ=255)
                                                                                                                 cACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCg                                                                 <  2:129192/114‑1 (MQ=255)
                                                                                                                 cACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCg                                                                 >  1:129192/1‑114 (MQ=255)
                                                                                                                 cACCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCActtctt                             >  3:605142/1‑150 (MQ=255)
                                                                                                                  aCCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACttcttc                            >  1:828484/1‑150 (MQ=255)
                                                                                                                  aCCCAGTTCAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTAATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACttcttc                            >  1:828461/1‑150 (MQ=255)
                                                                                                                    ccAGTTCAGTAGACACCGCCTGAACAAAAAAGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGc                          >  2:849352/1‑150 (MQ=255)
                                                                                                                       gTTCAGTAGACACCGCCTGAACACACGCGCCCCGCCTACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCgc                       <  4:100388/150‑2 (MQ=255)
                                                                                                                          cAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGAt                                                                                                             >  3:981187/1‑61 (MQ=255)
                                                                                                                          cAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAttt                    >  1:78784/1‑150 (MQ=255)
                                                                                                                          cAGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGAttt                    >  3:744725/1‑150 (MQ=255)
                                                                                                                           aGTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGTACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTc                   >  1:434549/1‑150 (MQ=255)
                                                                                                                            gTAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCt                  <  2:692365/150‑1 (MQ=255)
                                                                                                                             tAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTg                 >  4:117277/1‑150 (MQ=255)
                                                                                                                             tAGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTg                 <  2:347160/150‑1 (MQ=255)
                                                                                                                              aGACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTgg                <  1:1185244/150‑1 (MQ=255)
                                                                                                                               gACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCAc                                                                          <  2:1459960/91‑1 (MQ=255)
                                                                                                                               gACACCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGc               <  2:1047316/150‑1 (MQ=255)
                                                                                                                                acacCGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCGACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCa              >  4:896981/1‑150 (MQ=255)
                                                                                                                                   ccGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt           <  4:129034/150‑1 (MQ=255)
                                                                                                                                   ccGCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACtt           <  1:1078426/150‑1 (MQ=255)
                                                                                                                                    cGCCTGAACACACGCGCCACGCATACCTACGCCAGCACCTTCCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTc          <  1:954662/150‑1 (MQ=255)
                                                                                                                                     gCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCa         <  2:278148/150‑1 (MQ=255)
                                                                                                                                     gCCTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCa         >  3:883837/1‑150 (MQ=255)
                                                                                                                                       cTGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAAt       >  4:449977/1‑150 (MQ=255)
                                                                                                                                        tGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATa      <  1:1066120/150‑1 (MQ=255)
                                                                                                                                        tGAACACACGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATa      <  1:637940/150‑1 (MQ=255)
                                                                                                                                           acacacGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCACTCACTTCTTCGCCGATTTCTGGCACTTCAATACgg   <  4:544090/150‑1 (MQ=255)
                                                                                                                                           acacGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACgg   <  1:357794/150‑1 (MQ=255)
                                                                                                                                            cacaCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGa  >  2:1002354/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CATGGTGTGTTTATCTTCATCCACCACGATAGAAGCAACGTCTGCCGGTGCCATTGCGTTAATCACGAACTGCGCCGGGTTATCATCCCACAGGACGATATCGATACGCTCGCCACCCAGTTCAGTAGACACCGCCTGAACACGCGCGCCACGCATACCTACGCAAGCACCTACCGGATCGATACGTTTATCGTTGGTTTTCACCGCGATTTTCGCACGAGAACCCGGATCGCGAGCCGCTGCTTTAATTTCAATCACTTCTTCGCCGATTTCTGGCACTTCAATACGGA  >  NC_000913_3_pae_pgi/3316612‑3316901

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