Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,953,437 A→G 6.8% intergenic (+24/‑41) ilvE → / → ilvD branched‑chain amino acid aminotransferase/dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,953,4370AG 6.8% 230.7 / 6.0 73intergenic (+24/‑41)ilvE/ilvDbranched‑chain amino acid aminotransferase/dihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base G (2/3);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

TATTCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGC  >  NC_000913/3953293‑3953580
                                                                                                                                                |                                                                                                                                               
tATTCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGt                                                                                                                                             <  4:282418/149‑1 (MQ=255)
 aTTCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTc                                                                                                                                            <  3:236959/149‑1 (MQ=255)
 aTTCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTc                                                                                                                                            >  1:450909/1‑149 (MQ=255)
  ttCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTcc                                                                                                                                           >  3:563725/1‑149 (MQ=255)
   tCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTccc                                                                                                                                          >  1:72131/1‑149 (MQ=255)
      ggTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCAtt                                                                                                                                       <  4:532352/149‑1 (MQ=255)
       gTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAttt                                                                                                                                      <  4:99795/149‑1 (MQ=255)
           gCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACga                                                                                                                                  <  4:42730/149‑1 (MQ=255)
           gCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACga                                                                                                                                  <  4:92538/149‑1 (MQ=255)
               aGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGTGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGagaca                                                                                                                              <  3:405249/149‑1 (MQ=255)
               aGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGagaca                                                                                                                              <  3:501357/149‑1 (MQ=255)
                         ggCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGa                                                                                                                    <  2:450909/149‑1 (MQ=255)
                          gCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAg                                                                                                                   >  1:415113/1‑149 (MQ=255)
                                  aCCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaaataaa                                                                                                           <  4:322365/149‑1 (MQ=255)
                                   ccAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                                                                      >  2:468838/1‑137 (MQ=255)
                                   ccAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                                                                      <  1:468838/137‑1 (MQ=255)
                                   ccAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAg                                                                                                          <  1:270847/149‑1 (MQ=255)
                                    cAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTg                                                                                                                       <  4:146303/135‑1 (MQ=255)
                                    cAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTg                                                                                                                       >  3:146303/1‑135 (MQ=255)
                                     aaaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTa                                                                                                        >  3:55595/1‑149 (MQ=255)
                                     aaaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTa                                                                                                        >  1:291240/1‑149 (MQ=255)
                                       aCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATg                                                                                                      <  2:130752/149‑1 (MQ=255)
                                           aTTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAtttttt                                                                                                                                   <  4:214532/116‑4 (MQ=39)
                                           aTTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAtttttt                                                                                                                                   >  3:214532/1‑113 (MQ=39)
                                           aTTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAg                                                                                                          <  2:155931/141‑1 (MQ=255)
                                           aTTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAg                                                                                                          >  1:155931/1‑141 (MQ=255)
                                            ttCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTaa                                                                                                 >  3:581419/1‑149 (MQ=255)
                                                gCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                                                             >  4:427042/1‑149 (MQ=255)
                                                      cttcttCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTcc                                                                                       <  1:406863/149‑1 (MQ=255)
                                                       ttcttcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCg                                                                                      >  1:525954/1‑149 (MQ=255)
                                                             ggCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacca                                                                                >  3:24926/1‑149 (MQ=255)
                                                                 tctTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCcaccaccac                                                                            >  3:575457/1‑149 (MQ=255)
                                                                  ctTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacc                                                                                 >  1:38902/1‑143 (MQ=255)
                                                                  ctTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGccacc                                                                                 <  2:38902/143‑1 (MQ=255)
                                                                    tCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCa                                                                         <  3:427042/149‑1 (MQ=255)
                                                                     cACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCAt                                                                        <  2:158714/149‑1 (MQ=255)
                                                                      aCTGGAGAAACCGAAGATAAATGGGGCTGGTTAGATCAAATTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATg                                                                       <  3:216885/149‑1 (MQ=255)
                                                                         ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGCCGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTc                                                                    <  1:408849/149‑1 (MQ=255)
                                                                         ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGTGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTc                                                                    <  4:120081/149‑1 (MQ=255)
                                                                              aaCCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTaa                                                                >  4:579578/1‑148 (MQ=255)
                                                                              aaCCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTaa                                                                <  3:579578/148‑1 (MQ=255)
                                                                               aCCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAata                                                              <  2:197706/149‑1 (MQ=255)
                                                                                 cGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATg                                                            >  4:562558/1‑149 (MQ=255)
                                                                                  gAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATgg                                                           >  2:383980/1‑149 (MQ=255)
                                                                                  gAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATgg                                                           >  3:99816/1‑149 (MQ=255)
                                                                                     gATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCgg                                                        <  2:415113/149‑1 (MQ=255)
                                                                                     gATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCgg                                                        >  2:28577/1‑149 (MQ=255)
                                                                                     gATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCgg                                                        <  2:72131/149‑1 (MQ=255)
                                                                                                    ttAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCtgt                                         <  1:174556/149‑1 (MQ=255)
                                                                                                      aGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTgg                                       <  3:581926/149‑1 (MQ=255)
                                                                                                       gATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGgc                                      <  1:53772/149‑1 (MQ=255)
                                                                                                         tCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGgcgc                                    <  3:352270/149‑1 (MQ=255)
                                                                                                         tCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGgcgc                                    >  3:91310/1‑149 (MQ=255)
                                                                                                         tCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGgcgc                                    >  4:229072/1‑149 (MQ=255)
                                                                                                            aGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATgg                                                                      >  2:210022/1‑112 (MQ=255)
                                                                                                            aGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATgg                                                                      <  1:210022/112‑1 (MQ=255)
                                                                                                             gTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCt                                           >  3:259317/1‑138 (MQ=255)
                                                                                                             gTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCt                                           <  4:259317/138‑1 (MQ=255)
                                                                                                             gTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGGGCTCGTGCGCTGTGGCGCGCCa                                >  2:503020/1‑149 (MQ=255)
                                                                                                                   cAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGaa                          <  2:525954/149‑1 (MQ=255)
                                                                                                                    aataaataCAAAAAATGGGCCGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGCGTGCTCGTTGGCTTTGGCGGGCCACCGGata                         >  3:243168/1‑147 (MQ=255)
                                                                                                                    aataaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAAt                         >  2:112992/1‑149 (MQ=255)
                                                                                                                       aaataCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGAc                      <  4:357358/149‑1 (MQ=255)
                                                                                                                            tcgtaaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATgg                                                                      >  2:203913/5‑96 (MQ=37)
                                                                                                                            tcgtaaaTGGGACGGTGCGTGCCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATgg                                                                      <  1:203913/92‑1 (MQ=37)
                                                                                                                             aaaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGc                <  3:283768/149‑1 (MQ=255)
                                                                                                                                     ggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCg        >  1:554496/1‑149 (MQ=255)
                                                                                                                                      gACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGaa                          <  2:99746/130‑1 (MQ=255)
                                                                                                                                      gACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGaa                          >  1:99746/1‑130 (MQ=255)
                                                                                                                                       aCGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGt      >  2:417140/1‑149 (MQ=255)
                                                                                                                                         ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTaa    <  4:575457/149‑1 (MQ=255)
                                                                                                                                         ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTaa    <  1:28577/149‑1 (MQ=255)
                                                                                                                                         ggCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCATGGGCTGTGGGGCGCCACCGGGAAGACCGACGCCGATTTCCGTaa    >  1:104951/1‑149 (MQ=255)
                                                                                                                                           cacgcacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGc  >  1:140979/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TATTCAGGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGC  >  NC_000913/3953293‑3953580

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