Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,959,480 A→T 18.9% intergenic (+35/+52) ilvC → / ← ppiC ketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,959,4800AT18.9% 225.3 / 47.9 90intergenic (+35/+52)ilvC/ppiCketol‑acid reductoisomerase, NAD(P)‑binding/peptidyl‑prolyl cis‑trans isomerase C (rotamase C)
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base T (9/8);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 7.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

AACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACT  >  NC_000913/3959332‑3959622
                                                                                                                                                    |                                                                                                                                              
aaCGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCca                                                                                                                                                >  4:682553/1‑149 (MQ=255)
  cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                              <  2:21468/149‑1 (MQ=255)
  cGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCcaca                                                                                                                                              <  3:682553/149‑1 (MQ=255)
      cAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                          <  4:1087297/149‑1 (MQ=255)
      cAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                          <  2:601940/149‑1 (MQ=255)
       aaCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggag                                                                                                                                         >  4:27026/1‑149 (MQ=255)
        aCTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                        <  2:82642/149‑1 (MQ=255)
         cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagag                                                                                                                                       <  3:770531/149‑1 (MQ=255)
         cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagag                                                                                                                                       <  4:241689/149‑1 (MQ=255)
         cTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagag                                                                                                                                       <  3:959285/149‑1 (MQ=255)
          tGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                      <  1:1454371/149‑1 (MQ=255)
          tGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                      <  1:550165/149‑1 (MQ=255)
             gTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                      <  3:1401290/146‑1 (MQ=255)
             gTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                      >  4:1401290/1‑146 (MQ=255)
             gTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                      >  4:1464393/1‑146 (MQ=255)
                   tgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagag                                                                                                                                       <  3:1400284/139‑1 (MQ=255)
                   tgAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagag                                                                                                                                       >  4:1400284/1‑139 (MQ=255)
                    gAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                                             <  4:340193/132‑1 (MQ=255)
                    gAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACag                                                                                                                                             >  3:340193/1‑132 (MQ=255)
                        gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                      >  2:393312/1‑135 (MQ=255)
                        gAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                      <  1:393312/135‑1 (MQ=255)
                           gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTTCGGGTTAAGTGCGCGCTGATGTCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAaca                                                                                                                             <  1:522949/141‑1 (MQ=255)
                           gCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAaca                                                                                                                             >  2:522949/1‑141 (MQ=255)
                            cGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGcc                                                                                                                    <  3:27026/149‑1 (MQ=255)
                             gATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                      <  1:1015579/130‑1 (MQ=255)
                             gATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                      >  2:1015579/1‑130 (MQ=255)
                                  gCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTcc                                                                                                              >  4:227314/1‑149 (MQ=255)
                                      ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                                     >  3:1204977/1‑122 (MQ=255)
                                      ccATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                                     <  4:1204977/122‑1 (MQ=255)
                                        aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCt                                                                                                                   >  3:224784/1‑138 (MQ=255)
                                        aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCt                                                                                                                   <  4:224784/138‑1 (MQ=255)
                                        aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGga                                                                                                          <  2:398052/147‑1 (MQ=255)
                                        aTGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCAGga                                                                                                          >  1:398052/1‑147 (MQ=255)
                                           cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                        <  2:149803/130‑1 (MQ=255)
                                           cGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                        >  1:149803/1‑130 (MQ=255)
                                            gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                      <  2:212633/115‑1 (MQ=255)
                                            gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                      >  1:212633/1‑115 (MQ=255)
                                            gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggaga                                                                                                                                        <  3:522069/113‑1 (MQ=255)
                                            gATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggaga                                                                                                                                        >  4:522069/1‑113 (MQ=255)
                                             aTTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                      >  4:583900/1‑114 (MQ=255)
                                              ttGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                  <  3:227314/149‑1 (MQ=255)
                                                     ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                       <  4:119928/121‑1 (MQ=255)
                                                     ggtaggtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAg                                                                                                                       >  3:119928/1‑121 (MQ=255)
                                                          gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagggt                                                                                                                                    >  3:1419891/1‑102 (MQ=255)
                                                          gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagggt                                                                                                                                    <  4:1419891/103‑2 (MQ=38)
                                                          gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                          <  2:1244116/97‑1 (MQ=255)
                                                          gtaAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                          >  1:1244116/1‑97 (MQ=255)
                                                             aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                          >  3:549343/1‑94 (MQ=255)
                                                             aGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                          <  4:549343/94‑1 (MQ=255)
                                                                           tataTGACATATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                  <  4:109226/120‑1 (MQ=255)
                                                                           tataTGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                  >  3:109226/1‑120 (MQ=255)
                                                                               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                        >  2:923828/1‑94 (MQ=255)
                                                                               tGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCaa                                                                                                                        <  1:923828/94‑1 (MQ=255)
                                                                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                      <  3:1427092/78‑1 (MQ=255)
                                                                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagagg                                                                                                                                      >  4:1427092/1‑78 (MQ=255)
                                                                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                          >  4:394491/1‑74 (MQ=255)
                                                                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACaggga                                                                                                                                          <  3:394491/74‑1 (MQ=255)
                                                                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                    >  3:122851/1‑112 (MQ=255)
                                                                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                    >  1:383116/1‑112 (MQ=255)
                                                                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                    <  2:383116/112‑1 (MQ=255)
                                                                                 aCAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCACTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                    <  4:122851/112‑1 (MQ=255)
                                                                                   aGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagagg                                                                                                                                      >  3:501232/1‑76 (MQ=37)
                                                                                       atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                  <  1:1259832/108‑1 (MQ=255)
                                                                                       atGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                  >  2:1259832/1‑108 (MQ=255)
                                                                                        tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                      >  1:190042/1‑103 (MQ=255)
                                                                                        tGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAg                                                                                                      <  2:190042/103‑1 (MQ=255)
                                                                                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                    >  3:617020/1‑104 (MQ=255)
                                                                                         gAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGc                                                                                                    <  4:617020/104‑1 (MQ=255)
                                                                                          aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                                     <  1:1225198/70‑1 (MQ=255)
                                                                                          aaaCGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggagaggg                                                                                                                                     >  2:1225198/1‑70 (MQ=255)
                                                                                                aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                        <  3:465979/61‑1 (MQ=255)
                                                                                                aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                        >  4:465979/1‑61 (MQ=255)
                                                                                                aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                        >  4:676892/1‑61 (MQ=255)
                                                                                                aTTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACagggaga                                                                                                                                        <  3:676892/61‑1 (MQ=255)
                                                                                                    cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGcc                                                                                                   >  4:1296423/1‑94 (MQ=255)
                                                                                                    cTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACGCAAGGCCTTCTCCAGGAGAAGGcc                                                                                                   <  3:1296423/94‑1 (MQ=255)
                                                                                                     tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                              >  1:590090/1‑130 (MQ=255)
                                                                                                     tGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGa                                                              <  2:590090/130‑1 (MQ=255)
                                                                                                             gTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGGAAACACTCAAGGc                                                                                                                     <  2:1256611/67‑1 (MQ=255)
                                                                                                             gTTAAGTGCACGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGGAAACACTCAAGGc                                                                                                                     >  1:1256611/1‑67 (MQ=255)
                                                                                                                  gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagaggg                                                                                                                                     >  1:857350/1‑46 (MQ=11)
                                                                                                                  gTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCTGtgggagaggg                                                                                                                                     <  2:857350/46‑1 (MQ=11)
                                                                                                                        gcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                       >  2:731122/1‑70 (MQ=255)
                                                                                                                        gcTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGaa                                                                                                       <  1:731122/70‑1 (MQ=255)
                                                                                                                                      cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                           <  2:150930/132‑1 (MQ=255)
                                                                                                                                      cccGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTc                           >  1:150930/1‑132 (MQ=255)
                                                                                                                                          tcccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTata                                                                               <  1:1212121/75‑3 (MQ=38)
                                                                                                                                          tcccTCTCCCTGTGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTata                                                                               >  2:1212121/2‑74 (MQ=38)
                                                                                                                                           cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAc                     >  2:837488/1‑133 (MQ=255)
                                                                                                                                           cccTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTAc                     <  1:837488/133‑1 (MQ=255)
                                                                                                                                              tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                  <  1:289341/53‑1 (MQ=255)
                                                                                                                                              tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCt                                                                                                  >  2:289341/1‑53 (MQ=255)
                                                                                                                                              tctcCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACt  >  1:485822/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AACGGGCAACTGCGTGATGTGAACGAAGCGATTCGCAGCCATGCGATTGAGCAGGTAGGTAAGAAACTGCGCGGCTATATGACAGATATGAAACGTATTGCTGTTGCGGGTTAAGTGCGCGCTGATGCCCTCACCCCGACCCTCTCCCACAGGGAGAGGGAGAAAACACTCAAGGCCTTCTCCTGGAGAAGGCCTTGCTATTAGTTGCGGTACAGCACCTTAATGATGTGATATCCGAACTGGGTGTGCAGCGGGCCGGTCGGCTCCAGTACCGGACAAGAGAAAACCACT  >  NC_000913/3959332‑3959622

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