Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,127,164 C→G 7.4% V241L (GTG→CTG) ‡ PP_5375 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,127,1640CG 7.4% 257.7 / 8.9 81V241L (GTG→CTG) ‡PP_5375LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (41/34);  new base G (4/2);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.90e-01

AACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCAAGCCAGAGCCAG  >  NC_002947/6127016‑6127312
                                                                                                                                                    |                                                                                                                                                    
aaCAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCAc                                                                                                                                                      >  5:341287/1‑149 (MQ=255)
        cACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGTGGCCCTGGc                                                                                                                                              >  6:33130/1‑149 (MQ=255)
        cACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGTGGCCCTGGc                                                                                                                                              <  5:33130/149‑1 (MQ=255)
        cACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCGTGGCCCTGGc                                                                                                                                              <  6:257547/149‑1 (MQ=255)
        cACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGc                                                                                                                                              <  3:241364/149‑1 (MQ=255)
         aCGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGcc                                                                                                                                             <  4:379896/149‑1 (MQ=255)
           gCTTTCACGAGGACCAGGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGTAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCAc                                                                                                                                           <  3:201640/149‑1 (MQ=255)
                    aGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCg                                                                                                                                  >  3:116877/1‑149 (MQ=255)
                      gACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGc                                                                                                                                              <  1:172911/135‑1 (MQ=255)
                      gACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGc                                                                                                                                              >  2:172911/1‑135 (MQ=255)
                          aCGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCCGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGa                                                                                                                            >  1:356001/1‑149 (MQ=255)
                            gCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCTGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATc                                                                                                                          <  2:103428/149‑1 (MQ=255)
                             cTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCa                                                                                                                         <  6:341287/149‑1 (MQ=255)
                                 gTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCCCTGATCGCCGACTGGATCAGCa                                                                                                                      >  8:325162/1‑148 (MQ=255)
                                  tGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATa                                                                                                                    <  5:266602/149‑1 (MQ=255)
                                  tGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATa                                                                                                                    >  7:94531/1‑149 (MQ=255)
                                   ggACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATAc                                                                                                                   >  7:11913/1‑149 (MQ=255)
                                     acgacTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCTGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCg                                                                                                                 >  7:155709/1‑149 (MQ=255)
                                         cTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGa                                                                                                                                        <  6:17998/122‑1 (MQ=255)
                                         cTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGa                                                                                                                                        >  5:17998/1‑122 (MQ=255)
                                         cTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGc                                                                                                             >  4:54674/1‑149 (MQ=255)
                                          tCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCt                                                                                                            >  6:371484/1‑149 (MQ=255)
                                          tCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCt                                                                                                            >  2:161450/1‑149 (MQ=255)
                                              gcaggcagAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTg                                                                                                        <  1:935/149‑1 (MQ=255)
                                                    agAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCAc                                                                                                                                           <  4:243992/108‑1 (MQ=255)
                                                    agAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCAc                                                                                                                                           >  3:243992/1‑108 (MQ=255)
                                                     gAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGa                                                                                                 >  5:19431/1‑149 (MQ=255)
                                                      aaCTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGAc                                                                                                <  3:54674/149‑1 (MQ=255)
                                                      aaCTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGAc                                                                                                >  3:244179/1‑149 (MQ=255)
                                                        cTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGAcgc                                                                                              <  4:116877/149‑1 (MQ=255)
                                                           ggCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGa                                                                                            <  4:244179/148‑1 (MQ=255)
                                                                 ggTCCGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTg                                                                                      <  2:356001/148‑1 (MQ=255)
                                                                   tCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTgggg                                                                                   <  7:155931/149‑1 (MQ=255)
                                                                      ggAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGt                                                                                 >  3:89854/1‑148 (MQ=255)
                                                                       gAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGGTCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTggg                                                                                    <  1:83092/144‑1 (MQ=255)
                                                                       gAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGGTCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTggg                                                                                    >  2:83092/1‑144 (MQ=255)
                                                                       gAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGTATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGGCa                                                                               >  8:34739/1‑149 (MQ=255)
                                                                            ccACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGc                                                                          <  2:177779/149‑1 (MQ=255)
                                                                             cacaTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTgg                                                                                     <  1:390362/137‑1 (MQ=255)
                                                                             cacaTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTgg                                                                                     >  2:390362/1‑137 (MQ=255)
                                                                             cacaTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCt                                                                         >  7:366140/1‑149 (MQ=255)
                                                                             cacaTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCt                                                                         >  4:346160/1‑149 (MQ=255)
                                                                              acaTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTc                                                                        >  2:167766/1‑149 (MQ=255)
                                                                                aTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGt                                                                                                          >  7:262765/1‑113 (MQ=255)
                                                                                aTCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGt                                                                                                          <  8:262765/113‑1 (MQ=255)
                                                                                   gTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTggg                                                                                    <  7:313127/132‑1 (MQ=255)
                                                                                   catGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGACCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGAACACCATACCGAACCTCTTGATATTGACGCGATCCTTtgg                                                                                    >  8:313127/3‑132 (MQ=255)
                                                                                      gggCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATccc                                                                 <  8:94531/148‑1 (MQ=255)
                                                                                       ggCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCa                                                                <  8:11913/148‑1 (MQ=255)
                                                                                        gCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAgg                                                              >  4:126853/1‑149 (MQ=255)
                                                                                         cTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGgc                                                             >  3:366901/1‑149 (MQ=255)
                                                                                            gACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGgcgc                                                           <  2:36481/148‑1 (MQ=255)
                                                                                              ccAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCtt                                                         >  5:33122/1‑148 (MQ=255)
                                                                                               cAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTg                                                        >  5:182502/1‑148 (MQ=255)
                                                                                                aGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGaa                                                      <  6:84598/149‑1 (MQ=255)
                                                                                                             tCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGGTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCaa                                         <  8:155709/149‑1 (MQ=255)
                                                                                                              cAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAAc                                        >  1:146434/1‑149 (MQ=255)
                                                                                                              cAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAAc                                        >  3:390728/1‑149 (MQ=255)
                                                                                                               aaCGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAAc                                        <  6:19431/148‑1 (MQ=255)
                                                                                                                  gCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTAc                                    <  2:146434/149‑1 (MQ=255)
                                                                                                                        tCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTc                              >  8:202689/1‑149 (MQ=255)
                                                                                                                           accaGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTcgcg                           >  2:85652/1‑149 (MQ=255)
                                                                                                                           accaGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTcgcg                           >  1:40721/1‑149 (MQ=255)
                                                                                                                           accaGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCGCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTcgcg                           >  8:183719/1‑149 (MQ=255)
                                                                                                                             caGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGCGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCt                         >  1:65165/1‑149 (MQ=255)
                                                                                                                               ggTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTgg                       <  2:40721/149‑1 (MQ=255)
                                                                                                                               ggTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTgg                       >  8:26555/1‑149 (MQ=255)
                                                                                                                               ggTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTgg                       <  3:346160/149‑1 (MQ=255)
                                                                                                                                gTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTgg                       >  4:128623/1‑148 (MQ=255)
                                                                                                                                 ttGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGa                     <  3:128623/149‑1 (MQ=255)
                                                                                                                                    ccccAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAgac                  >  4:248735/1‑147 (MQ=255)
                                                                                                                                    ccccAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGt                   <  3:126853/148‑1 (MQ=255)
                                                                                                                                     cccAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCt                 >  1:164546/1‑149 (MQ=255)
                                                                                                                                     cccAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGGTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCt                 >  2:174250/1‑149 (MQ=255)
                                                                                                                                      ccaggccaggGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCtt                <  7:325162/149‑1 (MQ=255)
                                                                                                                                       gtggccaggGCGTGGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCtt                >  2:339351/3‑148 (MQ=255)
                                                                                                                                         ggccaggGCGTGGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTc               >  5:271320/1‑147 (MQ=255)
                                                                                                                                          gccaggGCGTGGCCCTGGCCACTGATCGCCGACTGGATCAGCATAACGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCAAAGGCGCTTGAACCATTCGGGCCAAACTACCCAGTAGCGCCGGGAGTCTTc               >  1:172310/1‑146 (MQ=255)
                                                                                                                                            caggGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCACCCTACCCAGTCGCGCTGGGAGTCtt                >  4:101927/1‑143 (MQ=255)
                                                                                                                                            caggGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCtt                <  3:101927/143‑1 (MQ=255)
                                                                                                                                               ggCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCAagccaga       <  1:161450/149‑1 (MQ=255)
                                                                                                                                                  cACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCAagccagagc     >  6:31274/1‑148 (MQ=255)
                                                                                                                                                    cGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCAagccagagccag  <  4:89854/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACAGCCGCACGCTTTCACGAGGACCACGCTCGGTGGACGACTCCAGCAGGCAGAACTTGGCGTCGGTCTGGAGTACCACATCGTTGGGCTTGACCAGCGCACCGGTGCTCAACGCATCGTCCACCAGGTTGCCCCAGGCCAGGGCCACGCCCTGGCCACTGATCGCCGACTGGATCAGCATACCGTAGCTGTTGATATTGACGCGATGCTTGGGGGTCATGGGCTCGTATCCCAGGCGCTTGAACCATTCGGGCCAACCTACCCAGTCGCGCTGGGAGTCTTCAAGCCAGAGCCAG  >  NC_002947/6127016‑6127312

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