Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,293,960 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate‑aspartate symport protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,293,9591.C99.1% 217.6 / ‑3.6 113intergenic (+586/+56)gltP/yjcOglutamate‑aspartate symport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (53/59);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01
*NC_0009134,293,9592.G99.1% 216.4 / ‑3.6 113intergenic (+586/+56)gltP/yjcOglutamate‑aspartate symport protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (53/59);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

Only 100 of 116 total aligned reads displayed.

GCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCT  >  NC_000913/4293929‑4294042
                               ||                                                                                   
gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                         <  1:6554344/93‑1 (MQ=1)
gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                         <  2:6167578/93‑1 (MQ=1)
gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                         <  1:4572914/93‑1 (MQ=1)
gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                         >  2:6547423/1‑93 (MQ=1)
gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                         <  2:4601932/93‑1 (MQ=1)
gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                         <  1:6181071/93‑1 (MQ=1)
gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                         >  1:1975887/1‑93 (MQ=1)
gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                         <  1:163956/93‑1 (MQ=1)
gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                         >  1:6612484/1‑93 (MQ=1)
 ggTGATTTTGTAGGTCGGAGAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                        <  2:6796206/92‑1 (MQ=11)
 cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                        >  2:8445983/1‑93 (MQ=11)
 cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                        >  1:1181999/1‑93 (MQ=11)
  gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                       >  1:5509626/1‑93 (MQ=11)
  gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                       >  2:1558026/1‑93 (MQ=11)
    gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGTCCTCAATGCCTGATGCGAGGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                     >  1:344660/1‑93 (MQ=11)
    gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                     <  1:1906157/93‑1 (MQ=11)
    gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                     <  1:7240079/93‑1 (MQ=11)
    gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                     >  1:503883/1‑93 (MQ=11)
    gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                     >  1:7915955/1‑93 (MQ=11)
    gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                     >  1:344662/1‑93 (MQ=11)
    gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGATTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                     <  1:1318474/93‑1 (MQ=11)
     aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                    >  2:5441234/1‑93 (MQ=11)
     aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                    <  2:563910/93‑1 (MQ=11)
     aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                    >  1:5949585/1‑93 (MQ=11)
     aTTTTGTAGGTCGGATAAGGCGTCTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                    <  2:4835841/93‑1 (MQ=11)
      ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                   >  2:3665076/1‑93 (MQ=11)
      ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                   <  1:2104548/93‑1 (MQ=11)
        ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                 <  1:7477197/93‑1 (MQ=14)
        ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                 <  1:2080648/93‑1 (MQ=14)
        ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                 >  1:2977150/1‑93 (MQ=14)
         tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                <  1:8096922/93‑1 (MQ=14)
         tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                >  1:5754083/1‑93 (MQ=14)
          gTAGGTCGGATAAGGCGTTTACGCCGCATCCGGCATCAATGTCTGATGCGACTCTTGTCGCGTCTGATCAGGCCTATCTTAACCGCTGGTTaa               >  2:7789333/1‑93 (MQ=14)
          gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCATTGCCTGATGCGACGCTTGCCGCGTCTTCTCAGGCCTATCTTAACCGTTGGTTaa               <  2:7808996/93‑1 (MQ=14)
          gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTat               >  1:3793036/1‑92 (MQ=14)
          gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa               <  1:3776464/93‑1 (MQ=14)
          gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa               <  1:6269467/93‑1 (MQ=14)
          gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa               <  1:7462381/93‑1 (MQ=14)
          gTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa               <  2:2502474/93‑1 (MQ=14)
          gTAGGTCGGATAAGGCGTTGACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa               <  2:600062/93‑1 (MQ=14)
          gTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCAATGCCATATGCGACGCTTGCCGCGTTTTATCAGGCCTATTTTAACCGTTGGCTaa               >  2:6674967/1‑93 (MQ=11)
           tAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt              <  2:5822256/93‑1 (MQ=14)
           tAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt              <  2:6319465/93‑1 (MQ=14)
            tggTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt             <  1:6582560/92‑1 (MQ=14)
            aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt             <  1:4572008/93‑1 (MQ=14)
            aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt             <  2:17847/93‑1 (MQ=14)
             ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt            >  2:5607949/1‑93 (MQ=14)
             ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt            >  1:3305290/1‑93 (MQ=14)
              gTCGGATAAGGCGTTTACGCCGCATCCGCCGTCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt           <  2:7785658/93‑1 (MQ=17)
              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt           >  2:409967/1‑93 (MQ=17)
              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt           <  2:2985542/93‑1 (MQ=17)
              gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt           >  1:151163/1‑93 (MQ=17)
               tCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt          <  1:919332/93‑1 (MQ=17)
               tCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt          <  1:2471758/93‑1 (MQ=17)
               gcGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt          <  2:6554985/92‑1 (MQ=17)
                cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc         <  1:3837506/93‑1 (MQ=17)
                cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc         <  1:3201739/93‑1 (MQ=17)
                cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc         <  2:1544507/93‑1 (MQ=17)
                cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc         <  1:4173204/93‑1 (MQ=17)
                cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc         >  2:267578/1‑93 (MQ=17)
                cGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc         <  2:5339956/93‑1 (MQ=17)
                 ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGGCGCTTGCCGCGTCTTATACACCCTATCTTAACCGTTGTTTAATTTTTCa        >  1:8654432/1‑93 (MQ=14)
                 ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa        <  1:3313344/93‑1 (MQ=17)
                 ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa        >  2:1723391/1‑93 (MQ=17)
                 ggATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa        >  1:2168603/1‑93 (MQ=17)
                  gATAAGGCGTTTACGCCGCATCCGACATCAATGCTTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       >  1:7350775/1‑93 (MQ=17)
                  gATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       >  2:7193032/1‑93 (MQ=17)
                  gATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       >  2:6425306/1‑93 (MQ=17)
                  gATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       >  2:7193031/1‑93 (MQ=17)
                  gATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       >  1:6754576/1‑93 (MQ=17)
                  gATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       <  2:3709526/93‑1 (MQ=17)
                  gATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       >  2:4432540/1‑93 (MQ=17)
                  gATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       <  2:4575317/93‑1 (MQ=17)
                   aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa      >  2:554923/1‑93 (MQ=17)
                   aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa      <  1:3740829/93‑1 (MQ=17)
                   aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa      <  1:3740830/93‑1 (MQ=17)
                   aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa      <  1:8379588/93‑1 (MQ=17)
                   aTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa      <  1:5603009/93‑1 (MQ=17)
                    tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc     >  2:909921/1‑93 (MQ=17)
                    tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc     <  2:2556698/93‑1 (MQ=17)
                    tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc     >  2:2131880/1‑93 (MQ=18)
                    tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc     >  2:1239016/1‑93 (MQ=17)
                    tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc     <  1:6567866/93‑1 (MQ=17)
                    tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc     >  1:4537978/1‑93 (MQ=18)
                    tAAGGCGTGTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc     <  1:4495298/93‑1 (MQ=17)
                     aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct    >  2:8674772/1‑93 (MQ=18)
                     aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct    >  1:2325369/1‑93 (MQ=18)
                     aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct    >  1:8198733/1‑93 (MQ=18)
                     aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct    <  1:6478939/93‑1 (MQ=18)
                     aaGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct    >  1:5281666/1‑93 (MQ=18)
                      aGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc   >  1:1624882/1‑93 (MQ=18)
                      aGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc   <  1:2833191/93‑1 (MQ=18)
                      aGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc   >  1:3756693/1‑93 (MQ=18)
                      aGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc   >  1:5094332/1‑93 (MQ=18)
                      aGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc   >  2:6903390/1‑93 (MQ=18)
                       ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAGGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct  <  2:55035/93‑1 (MQ=18)
                       ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct  <  1:7321277/93‑1 (MQ=18)
                       ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct  >  1:7687704/1‑93 (MQ=18)
                       ggCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct  >  2:2942343/1‑93 (MQ=18)
                       ggCGTTTACGCCGCATCCGACATCAATGCCTGAGGCGACGCTTGCCTCGCCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct  <  1:8485964/93‑1 (MQ=18)
                               ||                                                                                   
GCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCT  >  NC_000913/4293929‑4294042

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