Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,990,887 G→T 5.6% Q262K (CAA→AAA)  PP_3517 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,990,8870GT 5.6% 271.3 / 3.2 88Q262K (CAA→AAA) PP_3517hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/40);  new base T (3/2);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CGCTTGCACAGGGGCCGCAGCTGCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGC  >  NC_002947/3990738‑3991036
                                                                                                                                                     |                                                                                                                                                     
cGCTTGCACAGGGGCCGCAGCTGCGTTGGCTGTGCGGGCGCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTGGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTg                                                                                                                                                       >  1:392461/1‑150 (MQ=255)
 gCTTGCACAGGGGCCGCAGCTGCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTgg                                                                                                                                                      >  1:516520/1‑150 (MQ=255)
  cTTGCACAGGGGCCGCAGCTGCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGc                                                                                                                                                     >  2:771859/1‑150 (MQ=255)
       acaGGGGCCGCAGCTGCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCa                                                                                                                                                <  2:1344453/150‑1 (MQ=255)
         aGGGGCCGCAGCTGCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATg                                                                                                                                              <  1:2201508/150‑1 (MQ=255)
         aGGGGCCGCAGCTGCGTTGGCTGTGCGGGCTCAATCCCCCAACGGCTCACGCACCACGTTCCATCTGTCATACGGATCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATg                                                                                                                                              <  1:534498/150‑1 (MQ=255)
          ggggCCGCAGCTGCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATgg                                                                                                                                             <  1:675849/150‑1 (MQ=255)
           gggCCGCAGCTGCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGc                                                                                                                                            >  1:704862/1‑150 (MQ=255)
                     tGCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGaa                                                                                                                                  <  1:743184/150‑1 (MQ=255)
                      gCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAg                                                                                                                                 <  2:1302319/150‑1 (MQ=255)
                      gCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAg                                                                                                                                 <  2:1525571/150‑1 (MQ=255)
                       cGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTAGCCATGCGGCGCGAGGGCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGc                                                                                                                                >  2:2276838/1‑150 (MQ=255)
                                    ggCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGcc                                                                                                                                                    >  2:1446382/1‑117 (MQ=255)
                                    ggCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGcc                                                                                                                                                    <  1:1446382/117‑1 (MQ=255)
                                    ggCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGgt                                                                                                                   <  1:1167956/150‑1 (MQ=255)
                                       tCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTgtc                                                                                                                >  1:951853/1‑150 (MQ=255)
                                       tCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCACGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTgtc                                                                                                                <  1:103321/150‑1 (MQ=255)
                                         aaTCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTgtcgt                                                                                                              >  2:1308213/1‑150 (MQ=255)
                                         aaTCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTgtcgt                                                                                                              >  2:1114630/1‑150 (MQ=255)
                                           tCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTAGAAGCGTTTGCTGTTCTTCGCCAGCCTTTCTTCAGTCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGtcgtcg                                                                                                            >  1:1486727/1‑150 (MQ=255)
                                             ccccGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGACGACATTGCGTGTCATGAAGCTCATCACCCCGGTGTCGTCGTa                                                                                                          >  2:391106/1‑150 (MQ=255)
                                              cccGCCGCCACATCCACCACTTTGCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTAc                                                                                                         >  2:1196388/1‑150 (MQ=255)
                                              cccGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTAc                                                                                                         <  2:516520/150‑1 (MQ=255)
                                              cccGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTAc                                                                                                         >  1:2419021/1‑150 (MQ=255)
                                               ccgccgCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACt                                                                                                        >  1:1811563/1‑150 (MQ=255)
                                                  acgccACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTgg                                                                                                     >  2:1839932/2‑150 (MQ=255)
                                                    gccACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTggcc                                                                                                   <  2:106505/150‑1 (MQ=255)
                                                         aTCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCa                                                                                              <  1:1308213/150‑1 (MQ=255)
                                                            caccacTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGc                                                                                           >  2:595869/1‑150 (MQ=255)
                                                            caccacGTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGc                                                                                           >  1:1394666/1‑150 (MQ=255)
                                                             accacTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTTGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTgtcgtc                                                                                                             >  1:216527/1‑131 (MQ=255)
                                                             accacTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTTGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTgtcgtc                                                                                                             >  1:216492/1‑131 (MQ=255)
                                                             accacTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTTGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTgtcgtc                                                                                                             <  2:216527/131‑1 (MQ=255)
                                                             accacTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTTGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTgtcgtc                                                                                                             <  2:216492/131‑1 (MQ=255)
                                                              ccacTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCtc                                                                                                                              >  2:2015382/1‑113 (MQ=255)
                                                              ccacTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCtc                                                                                                                              <  1:2015382/113‑1 (MQ=255)
                                                              ccaaTTCCCAGCTCTCATCCGGCTCACAGCGTGTCATTTCGTTCGCCAGCCTGTCGCCATGAGGCTCGAGGTGTTTCTCGAGCACTTGGCCGTCATCGGTGATCATGAAGCTCAGACCCCCGGTGTCGTCGTACTTGGCCGTACAGGCGa                                                                                         <  2:538327/150‑1 (MQ=255)
                                                                  tttCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTATTGAGCACTTGGCCGTCATAGCTAATGCTGGATTATTGGACCCCTTTGTCTTCGTACTTGGCCGGCCAGGCGACCAt                                                                                     >  1:2377704/1‑150 (MQ=255)
                                                                   ttCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATg                                                                                    <  1:595869/150‑1 (MQ=255)
                                                                     ccAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACAGCATCACGTGCTGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGc                                                                                  >  1:42627/1‑150 (MQ=255)
                                                                     ccAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGa                                                                                                                                         <  2:801697/95‑1 (MQ=255)
                                                                     ccAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGa                                                                                                                                         >  1:801697/1‑95 (MQ=255)
                                                                     ccAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGc                                                                                  <  1:1536680/150‑1 (MQ=255)
                                                                        gCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGaa                                                                               <  1:771859/150‑1 (MQ=255)
                                                                          tGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGc                                                                             >  1:2022133/1‑150 (MQ=255)
                                                                               tCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACCCTTGGCAGCCATGGCTGATCATTAATCTCATCACCCCGGTGTCGGCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACgg                                                                        <  1:291434/150‑1 (MQ=255)
                                                                                       gAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTg                                                                >  2:182549/1‑150 (MQ=255)
                                                                                             tttCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTggc                                                          >  2:2124684/1‑148 (MQ=255)
                                                                                               tcgtctcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGa                                                        <  1:2437846/147‑1 (MQ=255)
                                                                                                cATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGAt                                                       >  2:1401515/1‑150 (MQ=255)
                                                                                                 aTCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATc                                                      <  2:1811563/150‑1 (MQ=255)
                                                                                                 aTCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATc                                                      >  1:2255667/1‑150 (MQ=255)
                                                                                                  tctcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCa                                                     >  2:1457320/1‑150 (MQ=255)
                                                                                                   ctcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAg                                                    >  2:2239333/1‑150 (MQ=255)
                                                                                                    tcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAgg                                                   >  2:1657435/1‑150 (MQ=255)
                                                                                                    tcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAgg                                                   >  1:894131/1‑150 (MQ=255)
                                                                                                    tcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAgg                                                   <  1:1657435/150‑1 (MQ=255)
                                                                                                    tcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAgg                                                   >  1:621061/1‑150 (MQ=255)
                                                                                                    tcCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATAAgg                                                   >  2:2388559/1‑150 (MQ=255)
                                                                                                       ttCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCAAAGCTCATCTGGTTGACGAAAAGGTAg                                                >  1:1696156/1‑150 (MQ=255)
                                                                                                            cAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTa                                           <  2:894131/150‑1 (MQ=255)
                                                                                                            cAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTa                                           >  2:1629927/1‑150 (MQ=255)
                                                                                                            cAGCCTGTCGCCATGCGGCCAGAGGTCTTTCTCAAACAGTGGGCCGTCATGTCTTATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTa                                           <  2:1455285/150‑1 (MQ=255)
                                                                                                             aGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCCAACACTTGGCCGTCATGGCTGATCATGAAGCTCCGCGCCCCGGTGTCGTCGGACTTGGCCGGCCAGGCGACGATGGCGACGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAg                                          <  2:1860576/150‑1 (MQ=255)
                                                                                                               ccTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGc                                        >  1:1276762/1‑150 (MQ=255)
                                                                                                                cTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCa                                       <  1:274391/150‑1 (MQ=255)
                                                                                                                   tCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGc                                                                             <  2:772834/109‑1 (MQ=255)
                                                                                                                   tCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGc                                                                             >  1:772834/1‑109 (MQ=255)
                                                                                                                    cGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGa                                                        <  1:1539850/129‑1 (MQ=255)
                                                                                                                    cGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGa                                                        >  2:1539850/1‑129 (MQ=255)
                                                                                                                     gCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGcc                                  >  1:1189431/1‑150 (MQ=255)
                                                                                                                      ccATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGccc                                 >  2:1630013/1‑150 (MQ=255)
                                                                                                                       cATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCg                                >  1:172226/1‑150 (MQ=255)
                                                                                                                              cccGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGCTCCTGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGa                         <  1:391106/150‑1 (MQ=255)
                                                                                                                               ccGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGAt                        <  2:704862/150‑1 (MQ=255)
                                                                                                                                 gAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTTTAGGCACCGCCCGGGGCGGATgg                      <  2:983476/150‑1 (MQ=255)
                                                                                                                                 gAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATgg                      <  1:1629927/150‑1 (MQ=255)
                                                                                                                                    gTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACtt                                                                                                       <  1:169631/66‑1 (MQ=255)
                                                                                                                                    gTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACtt                                                                                                       >  2:169631/1‑66 (MQ=255)
                                                                                                                                    gTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGcc                   <  2:1885670/150‑1 (MQ=255)
                                                                                                                                    gTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGcc                   >  2:1488493/1‑150 (MQ=255)
                                                                                                                                       tttCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTTGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCtg                >  1:1966841/1‑150 (MQ=255)
                                                                                                                                        ttCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCtgt               <  1:1401515/150‑1 (MQ=255)
                                                                                                                                        ttCTCGAACACTTGGCCGGCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCtgt               <  2:2272584/150‑1 (MQ=255)
                                                                                                                                             gAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATc          >  2:1936638/1‑150 (MQ=255)
                                                                                                                                              aaCACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATTAAGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATcc         >  1:626483/1‑150 (MQ=255)
                                                                                                                                              aaCACTTGGCCGTCATGGCTGATCATGAAGCACATCACCCCGGTGTCGTCCTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACTGCCCGGGGCGGATGGGCCCTGTGCATcc         <  2:1250174/150‑1 (MQ=255)
                                                                                                                                                  cTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGa     <  1:1196388/150‑1 (MQ=255)
                                                                                                                                                   ttGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCAc                                      >  2:2222929/1‑116 (MQ=255)
                                                                                                                                                   ttGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCAc                                      <  1:2222929/116‑1 (MQ=255)
                                                                                                                                                     ggCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGc  <  2:951853/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCTTGCACAGGGGCCGCAGCTGCGTTGGCTGTGCGGGCTCAATCCCCCGCCGCCACATCCACCACTTTCCAGCTGTCATCCGGGTCGAAGCGTTTCATCTCCTTCGCCAGCCTGTCGCCATGCGGCCCGAGGTCTTTCTCGAACACTTGGCCGTCATGGCTGATCATGAAGCTCATCACCCCGGTGTCGTCGTACTTGGCCGGCCAGGCGACCATGGCGAAGCCACGGCTCATCTGGTTGCCGATCAGGTAGCTGTAGGCACCGCCCGGGGCGGATGGGCCCTGTGCATCCAGGATGC  >  NC_002947/3990738‑3991036

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