Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,976,324 C→G 9.0% intergenic (+109/‑56) pheA → / → PP_1770 bifunctional chorismate mutase/prephenate dehydratase/3‑phosphoshikimate 1‑carboxyvinyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,976,3240CG 9.0% 280.8 / 15.6 89intergenic (+109/‑56)pheA/PP_1770bifunctional chorismate mutase/prephenate dehydratase/3‑phosphoshikimate 1‑carboxyvinyltransferase
Reads supporting (aligned to +/- strand):  ref base C (38/43);  new base G (5/3);  total (43/46)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCTCGTTCGCCAAGGGC  >  NC_002947/1976181‑1976472
                                                                                                                                               |                                                                                                                                                    
aGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGc                                                                                                                                                <  1:147230/150‑1 (MQ=255)
  gTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGccc                                                                                                                                              <  1:920229/150‑1 (MQ=255)
      tGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAAt                                                                                                                                          <  2:771471/150‑1 (MQ=255)
       gggCTCTTATCTGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATg                                                                                                                                         <  1:772184/150‑1 (MQ=255)
       gggCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCTTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATg                                                                                                                                         <  1:46650/150‑1 (MQ=255)
       gggCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGCGCTGCCTTGCAGCCCTTCGCGGGCATACGCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCAGCGAAGGGCCGCGCCGCGGCCCCAATg                                                                                                                                         <  1:165496/150‑1 (MQ=255)
        ggCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCCTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATgg                                                                                                                                        <  2:952361/150‑1 (MQ=255)
         gCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGc                                                                                                                                       <  1:666561/150‑1 (MQ=255)
         gCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCACTTCGCGGGCATGCCCGCTCCCACAGTGAGAATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCGGCGCCGCGGCCCCAATGGc                                                                                                                                       >  2:892001/1‑150 (MQ=255)
         gCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTGGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGGCGCGGcccttcgcgg                                                                                                                                       <  1:80994/150‑8 (MQ=255)
          gatctATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGAAGCCCTTCGCGGGCATGCCCGCTCCCCCAGTGAGCATGTGATGACGATGGTATTTGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGcc                                                                                                                                      <  1:1106240/146‑1 (MQ=255)
           tctTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATTGccc                                                                                                                                     >  1:1008990/1‑150 (MQ=255)
                tCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcaca                                                                                                                                <  2:1206619/150‑1 (MQ=255)
                tCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTcaca                                                                                                                                <  2:638182/150‑1 (MQ=255)
                 ccGCAGGCGCTGCTTTGATCCCGGCGCCTCGGGGCTGCTGGTCAGCCCTTCGCGGGCATGCCCGCTTCCACAGTGGGCATGTGACGCCGATCTCTGCGGGATCGGGCTAGCCCGCGACAGGCCGCGCCGCGGCCCCAATGGCCCTCACAg                                                                                                                               <  1:654642/150‑1 (MQ=255)
                   gAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGc                                                                                                                             >  2:621445/1‑150 (MQ=255)
                       gCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTAcc                                                                                                                         >  1:338995/1‑150 (MQ=255)
                            gCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAt                                                                                                                    >  2:159557/1‑150 (MQ=255)
                              tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGGCGCGGccctt                                                                                                                                            >  1:479774/1‑122 (MQ=255)
                              tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGGCGCGGccctt                                                                                                                                            <  2:479774/124‑3 (MQ=255)
                              tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAttt                                                                                                                  >  1:57197/1‑150 (MQ=255)
                              tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAttt                                                                                                                  >  2:1366987/1‑150 (MQ=255)
                              tttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGAttt                                                                                                                  <  1:334895/150‑1 (MQ=255)
                               ttGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTc                                                                                                                 <  1:1334429/150‑1 (MQ=255)
                                tGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCGCCCACAGTGAGCATGTGATGCCGATCCCTGTGGGAGCGGGCCAGCCCGCGAAGGGTCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCt                                                                                                                <  2:1031497/150‑1 (MQ=255)
                                      aGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAg                                                                                                          <  2:30751/150‑1 (MQ=255)
                                      aGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAg                                                                                                          >  2:1352121/1‑150 (MQ=255)
                                       ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                                                                         <  2:445034/150‑1 (MQ=255)
                                       ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                                                                         <  1:514007/150‑1 (MQ=255)
                                       ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAg                                                                                                          >  2:200715/1‑149 (MQ=255)
                                       ggggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAg                                                                                                          <  1:200715/149‑1 (MQ=255)
                                        gggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAggg                                                                                                        >  1:472252/1‑150 (MQ=255)
                                        gggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAggg                                                                                                        >  2:245805/1‑150 (MQ=255)
                                         ggCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGt                                                                                                       <  2:369483/150‑1 (MQ=255)
                                         ggCCTTGGGGCTGCCTTGCAGCCCTTCGCGCGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGt                                                                                                       >  1:422692/1‑150 (MQ=255)
                                            cTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGGCGCGGccc                                                                                                                                              >  2:924855/1‑108 (MQ=255)
                                            cTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGGCGCGGccc                                                                                                                                              <  1:924855/108‑1 (MQ=255)
                                             ttGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACCGTGAGCATGTTATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGc                                                                                                   >  1:123845/1‑150 (MQ=255)
                                                gggCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTAc                                                                                                                          >  2:1158071/1‑124 (MQ=255)
                                                gggCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTAc                                                                                                                          <  1:1158071/124‑1 (MQ=255)
                                                 tgCTGCCTGGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCa                                                                                               <  1:1092347/149‑1 (MQ=255)
                                                 ggCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCa                                                                                               <  2:57197/150‑1 (MQ=255)
                                                  gCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAg                                                                                              <  1:1178370/150‑1 (MQ=255)
                                                   cTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGt                                                                                             <  2:422692/150‑1 (MQ=255)
                                                   cTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGt                                                                                             >  2:502154/1‑150 (MQ=255)
                                                    tGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGc                                                                                                                                       >  2:785820/1‑107 (MQ=255)
                                                    tGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGc                                                                                                                                       <  1:785820/107‑1 (MQ=255)
                                                     gCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTgg                                                                                           >  1:445735/1‑150 (MQ=255)
                                                        ttGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTaa                                                                                        >  1:1273908/1‑150 (MQ=255)
                                                          gCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAAt                                                                                      >  2:1024051/1‑150 (MQ=255)
                                                           cAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATg                                                                                     >  2:177251/1‑150 (MQ=255)
                                                           cAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATg                                                                                     <  2:333265/150‑1 (MQ=255)
                                                              cccTTCGCGGGCATGCCGGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATgcag                                                                                  >  1:1116981/1‑150 (MQ=255)
                                                              cccTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATgcag                                                                                  <  2:1273908/150‑1 (MQ=255)
                                                                          aTGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAaaaccaaac                                                                      >  1:928910/1‑150 (MQ=255)
                                                                              ccGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCcgc                                                                  <  2:338995/150‑1 (MQ=255)
                                                                                  tCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGccc                                                              <  1:1024051/150‑1 (MQ=255)
                                                                                   cccACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCca                                                             <  1:159557/150‑1 (MQ=255)
                                                                                     cacaGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcatc                                                           >  2:1170811/1‑150 (MQ=255)
                                                                                     cacaGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCcatc                                                           >  1:1211873/1‑150 (MQ=255)
                                                                                      acaGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGTAAAAACCAAACCCGCGCCcatca                                                          >  2:228202/1‑150 (MQ=255)
                                                                                        aGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCatcatc                                                        <  1:841092/150‑1 (MQ=255)
                                                                                           gAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGAc                                                     <  1:245805/150‑1 (MQ=255)
                                                                                           gAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCAGCGCCCATCATCGAc                                                     >  1:733925/1‑150 (MQ=255)
                                                                                                tgtgATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGCAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCt                                                >  1:1152980/1‑150 (MQ=255)
                                                                                                 gtgATGCCGATCTCTTTTGGAGCTAGCTAGCCCGCGAACGGCCGCGCCGCCGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTTGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtg                                               <  1:892001/150‑1 (MQ=255)
                                                                                                  tgATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                                                                         <  2:1368252/91‑1 (MQ=255)
                                                                                                  tgATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAgg                                                                                                         >  1:1368252/1‑91 (MQ=255)
                                                                                                    aTGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTaaa                                                                                       <  2:1392001/107‑1 (MQ=255)
                                                                                                    aTGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtggtc                                            <  1:1170811/150‑1 (MQ=255)
                                                                                                    aTGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAAGGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTaaa                                                                                       >  1:1392001/1‑107 (MQ=255)
                                                                                                       atgctcactGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtggtcgtg                                         <  1:1282550/143‑1 (MQ=255)
                                                                                                        cGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCtggtcgtgg                                        >  1:934053/1‑150 (MQ=255)
                                                                                                              ctGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAaaaccaaa                                                                       <  2:522217/113‑1 (MQ=255)
                                                                                                              ctGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAaaaccaaa                                                                       >  1:522217/1‑113 (MQ=255)
                                                                                                               tgtgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCt                                 >  2:995795/1‑150 (MQ=255)
                                                                                                               tgtgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCt                                 >  2:363099/1‑150 (MQ=255)
                                                                                                                 tgGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCCACCGACTGGTCGTGGTCGGTCTaa                               >  2:298240/1‑149 (MQ=255)
                                                                                                                  gggAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAgg                              <  2:445735/150‑1 (MQ=255)
                                                                                                                    gAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCcgcg                                                                 <  1:795141/113‑1 (MQ=255)
                                                                                                                    gAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCcgcg                                                                 >  2:795141/1‑113 (MQ=255)
                                                                                                                     aGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCt                           <  1:927013/150‑1 (MQ=255)
                                                                                                                          gCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATcg                      >  1:38841/1‑150 (MQ=255)
                                                                                                                              gCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATcggcgg                  <  1:494095/150‑1 (MQ=255)
                                                                                                                                ccGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCt                >  1:601643/1‑150 (MQ=255)
                                                                                                                                 cgcgAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCTc               >  2:1128625/1‑150 (MQ=255)
                                                                                                                                  cattggggCCGCGGCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCtgg              >  2:865020/6‑148 (MQ=255)
                                                                                                                                      ggggCCGCGGCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCTCGTTCg          >  1:1359907/2‑150 (MQ=255)
                                                                                                                                      ggggCCGCGGCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCTCGTTCg          >  1:1360205/2‑150 (MQ=255)
                                                                                                                                          ccgcgccgcgGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCTCGTTCGCCaa      <  2:123845/150‑1 (MQ=255)
                                                                                                                                            gcgccgcgGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCTCGTTCGCCAAgg    <  2:1152980/150‑1 (MQ=255)
                                                                                                                                            gcgccgcgGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCTCGTTCGCCAAgg    <  1:1128625/150‑1 (MQ=255)
                                                                                                                                              gccgcgGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCTCGTTCGCCAAGGGc  <  1:177251/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AGGTGCTGGGCTCTTATCCGAAGGCGGTGCTTTGATTCAGGGGCCTTGGGGCTGCCTTGCAGCCCTTCGCGGGCATGCCCGCTCCCACAGTGAGCATGTGATGCCGATCTCTGTGGGAGCGGGCTAGCCCGCGAAGGGCCGCGCCGCGGCCCCAATGGCCCTCACAGGCTACCTAGATTTCTGAACAGGGTTCGCACCAGTGGTAAATGCAGCAAAAACCAAACCCGCGCCCATCATCGACCGCCTGGTCGTGGTCGGTCTAGGCCTGATCGGCGGCTCGTTCGCCAAGGGC  >  NC_002947/1976181‑1976472

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