Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,156,036 C→T 7.5% A149A (GCC→GCT)  PP_4539 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,156,0360CT 7.5% 179.1 / 2.2 62A149A (GCC→GCT) PP_4539LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (32/25);  new base T (2/3);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTGGTCGATGCCGCCTACCCCAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCAGGTGGTGATTCG  >  NC_002947/5155890‑5156185
                                                                                                                                                  |                                                                                                                                                     
gtggtCGATGCCGCCTACCCCAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCCGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATc                                                                                                                                                    <  1:712667/150‑1 (MQ=255)
  ggtCGATGCCGCCTACCCCAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATcag                                                                                                                                                  >  1:435135/1‑150 (MQ=255)
    tCGATGCCGCCTACCCCAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCTATcagca                                                                                                                                                >  1:477028/1‑150 (MQ=255)
          ccgccTACCCCAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGGGCGCCTGCGTGAAGAAGTGCTGTCTCGCGTGGAGGAAGTACTGCCCGAAGGCATCGCCGACCTGGCCATCAGCAGCTACa                                                                                                                                          <  2:342160/150‑1 (MQ=255)
           cgccTACCCCAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAg                                                                                                                                         >  2:1126155/1‑150 (MQ=255)
                aCCCCAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATcg                                                                                                                                    >  1:1330554/1‑150 (MQ=255)
                  cccAGGGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCTGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAATGCATCGCCGACCTGGCCATCAGCAGCTACAGCATcggc                                                                                                                                  <  2:239825/150‑1 (MQ=255)
                  cccAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATcggc                                                                                                                                  >  1:938643/1‑150 (MQ=255)
                  cccAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATcggc                                                                                                                                  >  2:153049/1‑150 (MQ=255)
                   ccAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATcggcg                                                                                                                                 <  1:153049/150‑1 (MQ=255)
                    cAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATcggcgg                                                                                                                                <  2:938643/150‑1 (MQ=255)
                    cAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATcggc                                                                                                                                  >  2:1431949/1‑148 (MQ=255)
                    cAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATcggc                                                                                                                                  <  1:1431949/148‑1 (MQ=255)
                         cccgccTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTAcc                                                                                                                           <  1:630758/150‑1 (MQ=255)
                            gccTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGAGGCACGAAGGCATCGCCGACCTGGCCATCAGCAGTTACAGCATCGGCGGCTACCTgg                                                                                                                        <  1:810060/150‑1 (MQ=255)
                            gacTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTgg                                                                                                                        >  1:1373428/3‑150 (MQ=255)
                                    cgcgcgcTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGaa                                                                                                                <  2:367697/150‑1 (MQ=255)
                                    cgcgcgcTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGaa                                                                                                                <  1:761543/150‑1 (MQ=255)
                                        cgcTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCTATCAGCAGCTCTAACAACTGCGGCTACCTGGGCGCCAAATTgt                                                                                                            <  2:477028/150‑2 (MQ=255)
                                        cgcTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTgt                                                                                                            <  2:1124093/150‑2 (MQ=255)
                                        cgcTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGa                                                                                                            <  2:435135/150‑1 (MQ=255)
                                            ggctgCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAgcgc                                                                                                        <  2:1471776/150‑1 (MQ=255)
                                                gcgcTCATGCCACAAAGCCGTGGCGGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATNANCNGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTg                                                                                                    <  1:1126155/150‑1 (MQ=255)
                                                   ttCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCTATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGaa                                                                                                 <  1:40687/150‑1 (MQ=255)
                                                   ttCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGaa                                                                                                 <  2:1330554/150‑1 (MQ=255)
                                                      aTGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTc                                                                                              >  1:178303/1‑150 (MQ=255)
                                                          cacaAAGCCGTGGCTGCCGGGTGCGCCGGCTTGAAGAAGTGCGGGCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCTATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCg                                                                                          >  1:856344/1‑150 (MQ=255)
                                                           acaAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAAt                                                                                                >  1:1229152/1‑143 (MQ=255)
                                                           acaAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAAt                                                                                                <  2:1229152/143‑1 (MQ=255)
                                                             aaaGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGGGGTGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCg                                                                                       >  2:17398/1‑150 (MQ=255)
                                                                  cGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGgc                                                                                                                      >  2:967273/1‑114 (MQ=255)
                                                                  cGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGgc                                                                                                                      <  1:967273/114‑1 (MQ=255)
                                                                   gTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGGGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGccc                                                                                 >  2:764537/1‑150 (MQ=255)
                                                                        tGCCGGGTGAGCAAGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCAcccc                                                                            >  2:100147/1‑150 (MQ=255)
                                                                              gTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGcc                                                                                        <  1:746812/132‑1 (MQ=255)
                                                                              gTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGcc                                                                                        >  2:746812/1‑132 (MQ=255)
                                                                                     tGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGc                                                               >  2:617010/1‑150 (MQ=255)
                                                                                          gaagaaGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGc                                                          >  2:1383030/1‑150 (MQ=255)
                                                                                              aagTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTgg                                                      >  1:650284/1‑150 (MQ=255)
                                                                                                gTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGTCGCCGAATTGa                                                                                                            <  2:483187/94‑1 (MQ=255)
                                                                                                gTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGa                                                                                                            >  1:483187/1‑94 (MQ=255)
                                                                                                        tGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCTATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACt                                            <  2:856344/150‑1 (MQ=255)
                                                                                                             tGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCACCCCCTAACACAGCCTGCACCGCCTGGGACGCGAACTGACCt                                       >  1:1249814/1‑150 (MQ=255)
                                                                                                                 ggaAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCa                                   >  1:339894/1‑150 (MQ=255)
                                                                                                                  gaAGTGATGCACGACGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAg                                  >  1:470679/1‑150 (MQ=255)
                                                                                                                    aGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGc                                                                                   <  2:792747/99‑1 (MQ=255)
                                                                                                                    aGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGc                                                                                   >  1:792747/1‑99 (MQ=255)
                                                                                                                     gTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGAc                               >  2:966511/1‑150 (MQ=255)
                                                                                                                     gTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGAc                               <  1:617010/150‑1 (MQ=255)
                                                                                                                       gATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATACCTCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCt                             >  2:464071/1‑150 (MQ=255)
                                                                                                                        aTGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTg                            >  2:900014/1‑150 (MQ=255)
                                                                                                                          gCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAATGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGTCTGGGACGCGGACGGACCCTCCAGAACCTGGa                          >  2:308956/1‑150 (MQ=255)
                                                                                                                            aCGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATc                                                                                           <  1:1212012/83‑1 (MQ=255)
                                                                                                                            aCGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATc                                                                                           >  2:1212012/1‑83 (MQ=255)
                                                                                                                             cGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGc                                                                   <  2:1428889/106‑1 (MQ=255)
                                                                                                                             cGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGc                                                                   >  1:1428889/1‑106 (MQ=255)
                                                                                                                              gAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGcc                                                                  >  1:978717/1‑106 (MQ=255)
                                                                                                                              gAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGcc                                                                  <  2:978717/106‑1 (MQ=255)
                                                                                                                              gAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGc                      >  1:381497/1‑150 (MQ=255)
                                                                                                                                      cGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCa              <  2:650284/150‑1 (MQ=255)
                                                                                                                                           aCCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCAggtgg         >  2:1385566/1‑150 (MQ=255)
                                                                                                                                                gCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCAGGTGGTGAtt    >  1:980476/1‑150 (MQ=255)
                                                                                                                                                 ccATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCAGGTGGTGATTc   <  2:339894/150‑1 (MQ=255)
                                                                                                                                                 ccATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCAGGTGGTGATTc   >  2:1631699/1‑150 (MQ=255)
                                                                                                                                                 ccATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCAGGTGGTGATTc   <  1:1383030/150‑1 (MQ=255)
                                                                                                                                                  cATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCAGGTGGTGATTCg  >  2:83498/1‑150 (MQ=255)
                                                                                                                                                  cATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCAGGTGGTGATTCg  >  1:242810/1‑150 (MQ=255)
                                                                                                                                                  cATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGAGCGGTGGGTTGCATCGCCGACGCCAACGCCGAACACAGCCTGCACCGCCTGGGGTGCGAAATGACAATCCCGGAACGGGAAAGCCAGCGCCAGGGGGAGAttgg  >  2:684984/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTGGTCGATGCCGCCTACCCCAGCGCCCGCCTGGTGCGCGCGCTGGCTGCGTTCATGCCACAAAGCCGTGGCTGCCGGGTGCGCCTGCGTGAAGAAGTGCTGTCTGGCGTGGAGGAAGTGATGCACGAAGGCATCGCCGACCTGGCCATCAGCAGCTACAGCATCGGCGGCTACCTGGGCGCCGAATTGAGCGCAGTGGAATTCATCGCCGTCGCCCACCCCGAACACAGCCTGCACCGCCTGGGACGCGAACTGACCTTCCAGGACCTGGAAAGCCAGCTGCAGGTGGTGATTCG  >  NC_002947/5155890‑5156185

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