Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,840,288	C→T	7.6%	D411N (GAT→AAT) ‡	PP_3390 ←	porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,840,288	0	C	T	7.6%	194.1 / 4.5	66	D411N (GAT→AAT) ‡	PP_3390	porin
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base T (2/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

CTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGCCAATGCC  >  NC_002947/3840140‑3840421                                                                                                                                                     |                                                                                                                                      cTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATc                                                                                                                                       >  6:117655/1‑149 (MQ=255) cTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATc                                                                                                                                       <  5:248877/149‑1 (MQ=255)  tCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCg                                                                                                                                      >  5:27739/1‑149 (MQ=255)  tCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCg                                                                                                                                      >  7:300854/1‑149 (MQ=255)      gCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGcg                                                                                                                                   >  5:368328/1‑147 (MQ=255)        ccGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCaaaa                                                                                                                                >  4:175939/1‑149 (MQ=255)             aaGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGa                                                                                                                           <  1:19737/149‑1 (MQ=255)             aaGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGa                                                                                                                           <  8:300854/149‑1 (MQ=255)              aGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAg                                                                                                                          <  3:9677/149‑1 (MQ=255)               gCGGCCGTCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAAGCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCTGAATTCGAACGTGTACCCCTCGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAgc                                                                                                                         <  4:257621/149‑1 (MQ=255)                     ggcGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGCGCAATGCGAgcatgcgc                                                                                                                   >  1:381047/1‑143 (MQ=255)                       cGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTc                                                                                                                 <  1:323382/149‑1 (MQ=255)                          aaaaGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCa                                                                                                              >  7:25561/1‑149 (MQ=255)                              gCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgc                                                                                                          >  3:6261/1‑149 (MQ=255)                               ccACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgcg                                                                                                         >  5:242546/1‑149 (MQ=255)                               ccACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgcg                                                                                                         <  1:360213/149‑1 (MQ=255)                               ccACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgcg                                                                                                         <  7:151457/149‑1 (MQ=255)                               ccACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgc                                                                                                          <  6:242546/148‑1 (MQ=255)                                cACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTgcgc                                                                                                        <  8:249319/149‑1 (MQ=255)                                   gCAAAACGCGCAGATCCAAAGCGCAAGCCAAACCGGCTGCGCCCAACGAAGATTGTGTCTCAATTCAGAATTCGAACGAGTACCCCAAGATCACCCGCTTCTCATCCCCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCaa                                                                                                     >  5:353835/1‑149 (MQ=255)                                    cAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAAc                                                                                                    <  5:117655/149‑1 (MQ=255)                                     aaaaCGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACg                                                                                                   >  5:45986/1‑149 (MQ=255)                                      aaaCGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGc                                                                                                  <  3:345525/149‑1 (MQ=255)                                        aCGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAg                                                                                                <  4:6261/149‑1 (MQ=255)                                          gcgcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCAt                                                                                                             >  6:206679/1‑134 (MQ=255)                                          gcgcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCAt                                                                                                             <  5:206679/134‑1 (MQ=255)                                            gcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGCGCAATGCGAGCATGCGCTCGAATTTTGCCGATGCAGCgga                                                                                            <  6:368328/149‑3 (MQ=21)                                            gcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGCGCAATGCGAGCATGCGCTCGAATTTTGCCGATGCAGCgga                                                                                            <  2:381047/149‑3 (MQ=21)                                            gcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGac                                                                                            <  3:26959/149‑1 (MQ=255)                                            gcAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGa                                                                                             >  4:302602/1‑148 (MQ=255)                                               gATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTTGGCGCAATGCGAGCATGCGCTCGAATTTTGCCGATGCAGCggatga                                                                                         <  7:316475/149‑6 (MQ=21)                                                aTCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACAt                                                                                        <  3:302602/149‑1 (MQ=255)                                                    aaaGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGCCACATTcc                                                                                     >  1:347236/1‑148 (MQ=255)                                                            aGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGccc                                                                            <  5:308344/149‑1 (MQ=255)                                                             gCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGcccc                                                                           >  4:382533/1‑149 (MQ=255)                                                              cGAAACCGGCTGCGCCCAACGAAGATTGTGGCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGccccc                                                                          <  8:306939/149‑1 (MQ=255)                                                               gAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGccccc                                                                          >  1:65440/1‑148 (MQ=255)                                                                   ccGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTg                                                                     >  5:407777/1‑149 (MQ=255)                                                                           gcCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGATCGCCCCCGACTGCACcaca                                                              >  6:50023/1‑148 (MQ=255)                                                                               aaCGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGca                                                                   >  5:146430/1‑139 (MQ=255)                                                                               aaCGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGca                                                                   <  6:146430/139‑1 (MQ=255)                                                                               aaCGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGt                                                         <  1:272329/149‑1 (MQ=255)                                                                                aCGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCGTCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTg                                                        <  6:45986/149‑1 (MQ=255)                                                                                   aaGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGata                                                     >  7:344623/1‑149 (MQ=255)                                                                                       ttGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCg                                                 <  6:27739/149‑1 (MQ=255)                                                                                          tgtCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTa                                               <  2:65440/148‑1 (MQ=255)                                                                                              tGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTg                                                                     >  8:363246/1‑122 (MQ=255)                                                                                               gAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGc                                         >  2:70190/1‑149 (MQ=255)                                                                                                  ttCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTcc                                                                                     <  8:154620/102‑1 (MQ=255)                                                                                                  ttCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTcc                                                                                     >  7:154620/1‑102 (MQ=255)                                                                                                  ttCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTcc                                      >  7:123213/1‑149 (MQ=255)                                                                                                  ttCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTAGCATTGCGCCAACGAAGCGACACATTACTGAGCGCACCCGAATGAACCACATAGGTGATATCCGTAGTCCGCTcc                                      >  1:14580/1‑149 (MQ=255)                                                                                                            gAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGa                                                                                  >  6:366645/1‑95 (MQ=255)                                                                                                            gAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGa                                                                                  <  5:366645/95‑1 (MQ=255)                                                                                                             aaCGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGccc                           <  6:179000/149‑1 (MQ=255)                                                                                                                 tgtACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCctct                       <  3:175939/149‑1 (MQ=255)                                                                                                                    aCCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTACACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCg                    >  2:357900/1‑149 (MQ=255)                                                                                                                      cccACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTc                  >  7:270807/1‑149 (MQ=255)                                                                                                                      cccACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTACACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTc                  <  1:357900/149‑1 (MQ=255)                                                                                                                       ccACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCg                 >  2:403468/1‑149 (MQ=255)                                                                                                                       ccACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCg                 >  1:380871/1‑149 (MQ=255)                                                                                                                        cACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGc                <  2:380871/149‑1 (MQ=255)                                                                                                                                tctggTGCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCTCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGc         <  7:256880/145‑1 (MQ=255)                                                                                                                                cccGGTTCTCATCAGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGc         >  5:173543/1‑148 (MQ=255)                                                                                                                                 aaggTTCTCAACCGCTTCATCGGCAAAATTAGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTACTGAGCGCCCCCGACTGCACCAAATAGGTGAAATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGcc        >  2:236031/3‑148 (MQ=255)                                                                                                                                  cGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGCCaa      >  1:160203/1‑149 (MQ=255)                                                                                                                                    gTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCgac         >  8:256880/1‑142 (MQ=255)                                                                                                                                       ctcATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGCCAAtgcc  <  6:407777/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                      CTCAAGCCCGGCAAGCGGCCGGCGAAAAAGCCACGCAAAACGCGCAGATCCAAAGCGCAAGCGAAACCGGCTGCGCCCAACGAAGATTGTGTCTGAATTCAGAATTCGAACGTGTACCCCACGATCACCCGGTTCTCATCCGCTGCATCGGCAAAATTCGAGCGCATGCTCGCATTGCGCCAACGCAGCGACACATTCCTGAGCGCCCCCGACTGCACCACATAGGTGATATCCGTAGTCCGCTCCCATTCCCGCCCCTCTCCGTCGCCAGCGGCCAATGCC  >  NC_002947/3840140‑3840421

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