Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,476 (T)8→9 intergenic (+56/‑31) ilvL → / → ilvX ilvG operon leader peptide/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,4681.T94.0% 211.4 / 5.2 67intergenic (+48/‑39)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (30/33);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCAC  >  NC_000913/3950333‑3950599
                                                                                                                                        |                                                                                                                                    
tCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑tt‑t                                                                                                                                  >  1:728278/1‑139 (MQ=255)
       gTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑tt‑tttttttg                                                                                                                           <  1:439051/139‑3 (MQ=255)
          aatAGCCTGCTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑tt‑tttttttgacc                                                                                                                        <  2:728278/137‑6 (MQ=255)
          aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑tt‑tttttttgacc                                                                                                                        >  2:134921/1‑134 (MQ=255)
              gCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTaa                                                                                                                     <  1:617617/139‑1 (MQ=255)
                        gATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAAcc                                                                                                           <  1:37009/139‑1 (MQ=255)
                         aTTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCg                                                                                                          >  1:334183/1‑139 (MQ=255)
                          ttAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCga                                                                                                         <  1:58407/139‑1 (MQ=255)
                              cGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCgaggag                                                                                                     >  1:926917/1‑139 (MQ=255)
                                 ggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                                                                 <  1:846371/139‑1 (MQ=255)
                                            ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATaa                                                                                       <  2:1143965/139‑1 (MQ=255)
                                               tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                                                    >  2:54744/1‑139 (MQ=255)
                                               tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑t                                                                                                                                     >  2:566103/1‑90 (MQ=255)
                                               tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑t                                                                                                                                     <  1:566103/90‑1 (MQ=255)
                                                  cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcac                                                                                 >  2:558182/1‑139 (MQ=255)
                                                    ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCga                                                                                                         <  2:701231/113‑1 (MQ=255)
                                                    ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCga                                                                                                         >  1:701231/1‑113 (MQ=255)
                                                    ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaa                                                                               >  1:405270/1‑139 (MQ=255)
                                                    ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaa                                                                               >  1:845739/1‑139 (MQ=255)
                                                    ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaa                                                                               >  2:366844/1‑139 (MQ=255)
                                                    ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaa                                                                               >  2:751552/1‑139 (MQ=255)
                                                     cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCgagga                                                                                                      >  1:408884/1‑115 (MQ=255)
                                                     cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCgagga                                                                                                      <  2:408884/115‑1 (MQ=255)
                                                      gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                                                             >  2:183473/1‑139 (MQ=255)
                                                       tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTa                                                                                                                      <  1:864387/97‑1 (MQ=255)
                                                       tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTa                                                                                                                      >  2:864387/1‑97 (MQ=255)
                                                       tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCgaggag                                                                                                     <  2:1040421/114‑1 (MQ=255)
                                                       tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCgaggag                                                                                                     >  1:1040421/1‑114 (MQ=255)
                                                       tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                                                   >  1:482387/1‑115 (MQ=255)
                                                       tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                                                   <  2:482387/115‑1 (MQ=255)
                                                         cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                                                                  <  1:663006/115‑1 (MQ=255)
                                                         cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                                                                  >  2:663006/1‑115 (MQ=255)
                                                                 cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                                                                  >  1:984103/1‑123 (MQ=255)
                                                                 cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                                                                  <  2:984103/123‑1 (MQ=255)
                                                                 cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                                                        <  2:736395/133‑1 (MQ=255)
                                                                 cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtctc                                                                  >  2:726850/1‑139 (MQ=255)
                                                                 cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtctc                                                                  >  1:450449/1‑139 (MQ=255)
                                                                 cACTTCGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTTTg                                                                        >  1:736395/1‑133 (MQ=255)
                                                                  aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCGGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                                                                 <  2:450449/139‑1 (MQ=255)
                                                                      ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAt                                                             >  1:301958/1‑139 (MQ=255)
                                                                               aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgg                                                    <  1:134921/139‑1 (MQ=255)
                                                                               aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgg                                                    <  1:726850/139‑1 (MQ=255)
                                                                                aGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAATAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg                                                   <  2:845739/139‑1 (MQ=255)
                                                                                 ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg                                                  <  1:366844/139‑1 (MQ=255)
                                                                                          aTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAactaacta                                         >  1:1203182/1‑139 (MQ=255)
                                                                                             aaGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac                                               >  1:1030618/1‑130 (MQ=255)
                                                                                             aaGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac                                               <  2:1030618/130‑1 (MQ=255)
                                                                                             aaGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaa                                                >  2:926895/1‑129 (MQ=255)
                                                                                             aaGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaa                                                <  1:926895/129‑1 (MQ=255)
                                                                                             aaGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGTACTAACTATGa                                      <  1:853577/139‑1 (MQ=255)
                                                                                              aGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa                                     <  1:558182/139‑1 (MQ=255)
                                                                                                    aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                                                   >  2:873692/1‑71 (MQ=255)
                                                                                                    aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                                                   <  1:873692/71‑1 (MQ=255)
                                                                                                       gAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGcaca                            <  2:1203182/139‑1 (MQ=255)
                                                                                                       gAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGcaca                            <  1:751552/139‑1 (MQ=255)
                                                                                                          cTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCtt                                                                                                                       <  1:106115/45‑1 (MQ=255)
                                                                                                          cTAAGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCtt                                                                                                                       >  2:106115/1‑45 (MQ=255)
                                                                                                            aaGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAgtggg                       <  2:405270/139‑1 (MQ=255)
                                                                                                            aaGACCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAgtggg                       <  1:328571/139‑1 (MQ=255)
                                                                                                               aCCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                                                                   >  1:706796/1‑92 (MQ=255)
                                                                                                               aCCCCCGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                                                                   <  2:706796/92‑1 (MQ=255)
                                                                                                                   ccGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAGAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGca                              >  1:745114/1‑125 (MQ=255)
                                                                                                                   ccGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                                                          >  1:776159/1‑97 (MQ=255)
                                                                                                                   ccGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaa                                                >  1:1101618/1‑107 (MQ=255)
                                                                                                                   ccGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaa                                                <  2:1101618/107‑1 (MQ=255)
                                                                                                                   ccGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGca                              <  2:745114/125‑1 (MQ=255)
                                                                                                                   ccGCACCGAAAGGTCCGGGGG‑TTGTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAGAGCACAAAATTCTGTTTCTCAAGATTCa                                                          <  2:776159/97‑1 (MQ=255)
                                                                                                                    cGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCgaggaac                                                                                                    >  2:1060924/1‑52 (MQ=255)
                                                                                                                    cGCACCGAAAGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGc                                                                                                    <  1:1060924/54‑1 (MQ=255)
                                                                                                                           aaaGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg                                                   >  2:674975/1‑96 (MQ=255)
                                                                                                                           aaaGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg                                                   <  1:674975/96‑1 (MQ=255)
                                                                                                                            aaGGTCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTAGGAATGGCGCACAGTGGGTGGTACATGCGTTGCg       <  2:334183/139‑1 (MQ=255)
                                                                                                                                tCCGGGGGTTT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGca   <  2:301958/139‑1 (MQ=255)
                                                                                                                                  gggggt‑TT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGcac  >  1:18113/2‑137 (MQ=255)
                                                                                                                                  gggggt‑TT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGcac  <  2:18113/136‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                    
TCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGG‑TT‑TTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAATGGCGCACAGTGGGTGGTACATGCGTTGCGGGCAC  >  NC_000913/3950333‑3950599

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: ‑