Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,166,935 2 bp→CG coding (312‑313/483 nt) b3023 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,166,9350GC99.7% 1212.5 / NA 374L105V (CTT→GTT) b3023hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (174/199);  minor base A (0/1);  total (174/200)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0009133,166,9360CG100.0% 1217.0 / NA 375Q104H (CAG→CACb3023hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (173/202);  total (173/202)

Only 100 of 388 total aligned reads displayed.

GTGCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAAGCTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGCTGCC  >  NC_000913/3166848‑3167022
                                                                                       ||                                                                                      
gTGCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGc                                                                                    >  2:3968731/1‑93 (MQ=255)
  gCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                                  <  1:4928388/93‑1 (MQ=255)
  gCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                                  <  2:5974431/93‑1 (MQ=255)
  gCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                                  <  2:636017/93‑1 (MQ=255)
  gCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCgg                                                                                  >  1:2965930/1‑93 (MQ=255)
    tttACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGc                                                                                >  2:4957966/1‑93 (MQ=255)
     ttACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGcc                                                                               <  2:7585947/93‑1 (MQ=255)
      tACGCTTTTTCTCGCCGCTGGCAGGCAGCAAGTGGCGAATGATGCCCCATTCCGTGTGTGAAATATCGTCAAGCTCGCCAACGTGGCGGGCCa                                                                              <  2:1070397/93‑1 (MQ=255)
        cGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACg                                                                            <  2:5330002/93‑1 (MQ=255)
           aTTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCa                                                                         <  2:3555198/93‑1 (MQ=255)
            ttttCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCAt                                                                        <  1:8576517/93‑1 (MQ=255)
            ttttCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCAt                                                                        >  1:966938/1‑93 (MQ=255)
            ttttCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCAt                                                                        >  2:7810723/1‑93 (MQ=255)
              ttCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATag                                                                      <  1:4379240/93‑2 (MQ=255)
              ttCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATaa                                                                      <  1:7707504/93‑1 (MQ=255)
              ttCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATaa                                                                      <  2:1260680/93‑1 (MQ=255)
              ttCCCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATaa                                                                      <  1:1105427/93‑1 (MQ=255)
                 tcGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGt                                                                   <  1:6262940/93‑1 (MQ=255)
                     cgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAc                                                               >  1:7302961/1‑93 (MQ=255)
                     cgcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAc                                                               <  1:2504605/93‑1 (MQ=255)
                      gcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAcc                                                              >  1:1285513/1‑93 (MQ=255)
                      gcTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCAcc                                                              <  1:6900584/93‑1 (MQ=255)
                       cTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCg                                                             >  1:2176491/1‑93 (MQ=255)
                        ttGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCgg                                                            <  1:8425005/93‑1 (MQ=255)
                          gcaggcagCCAGCGGCGAATGATGCCCCATACCGTGTGTGCAATATCGTCCAGCTCGCCAACGTGGCGCGGCACGGCATTACGTCCACCggga                                                          >  2:5514933/1‑91 (MQ=255)
                           caggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTaa                                                         <  2:6965518/93‑1 (MQ=255)
                           caggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTaa                                                         <  1:2713983/93‑1 (MQ=255)
                            aggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAg                                                        >  2:8699250/1‑93 (MQ=255)
                            aggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAg                                                        <  1:587582/93‑1 (MQ=255)
                            aggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAg                                                        >  1:7895253/1‑93 (MQ=255)
                            aggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAg                                                        >  1:9760056/1‑93 (MQ=255)
                             ggcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGCGGCGGGGCCCGGCATAACGTCCACCGGTAAGt                                                       >  2:2053780/1‑93 (MQ=255)
                              gcagCCAGTGGCGAATGATGCCCCCTACCGTGTGTGAAATATCGTCCAGCTCGCCACCGTGGCGGGCCACGGCATAACGTCCACCGGTAAGtt                                                      <  1:808939/93‑1 (MQ=255)
                              gcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGtt                                                      <  2:9210451/93‑1 (MQ=255)
                              gcagCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGtt                                                      <  2:8354514/93‑1 (MQ=255)
                                agCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCa                                                    >  2:4786597/1‑93 (MQ=255)
                                agCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCa                                                    <  2:9357286/93‑1 (MQ=255)
                                 gCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCAc                                                   <  1:9557936/93‑1 (MQ=255)
                                 gCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCAc                                                   >  1:6731607/1‑93 (MQ=255)
                                 gCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCAc                                                   >  2:9635697/1‑93 (MQ=255)
                                  ccAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCAcc                                                  <  1:6157083/93‑1 (MQ=255)
                                  ccAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCAcc                                                  >  1:9189354/1‑93 (MQ=255)
                                   cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCa                                                 <  1:10023200/93‑1 (MQ=255)
                                   cAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCa                                                 <  1:5163948/93‑1 (MQ=255)
                                    aGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAt                                                <  2:8393234/93‑1 (MQ=255)
                                     gTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCTACGGCATAACGTCCACCGGTAAGTTCACCAtt                                               >  2:6485015/1‑93 (MQ=255)
                                     gTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCAtt                                               >  1:7159820/1‑93 (MQ=255)
                                      tGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTg                                              <  1:987241/93‑1 (MQ=255)
                                      tGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTg                                              <  1:2525664/93‑1 (MQ=255)
                                        gCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCt                                            >  2:4442759/1‑93 (MQ=255)
                                         cGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTc                                           >  1:6444896/1‑93 (MQ=255)
                                         cGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTc                                           <  1:5486140/93‑1 (MQ=255)
                                          gAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTca                                          <  2:9841642/93‑1 (MQ=255)
                                           aaTGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCAAGGCATAACGTCCAACGGTAAGTTCACCATTGCTcca                                         >  2:7971972/1‑91 (MQ=255)
                                              gatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTTCTCACAcc                                      >  2:8995818/1‑93 (MQ=255)
                                              gatgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAATTCACCATTGCTCACAcc                                      >  1:8233705/1‑93 (MQ=255)
                                               atgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCa                                     <  1:5075493/93‑1 (MQ=255)
                                                tgCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCAACCGGAAAGTTCACCCTTGCTCACACCaa                                    >  1:3321899/1‑92 (MQ=255)
                                                 gCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATa                                   <  1:10016997/93‑1 (MQ=255)
                                                  ccccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATaa                                  <  1:8664971/93‑1 (MQ=255)
                                                  ccccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATaa                                  >  1:3836174/1‑93 (MQ=255)
                                                    ccATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACg                                <  2:1782954/93‑1 (MQ=255)
                                                     cATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGa                               <  2:9459726/93‑1 (MQ=255)
                                                     cATACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGa                               >  2:4248940/1‑93 (MQ=255)
                                                      aTACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAt                              >  1:994688/1‑93 (MQ=255)
                                                      aTACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAt                              <  2:615106/93‑1 (MQ=255)
                                                      aTACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAt                              >  1:4165853/1‑93 (MQ=255)
                                                       tACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAtt                             >  2:2872737/1‑93 (MQ=255)
                                                       tACCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGAtt                             >  1:5339045/1‑93 (MQ=255)
                                                        aCCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTa                            <  2:1444183/93‑1 (MQ=255)
                                                        aCCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTa                            <  1:9419094/93‑1 (MQ=255)
                                                        aCCGTGTGTGAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTa                            >  1:3176084/1‑93 (MQ=255)
                                                          cGTGTGTGAAATATCGTCCAGCTCGCCAGCGTGGCGGGCTACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATc                          <  2:4087188/93‑1 (MQ=255)
                                                            tgtgttAAATATGGTCCAGCTCGCAAACGTGGCGGGCCACGGCATAACGGCCACCGCTAAGTTCGCCATTGCTCAGACCATAACGATTATcag                        >  2:7400962/1‑91 (MQ=255)
                                                               gtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa                     >  1:4633339/1‑93 (MQ=255)
                                                               gtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa                     <  1:1591364/93‑1 (MQ=255)
                                                               gtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa                     >  1:6389265/1‑93 (MQ=255)
                                                               gtgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGaa                     <  1:3649006/93‑1 (MQ=255)
                                                                tgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAAt                    >  2:3917901/1‑93 (MQ=255)
                                                                tgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAAt                    >  2:2393064/1‑93 (MQ=255)
                                                                tgAAATATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAAt                    <  2:7269449/93‑1 (MQ=255)
                                                                  aaaTATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCg                  <  2:4657954/93‑1 (MQ=255)
                                                                   aaTATCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCgg                 >  2:8038743/1‑93 (MQ=255)
                                                                    atatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGt                <  1:3819739/93‑1 (MQ=255)
                                                                    atatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGt                <  2:6634728/93‑1 (MQ=255)
                                                                    atatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGt                >  1:3492035/1‑93 (MQ=255)
                                                                     tatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGtt               <  2:2665697/93‑1 (MQ=255)
                                                                     tatCGTCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGtt               <  1:4008713/93‑1 (MQ=255)
                                                                          tCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTa          >  2:6845742/1‑93 (MQ=255)
                                                                          tCCAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTa          >  1:887652/1‑93 (MQ=255)
                                                                           ccAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTaa         <  2:5065298/93‑1 (MQ=255)
                                                                           ccAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTaa         >  1:1735866/1‑93 (MQ=255)
                                                                           ccAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTaa         <  2:5818551/93‑1 (MQ=255)
                                                                            cAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAAc        <  2:7123621/93‑1 (MQ=255)
                                                                            cAGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGGAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAAc        >  1:1987359/1‑93 (MQ=255)
                                                                             aGCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACg       >  1:7455407/1‑93 (MQ=255)
                                                                              gCTCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGc      >  2:231161/1‑93 (MQ=255)
                                                                                tCGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGCTg    <  1:8522849/93‑1 (MQ=255)
                                                                                 cGCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGCTGc   <  1:9501468/93‑1 (MQ=255)
                                                                                  gCCAACGTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGCTGcc  >  2:3693546/1‑93 (MQ=255)
                                                                                       ||                                                                                      
GTGCTTTACGCATTTTCTCGCCGCTTGCAGGCAGCCAGTGGCGAATGATGCCCCATACCGTGTGTGAAATATCGTCCAGCTCGCCAAGCTGGCGGGCCACGGCATAACGTCCACCGGTAAGTTCACCATTGCTCACACCATAACGATTATCGGGAATCGGTTCGCTAACGCTGCC  >  NC_000913/3166848‑3167022

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