Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,141,475 G→A 5.3% V132M (GTG→ATG)  PP_3645 → L‑proline amide hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,141,4750GA 5.3% 225.1 / 2.4 75V132M (GTG→ATG) PP_3645L‑proline amide hydrolase
Reads supporting (aligned to +/- strand):  ref base G (36/35);  new base A (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

CACCTGCAAATTGCACAAGCCACCGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAG  >  NC_002947/4141328‑4141624
                                                                                                                                                   |                                                                                                                                                     
caccTGCAAATTGCACAAGCCACCGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTg                                                                                                                                                     >  2:1647326/1‑150 (MQ=255)
 accTGCAAATTGCACAAGCCACCGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGa                                                                                                                                                    <  1:1647326/150‑1 (MQ=255)
     gCAAATTGCACAAGCCACCGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGc                                                                                                                                                <  1:926803/150‑1 (MQ=255)
           tGCACAAGCCACCGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCACgcgc                                                                                                                                          <  1:788799/150‑1 (MQ=255)
              acaAGCCACCGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCACCCGCACCCCCGAGCAACGTGCCGGGGTGATCGCCCGCGCTGc                                                                                                                                       >  2:677181/1‑150 (MQ=255)
                   ccACCGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAg                                                                                                                                  >  2:1673026/1‑150 (MQ=255)
                    cACCGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCGGGGGTGATCGCCCGCGCTGCCCAAgc                                                                                                                                 <  1:263159/150‑1 (MQ=255)
                       cGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAgcggc                                                                                                                              <  2:533548/150‑1 (MQ=255)
                       cGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAgcggc                                                                                                                              >  1:180783/1‑150 (MQ=255)
                               gCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCACACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTgg                                                                                                                      <  2:642014/150‑1 (MQ=255)
                                cTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTggg                                                                                                                     <  1:974237/150‑1 (MQ=255)
                                  tctcCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGcc                                                                                                                   >  2:285002/1‑150 (MQ=255)
                                           gcggcCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCa                                                                                                          >  1:1015290/1‑150 (MQ=255)
                                            cggcCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCaa                                                                                                         <  2:1686428/150‑1 (MQ=255)
                                                   gTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGAt                                                                                                  <  1:533777/150‑1 (MQ=255)
                                                        tgatggaTTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGCGCTGGGCCAAGACGCCAACGTCGAtgctgc                                                                                             >  1:67094/7‑150 (MQ=255)
                                                        ccGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGAtgctgt                                                                                             >  1:1898833/1‑149 (MQ=255)
                                                         cGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCATGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGcc                                                                                            >  2:1662050/1‑150 (MQ=255)
                                                               ttCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGAtgc                                                                                                >  1:479174/1‑140 (MQ=255)
                                                               ttCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGAtgc                                                                                                <  2:479174/140‑1 (MQ=255)
                                                                tCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGAcc                                                                                     <  2:1296655/150‑1 (MQ=255)
                                                                  gCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCCACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGa                                                                                                               >  2:1441632/1‑122 (MQ=255)
                                                                  gCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGa                                                                                                               <  1:1441632/122‑1 (MQ=255)
                                                                   cccTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAACGTGGCGGTGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCt                                                                                  <  2:66532/150‑1 (MQ=255)
                                                                   cccTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGc                                                                                                                                       <  1:1445525/97‑1 (MQ=255)
                                                                   cccTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGc                                                                                                                                       >  2:1445525/1‑97 (MQ=255)
                                                                   cccTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCt                                                                                  >  1:1456842/1‑150 (MQ=255)
                                                                   cccTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCATGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCt                                                                                  <  1:1662050/150‑1 (MQ=255)
                                                                    ccTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGACAACGTCGATGCTGCCCTGGACCGCTg                                                                                 <  2:1085573/150‑1 (MQ=255)
                                                                        aaTTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTc                                                                             <  2:1070828/150‑1 (MQ=255)
                                                                          ttACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAg                                                                           >  2:1142542/1‑150 (MQ=255)
                                                                             cccGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGAtg                                                                                                 <  2:560824/125‑1 (MQ=255)
                                                                             cccGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGAtg                                                                                                 >  1:560824/1‑125 (MQ=255)
                                                                              ccGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGt                                                                       >  1:82247/1‑150 (MQ=255)
                                                                              ccGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGt                                                                       <  2:16196/150‑1 (MQ=255)
                                                                               cGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTg                                                                      >  2:1268030/1‑150 (MQ=255)
                                                                                gCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGa                                                                     >  2:76367/1‑150 (MQ=255)
                                                                                  aGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCAGGGGTGATCGCCCGCGCTGCGCAAGCGGTGGGGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAAt                                                                   >  1:924206/1‑150 (MQ=255)
                                                                                     gcCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACa                                                                >  1:1775229/1‑150 (MQ=255)
                                                                                       cTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAg                                                              >  1:123687/1‑150 (MQ=255)
                                                                                        tGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAgg                                                             <  2:305431/150‑1 (MQ=255)
                                                                                          gCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCGCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCt                                                           >  1:1855214/1‑150 (MQ=255)
                                                                                              cACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGt                                                                                                                                                      <  1:1726894/55‑1 (MQ=255)
                                                                                              cACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGt                                                                                                                                                      >  2:1726894/1‑55 (MQ=255)
                                                                                                      gTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCcgcg                                               >  1:143201/1‑150 (MQ=255)
                                                                                                        gCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCa                                             >  1:847618/1‑150 (MQ=255)
                                                                                                         cTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAg                                            >  1:432319/1‑150 (MQ=255)
                                                                                                          tCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAgg                                           >  2:699324/1‑150 (MQ=255)
                                                                                                          tCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAgg                                           >  1:214885/1‑150 (MQ=255)
                                                                                                               aGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTcgcc                                      <  2:1236079/150‑1 (MQ=255)
                                                                                                               aGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTcgcc                                      >  1:1236079/1‑150 (MQ=255)
                                                                                                                gCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAAc                                                     <  2:49410/134‑1 (MQ=255)
                                                                                                                gCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAAc                                                     >  1:49410/1‑134 (MQ=255)
                                                                                                                 cGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTcgccgc                                    <  2:1022915/150‑1 (MQ=255)
                                                                                                                      tCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTc                               <  1:1673026/150‑1 (MQ=255)
                                                                                                                        aaCCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGc                             >  2:1177356/1‑150 (MQ=255)
                                                                                                                        aaCCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGc                             >  2:1177589/1‑150 (MQ=255)
                                                                                                                        aaCCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTAGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGc                             >  1:886520/1‑150 (MQ=255)
                                                                                                                          ccGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGCCACCTGCTCCAGCCGTTAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCa                           <  2:67094/150‑1 (MQ=255)
                                                                                                                          ccGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCa                           <  1:699324/150‑1 (MQ=255)
                                                                                                                          ccGCACCCCCGAGCAAAGTGCCGGGATGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTgg                                                                                >  2:6648/1‑97 (MQ=255)
                                                                                                                          ccGCACCCCCGAGCAAAGTGCCGGGATGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTgg                                                                                >  2:976812/1‑97 (MQ=255)
                                                                                                                          ccGCACCCCCGAGCAAAGTGCCGGGATGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTgg                                                                                <  1:6648/97‑1 (MQ=255)
                                                                                                                          ccGCACCCCCGAGCAAAGTGCCGGGATGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTgg                                                                                <  1:976812/97‑1 (MQ=255)
                                                                                                                            gCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAgg                         <  1:1868414/150‑1 (MQ=255)
                                                                                                                             cACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGt                        >  2:206642/1‑150 (MQ=255)
                                                                                                                                ccccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCa                                                                            >  1:1054394/1‑95 (MQ=255)
                                                                                                                                ccccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCa                                                                            <  2:1054394/95‑1 (MQ=255)
                                                                                                                                ccccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACt                     >  1:384544/1‑150 (MQ=255)
                                                                                                                                ccccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACt                     <  1:1177589/150‑1 (MQ=255)
                                                                                                                                ccccGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACt                     <  1:1177356/150‑1 (MQ=255)
                                                                                                                                   cGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGc                  <  1:206642/150‑1 (MQ=255)
                                                                                                                                     aGCAAAGTGCCGGGGTGATCGCCAGCGCGGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGa                <  2:245225/150‑1 (MQ=255)
                                                                                                                                         aaGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCaa            <  2:143201/150‑1 (MQ=255)
                                                                                                                                              ccGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTAGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGAcc       <  2:886520/150‑1 (MQ=255)
                                                                                                                                                ggggTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGAAGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCATGTACTGGCGAGCAACGACCCg     <  2:712748/150‑1 (MQ=255)
                                                                                                                                                 gggTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCTCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGc    >  1:856099/1‑150 (MQ=255)
                                                                                                                                                  ggTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCa   >  2:1805059/1‑150 (MQ=255)
                                                                                                                                                   gTGATCGCCCGCGCTGCCCAAGCGGCGGGGCTGGGCCCATACGCCAACGTCGATGCTGCCCTCAACTTCGAAATCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAg  <  1:1305000/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CACCTGCAAATTGCACAAGCCACCGTGATCGGCTTCTCCATGGGCGGCCTGGTTGCCCGTGCATTCGCCCTCAATTACCCGCAGCGCCTGGCTGCACTGGTCGTGCTCAACAGCGTGTTCAACCGCACCCCCGAGCAAAGTGCCGGGGTGATCGCCCGCGCTGCCCAAGCGGCGGAGCTGGGCCCAGACGCCAACGTCGATGCTGCCCTGGACCGCTGGTTCAGCCGTGAATACAAGGCTGCCAACCCCGCGCAGGTCGCCGCCATTCGCCAGGTACTGGCGAGCAACGACCCGCAG  >  NC_002947/4141328‑4141624

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