Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,381 +GC intergenic (+587/+55) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,3801.C97.8% 180.8 / ‑2.3 93intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base C (37/53);  minor base T (2/0);  total (40/53)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01
*NC_0009134,296,3802.G98.9% 189.5 / ‑2.9 93intergenic (+586/+56)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (39/53);  total (40/53)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296342‑4296470
                                       ||                                                                                            
cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTaa                                 >  2:3210390/1‑100 (MQ=1)
cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTaa                                 >  1:4087262/1‑100 (MQ=1)
 tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                >  1:3914235/1‑100 (MQ=11)
 tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                <  1:3751367/100‑1 (MQ=11)
  gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                               <  1:1099152/100‑1 (MQ=11)
  gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                               >  1:493784/1‑100 (MQ=11)
  gATTTAGCGTGATTCTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCATATCAGGCCTATCTTAAcc                               <  2:4110704/100‑1 (MQ=11)
  gATATAGCGTGATTTTGTAAGTCGGATAAGGCGTTTATGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCTGGCCTATCTTAAcc                               >  2:4112915/1‑100 (MQ=11)
   aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                              <  2:981885/100‑1 (MQ=11)
   aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                              <  2:1469902/100‑1 (MQ=11)
   aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTATCTTAACCg                              <  1:2424413/100‑1 (MQ=11)
    tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                             >  1:3574597/1‑100 (MQ=11)
     ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                            <  1:3173117/100‑1 (MQ=11)
     ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                            <  2:2600564/100‑1 (MQ=11)
     ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                            <  1:2689898/100‑1 (MQ=11)
      naGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                           >  1:377756/2‑100 (MQ=11)
       aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                          <  1:39825/100‑1 (MQ=11)
       aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                          >  2:819431/1‑100 (MQ=11)
       aGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                          <  1:40358/100‑1 (MQ=11)
        gagttaTTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         <  2:969135/95‑1 (MQ=11)
        gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         >  1:1432717/1‑100 (MQ=14)
        gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         >  1:3218523/1‑100 (MQ=14)
        gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                         >  1:975993/1‑100 (MQ=14)
        gCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCCTATCAGGCCTATCTTAACCGTTGGt                         >  1:439284/1‑100 (MQ=11)
         cGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        <  2:831310/100‑1 (MQ=14)
         cGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        >  2:2199636/1‑100 (MQ=14)
         cGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        >  2:3140499/1‑100 (MQ=14)
         cGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                        >  2:1109789/1‑100 (MQ=14)
          gTGATTTTGTAGGTCGGATAAGGCGTTTATGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  1:3388294/1‑100 (MQ=14)
          gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       >  2:342461/1‑100 (MQ=14)
          gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                       <  2:2650006/100‑1 (MQ=14)
          gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGCTGGTTa                       >  1:1825230/1‑100 (MQ=11)
           tGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                      <  2:3576828/100‑1 (MQ=14)
           tGATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                      >  2:1525370/1‑100 (MQ=14)
            gATTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                     <  2:3039060/100‑1 (MQ=14)
            gATTTTGCAGGCCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                     <  2:4127147/100‑1 (MQ=14)
            gATCTTGTAAGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                     <  2:1887531/100‑1 (MQ=14)
             aTTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                    >  2:741446/1‑100 (MQ=14)
             aTTTTGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                    <  1:3678459/100‑1 (MQ=14)
              ttttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATTCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                   <  1:3565432/100‑1 (MQ=14)
               tttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  <  1:226648/100‑1 (MQ=14)
               tttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  <  2:460381/100‑1 (MQ=14)
               tttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  <  1:3721190/100‑1 (MQ=17)
               tttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  <  1:2961619/100‑1 (MQ=17)
               tttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGAGGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  <  2:460402/100‑1 (MQ=14)
               tttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                  <  2:460513/100‑1 (MQ=14)
                ttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                 <  1:360245/100‑1 (MQ=17)
                ttGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                 <  1:2467844/100‑1 (MQ=17)
                 tgtgGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  1:924354/100‑1 (MQ=17)
                 tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                >  2:4019217/1‑100 (MQ=17)
                 tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  1:1366192/100‑1 (MQ=17)
                 tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                >  2:3593779/1‑100 (MQ=17)
                 tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                >  2:1763392/1‑100 (MQ=17)
                 tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                >  2:1818508/1‑100 (MQ=17)
                 tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  1:2225624/100‑1 (MQ=17)
                 tGTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  1:116605/100‑1 (MQ=17)
                 tGTAGGTCGGATAAGGCGGTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  2:1308278/100‑1 (MQ=17)
                 tGTAGGTCGGACAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                <  2:1508904/100‑1 (MQ=17)
                  gTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa               <  2:3788900/100‑1 (MQ=17)
                  gTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa               <  1:597378/100‑1 (MQ=17)
                  gTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa               <  1:3157700/100‑1 (MQ=17)
                  gTAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa               >  1:1350026/1‑100 (MQ=17)
                   tAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATTAATGCCTGATGCGGCGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa              <  2:4429372/100‑1 (MQ=17)
                   tAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa              <  2:2439546/100‑1 (MQ=17)
                   tAGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa              <  2:1702337/100‑1 (MQ=17)
                    aGGTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa             >  2:483721/1‑100 (MQ=17)
                    aGGTCGGATAAGGCGGTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa             <  1:1248510/100‑1 (MQ=17)
                     ggTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            <  1:3940566/100‑1 (MQ=17)
                     ggTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  1:1400100/1‑100 (MQ=17)
                     ggTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  2:2951425/1‑100 (MQ=17)
                     ggTCGGATAAGGCGTTTA‑‑CACAGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc            >  1:3902324/1‑100 (MQ=18)
                      gTCGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct           >  1:284798/1‑100 (MQ=17)
                        cGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct         <  2:3542544/100‑1 (MQ=18)
                        cGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct         <  2:2435828/100‑1 (MQ=18)
                        cGGATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct         >  1:1882595/1‑100 (MQ=18)
                         ggATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        <  1:1046096/100‑1 (MQ=18)
                         ggATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        >  2:3608438/1‑100 (MQ=18)
                         ggATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt        <  2:2154420/100‑1 (MQ=18)
                          gATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  2:1945944/100‑1 (MQ=18)
                          gATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  1:2295460/100‑1 (MQ=18)
                          gATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  1:3152260/100‑1 (MQ=18)
                          gATAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc       <  2:1457908/100‑1 (MQ=18)
                           aTAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca      <  1:257521/100‑1 (MQ=18)
                            tAAGGCGTTTACGC‑C‑GCATCTGACATAAATGCCTGACGCGACGCTGGACGCGTCCTATCATGACTATCTTAGCCGTTGGTTAATTTTTCAAACTCTTcac     >  2:575296/1‑100 (MQ=18)
                            tAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac     <  2:4438171/100‑1 (MQ=18)
                            tAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac     >  2:1236465/1‑100 (MQ=18)
                            tAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac     >  2:1821755/1‑100 (MQ=18)
                            tAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac     >  2:1170802/1‑100 (MQ=18)
                            tAAGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATCCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac     <  2:4075347/100‑1 (MQ=18)
                              aGGCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac   >  2:4405390/1‑100 (MQ=18)
                               ggCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTTTTTGCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc  >  1:1539128/1‑100 (MQ=21)
                               ggCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc  >  1:1551606/1‑100 (MQ=21)
                               ggCGTTTACGC‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAcc  >  1:1798419/1‑100 (MQ=21)
                                       ||                                                                                            
CTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑C‑GCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACC  >  NC_000913/4296342‑4296470

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