Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,259,759	Δ1 bp	6.7%	coding (347/2736 nt)	fimV →	motility protein FimV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,259,759	0	G	.	6.7%	190.5 / 6.9	60	coding (347/2736 nt)	fimV	motility protein FimV
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base . (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATT  >  NC_002947/2259615‑2259907                                                                                                                                                  |                                                                                                                                                     gCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCg                                                                                                                                                  >  2:147139/1‑149 (MQ=255)      gCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGCGAGAAGTGTGATGAGCGTGTCCTCACGTGTGACGCTGGCCGGTGCGCTGGTGAAGTTTATGGTGGTCGGGATGTGGCCGCact                                                                                                                                              >  1:261118/1‑147 (MQ=255)          gggTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGgccga                                                                                                                                         <  7:257955/149‑2 (MQ=255)            gTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCct                                                                                                                                       >  5:171520/1‑149 (MQ=255)                 gTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGAGCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCGTGGTGCAGGTGATGT‑GCCGCAGGGCCGCCTGCTGCg                                                                                                                                 >  8:207675/1‑149 (MQ=255)                 gTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGT‑GCCGCAGGGCCGCCTGCTGCg                                                                                                                                 >  8:207671/1‑149 (MQ=255)                   ttCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGt                                                                                                                                >  7:255414/1‑149 (MQ=255)                    tCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGt                                                                                                                                <  6:213862/148‑1 (MQ=255)                      cGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGa                                                                                                                              <  1:51110/148‑1 (MQ=255)                       gACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACt                                                                                                                            <  2:15359/149‑1 (MQ=255)                          ctacctCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGa                                                                                                                              <  4:298220/144‑1 (MQ=255)                          ctacctCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGa                                                                                                                              >  3:298220/1‑144 (MQ=255)                          ctacctCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACt                                                                                                                            <  7:207053/146‑1 (MQ=255)                          ctacctCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACt                                                                                                                            >  8:207053/1‑146 (MQ=255)                                 gAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGtgct                                                                                                                  >  3:37428/1‑149 (MQ=255)                                        tgacctTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGAc                                                                                                           >  4:204703/1‑149 (MQ=255)                                          acctTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGAc                                                                                                           >  6:271578/1‑147 (MQ=255)                                               cacaCCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGcc                                                                                                     >  4:234186/1‑148 (MQ=255)                                                 cacCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAg                                                                                                  <  2:202041/149‑1 (MQ=255)                                                  acCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAgg                                                                                                 <  6:112367/149‑1 (MQ=255)                                                      gTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACAATAGCGTGCTGCTCGACCAGGCCAAGGca                                                                                               >  5:219905/1‑147 (MQ=255)                                                         aTCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAggg                                                                                          >  2:184608/1‑149 (MQ=255)                                                         aTCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAgg                                                                                           <  6:171520/148‑1 (MQ=255)                                                          tCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGTTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGc                                                                                         >  1:235053/1‑149 (MQ=255)                                                              ccccAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGaa                                                                                      >  7:56456/1‑148 (MQ=255)                                                              ccccAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTAGACCAGGCCAAGGCACAGGGCGaa                                                                                      >  7:52074/1‑148 (MQ=255)                                                                ccAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAg                                                                                   >  8:27929/1‑149 (MQ=255)                                                                ccAACGGC‑AAAAGTGTGCTGAGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCAGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAg                                                                                   >  7:97934/1‑149 (MQ=255)                                                                       c‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAg                                                                                                  >  5:10489/1‑127 (MQ=255)                                                                       c‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAg                                                                                                  <  6:10489/127‑1 (MQ=255)                                                                          aaaGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCggca                                                                          <  1:184608/149‑1 (MQ=255)                                                                             gtgtGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACAATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCggcaggc                                                                       <  6:219905/149‑1 (MQ=255)                                                                                tgctgcGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAAt                                                                    <  4:37428/149‑1 (MQ=255)                                                                                    gcgcGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTca                                                                <  3:147569/149‑1 (MQ=255)                                                                                        gTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACAcc                                                            <  3:234186/149‑1 (MQ=255)                                                                                          gACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAgg                                                                                                       >  5:224089/1‑104 (MQ=255)                                                                                          gACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAgg                                                                                                       <  6:224089/104‑1 (MQ=255)                                                                                             cTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCa                                                       <  8:255414/149‑1 (MQ=255)                                                                                             cTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGT‑GCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCAt                                                      <  7:207671/149‑1 (MQ=255)                                                                                             cTCCAGACAGCCGCTGCCAGGTCCGGTGGTCAATTTCCTTGTGCAGGTGATGT‑GCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCTACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCac                                                      <  7:207675/149‑2 (MQ=255)                                                                                                   tCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCg                                                 >  3:131934/1‑149 (MQ=255)                                                                                                     aGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGt                                               >  8:30184/1‑149 (MQ=255)                                                                                                         gCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCca                                           <  1:42347/149‑1 (MQ=255)                                                                                                            gCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCcag                                          >  4:286113/1‑147 (MQ=255)                                                                                                             ccAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGTGATGGCCGGTgtga                                           >  1:306162/1‑142 (MQ=255)                                                                                                             ccAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGTGATGGCCGGTgtga                                           <  2:306162/145‑4 (MQ=255)                                                                                                               aGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCgg                                                          <  6:146526/128‑1 (MQ=255)                                                                                                               aGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCgg                                                          >  5:146526/1‑128 (MQ=255)                                                                                                                         gtCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACaccacc                             >  6:206482/1‑147 (MQ=255)                                                                                                                         gtCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCa                           <  7:27929/149‑1 (MQ=255)                                                                                                                          tCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAg                          >  5:193802/1‑149 (MQ=255)                                                                                                                           cAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGc                         >  2:230895/1‑149 (MQ=255)                                                                                                                            aaGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCg                        >  5:3190/1‑149 (MQ=255)                                                                                                                                 ccTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTcgc                    <  5:271578/148‑1 (MQ=255)                                                                                                                                  cTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGa                  >  6:291299/1‑149 (MQ=255)                                                                                                                                          ggTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCTTCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTtgt          <  1:230895/149‑1 (MQ=255)                                                                                                                                          ggTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTtgt          <  3:286113/149‑1 (MQ=255)                                                                                                                                           gTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTtgtg         <  8:52074/149‑1 (MQ=255)                                                                                                                                           gTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTtgtg         <  8:56456/149‑1 (MQ=255)                                                                                                                                           gTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACAGGACATCACCGGTGCCAGCACCTACACCACCCAGCGTCGCGACACCTtatg         >  8:136406/1‑149 (MQ=255)                                                                                                                                                gtgGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTgg        >  5:122019/1‑145 (MQ=255)                                                                                                                                                  ggCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGAtt  >  2:77989/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GCAAGGCCGGGGTGGCGTTTCCGACCTACCTCGAAGACCTGACCTTCACACCGGTGATCAACCCCAACGGC‑AAAAGTGTGCTGCGCGTGACCTCCAGTCAGCCGCTGCCAGGTCCGGTGGTCAAGTTCCTGGTGCAGGTGATGTGGCCGCAGGGCCGCCTGCTGCGTGACTATAGCGTGCTGCTCGACCAGGCCAAGGCACAGGGCGAACAGCCGGCGGCAGGCAATGTCACACCGGCCATCACCGGTGCCAGCAGCTACACCACCCAGCGTCGCGACACCTTGTGGCAGATT  >  NC_002947/2259615‑2259907

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