Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,076,095 T→A 6.0% N156I (AAT→ATT)  PP_5328 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,076,0950TA 6.0% 206.5 / 3.2 67N156I (AAT→ATT) PP_5328phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (31/32);  new base A (2/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

ACGCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGT  >  NC_002947/6075951‑6076236
                                                                                                                                                  |                                                                                                                                             
acgCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGa                                                                                                                                           <  2:206359/149‑1 (MQ=255)
 cgCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGAc                                                                                                                                          <  4:311557/149‑1 (MQ=255)
  gCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACa                                                                                                                                         <  7:291071/149‑1 (MQ=255)
   cTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGCGAAAACCAGCGCCTGGCCATTGGACa                                                                                                                                         <  1:3938/149‑1 (MQ=255)
        gTGGTCCAGCGCACGCCCCTGCTAGTCGAGCACGAACGATTGATCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGAAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACATGccaa                                                                                                                                   >  8:83322/1‑147 (MQ=255)
              cAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcacc                                                                                                                              <  1:3993/148‑1 (MQ=255)
                      gCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                            <  8:216496/142‑1 (MQ=255)
                      gCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                            >  7:216496/1‑142 (MQ=255)
                      gCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAgcgcg                                                                                                                     >  8:147950/1‑149 (MQ=255)
                      gCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAgcgcg                                                                                                                     <  8:60108/149‑1 (MQ=255)
                        cccTGCTCGTCGAGCACGAACGATTGACTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAgcgcgc                                                                                                                    >  1:4286/1‑149 (MQ=255)
                        cccTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCt                                                                                                                   >  2:280501/1‑149 (MQ=255)
                              tcgtcgAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCggg                                                                                                             >  4:285892/1‑149 (MQ=255)
                                gtcgAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCt                                                                                                           <  7:147950/149‑1 (MQ=255)
                                     gCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGt                                                                                                        >  8:174443/1‑147 (MQ=255)
                                        cgaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTa                                                                                                   <  7:174443/149‑1 (MQ=255)
                                         gaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAGGcgc                                                                                                                                    >  6:259203/1‑113 (MQ=255)
                                         gaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAGGcgc                                                                                                                                    <  5:259203/115‑3 (MQ=255)
                                         gaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTAc                                                                                                  >  2:11531/1‑149 (MQ=255)
                                         gaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGATGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACAACCGGCGCGCTGCCCGGCTGGTCAGTAc                                                                                                  >  7:331220/1‑149 (MQ=255)
                                                      gCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGt                                                                                        <  2:289702/147‑1 (MQ=255)
                                                               gTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGgtgt                                                                             >  4:286512/1‑149 (MQ=255)
                                                                   tCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGATCATATTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGCCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACTTGACCTGGGTGTCTGc                                                                         >  3:58559/1‑149 (MQ=255)
                                                                    cGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGAc                                                                                          >  3:345519/1‑131 (MQ=255)
                                                                    cGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGAc                                                                                          <  4:345519/131‑1 (MQ=255)
                                                                     gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGa                                                                                               >  7:217027/1‑125 (MQ=255)
                                                                     gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGa                                                                                               <  8:217027/125‑1 (MQ=255)
                                                                        gCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCa                                                                    >  8:167149/1‑149 (MQ=255)
                                                                            gggCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACggg                                                                <  5:37785/149‑1 (MQ=255)
                                                                            gggCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACggg                                                                <  1:280501/149‑1 (MQ=255)
                                                                              gCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCCGGGTGTCTGCCGGCACGGGTa                                                              <  1:11531/149‑1 (MQ=255)
                                                                                tGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGt                                                            >  8:59484/1‑149 (MQ=255)
                                                                                  aTGGTTTTCTTGTTGTCCGTGTAGGTGATCTTGTAGGGGTGCTGGAAAATCAGCGACGGGCCATTGGACAGGCCCAGCACCACCAGCGCGCGGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGtt                                                           <  8:331220/148‑1 (MQ=255)
                                                                                   tGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTc                                                          >  2:144708/1‑148 (MQ=255)
                                                                                     gTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGt                                                                                        <  5:47708/116‑1 (MQ=255)
                                                                                     gTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGt                                                                                        >  6:47708/1‑116 (MQ=255)
                                                                                       tttCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGacac                                                     >  8:301309/1‑149 (MQ=255)
                                                                                             gttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCtt                                                 >  6:246790/1‑147 (MQ=255)
                                                                                             gttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGa                                               >  3:95989/1‑149 (MQ=255)
                                                                                              ttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACgg                                                                 <  7:337489/130‑1 (MQ=255)
                                                                                              ttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACgg                                                                 >  8:337489/1‑130 (MQ=255)
                                                                                                gtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGtt                                            >  7:43139/1‑149 (MQ=255)
                                                                                                 tCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTc                                           <  7:59484/149‑1 (MQ=255)
                                                                                                         ggTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCg                                    <  7:301309/148‑1 (MQ=255)
                                                                                                          gTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGca                                                                                                                                <  8:109537/55‑1 (MQ=255)
                                                                                                          gTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGca                                                                                                                                >  7:109537/1‑55 (MQ=255)
                                                                                                                  gTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCCGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTGTGCCGGGACGGGTAGTTCGACACGCGTGAGTTCGTTACCGGTCTTGGTag                          >  4:41903/1‑148 (MQ=255)
                                                                                                                  gTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTg                              <  2:4286/145‑1 (MQ=255)
                                                                                                                    aGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGATGGTCAGTACTGATGGACGTGACCTGGGTGTATGCAGGCACTGGTAGTTCGAAACGCTTGAGTTCGTTAACGGTCGTGATATc                         >  8:331437/1‑148 (MQ=255)
                                                                                                                     ggTGTGCTGGAAAAGCAGCGCCTGGCCATTGGACAGGCCCAGCACCAGCAGCGCGCTGCCCGGGTGGTCAGTACTTATGGACGTGACCTGTGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCg                        <  4:58559/148‑1 (MQ=255)
                                                                                                                     ggTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCg                        <  1:144708/148‑1 (MQ=255)
                                                                                                                      gtgtGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCgaa                      <  8:43139/149‑1 (MQ=255)
                                                                                                                            tGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACaa                >  4:301755/1‑149 (MQ=255)
                                                                                                                             ggAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTc                                                          <  4:192221/106‑1 (MQ=255)
                                                                                                                             ggAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTc                                                          >  3:192221/1‑106 (MQ=255)
                                                                                                                               gAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATg              >  8:202713/1‑149 (MQ=255)
                                                                                                                                aaaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTAcc                                     >  4:59081/1‑125 (MQ=255)
                                                                                                                                aaaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTAcc                                     <  3:59081/125‑1 (MQ=255)
                                                                                                                                aaaaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGc             <  5:246790/149‑1 (MQ=255)
                                                                                                                                  aaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCTGTCTTGGTATCGAAGAACAATGCCt           <  3:286512/149‑1 (MQ=255)
                                                                                                                                    ccAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCgaag                     <  4:105879/137‑1 (MQ=255)
                                                                                                                                    ccAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCgaag                     >  3:105879/1‑137 (MQ=255)
                                                                                                                                     cAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGc        >  5:301412/1‑149 (MQ=255)
                                                                                                                                     cAGCGCCTGGACATTTGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGc        <  7:331437/149‑1 (MQ=255)
                                                                                                                                       gggcCTGTCCAATGGCCAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTa                           >  4:315014/3‑128 (MQ=255)
                                                                                                                                       gggcCTGTCCAATGGCCAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTa                           <  3:315014/126‑1 (MQ=255)
                                                                                                                                           cTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGt  <  7:167149/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
ACGCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGT  >  NC_002947/6075951‑6076236

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: