Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,420 G→A *33* (TAG→TAA ilvL → ilvG operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,4200GA98.8% 264.6 / NA 83*33* (TAG→TAAilvLilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (42/40);  minor base C (0/1);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950291‑3950554
                                                                                                                                 |                                                                                                                                      
ttCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGc                                                                                                                               >  1:600811/1‑139 (MQ=255)
ttCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGc                                                                                                                               >  1:171871/1‑139 (MQ=255)
  ccTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCt                                                                                                                             >  1:1270533/1‑139 (MQ=255)
  ccTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCt                                                                                                                             >  2:171463/1‑139 (MQ=255)
  ccTTCGAACAAGATGCAACAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATACCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCt                                                                                                                             >  1:926654/1‑139 (MQ=255)
     tCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTaa                                                                                                                          <  1:955258/139‑1 (MQ=255)
                   aGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGAcccc                                                                                                            <  2:886765/139‑1 (MQ=255)
                         gACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    <  1:404667/109‑1 (MQ=255)
                         gACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    >  2:404667/1‑109 (MQ=255)
                            aaaaTGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   <  1:623363/139‑1 (MQ=255)
                                     gCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCggg                                                                                          <  2:294115/139‑1 (MQ=255)
                                                    aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    <  2:947345/82‑1 (MQ=255)
                                                    aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    <  1:253771/82‑1 (MQ=255)
                                                    aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    >  1:947345/1‑82 (MQ=255)
                                                    aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    >  2:253771/1‑82 (MQ=255)
                                                    aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   >  2:102597/1‑115 (MQ=255)
                                                    aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   >  1:616336/1‑115 (MQ=255)
                                                    aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   <  2:616336/115‑1 (MQ=255)
                                                    aTTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGaa                                                                                                   <  1:102597/115‑1 (MQ=255)
                                                     ttAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          >  2:1116416/1‑139 (MQ=255)
                                                      tAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTa                                                                         >  2:218689/1‑139 (MQ=255)
                                                             gTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTACAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaa                                                                        <  1:724149/133‑1 (MQ=255)
                                                             gTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTaa                                                                        >  2:724149/1‑133 (MQ=255)
                                                                       gCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                       <  1:1129317/139‑1 (MQ=255)
                                                                        cGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                       <  2:600811/139‑1 (MQ=255)
                                                                         gtggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                      <  1:221792/139‑1 (MQ=255)
                                                                           ggtggtgATTATTATCCCACCGTGCGGGTCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCCCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  2:926654/139‑1 (MQ=255)
                                                                           ggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt                                                                  >  2:269217/1‑125 (MQ=255)
                                                                           ggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt                                                                  <  1:269217/125‑1 (MQ=255)
                                                                           ggtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:171463/139‑1 (MQ=255)
                                                                            gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGa                                                                                                                       >  2:1029298/1‑71 (MQ=255)
                                                                            gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGa                                                                                                                       <  1:1029298/71‑1 (MQ=255)
                                                                            gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  2:8179/1‑138 (MQ=255)
                                                                            gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  1:8179/138‑1 (MQ=255)
                                                                            gtggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGAc                                                   <  2:171871/139‑1 (MQ=255)
                                                                                  attattATCCCACCGTGCGGGGCTGCACTTGGACGATGAAAGGCTTAAAGATCAAGCCTTAACg                                                                                                                        >  2:1123263/1‑64 (MQ=255)
                                                                                  attattATCCCACCGTGCGGGGCTGCACTTGGACGATGAAAGGCTTAAAGATCAAGCCTTAACg                                                                                                                        <  1:1123263/64‑1 (MQ=255)
                                                                                       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         >  1:313450/1‑58 (MQ=255)
                                                                                       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         <  2:313450/58‑1 (MQ=255)
                                                                                          cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                     <  1:218689/139‑1 (MQ=255)
                                                                                           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  >  2:406046/1‑61 (MQ=255)
                                                                                           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTaa                                                                                                                  <  1:406046/61‑1 (MQ=255)
                                                                                              ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt                                                                  >  1:426402/1‑106 (MQ=255)
                                                                                              ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACAt                                                                  <  2:426402/106‑1 (MQ=255)
                                                                                              ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                 >  2:583644/1‑139 (MQ=255)
                                                                                               cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          >  1:782890/1‑97 (MQ=255)
                                                                                               cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          <  2:782890/97‑1 (MQ=255)
                                                                                                gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                               >  2:509688/1‑139 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    >  1:295559/1‑37 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGAt                                                                                                                                    <  2:295559/37‑1 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         >  2:326204/1‑48 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAAc                                                                                                                         <  1:326204/48‑1 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGATCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                   <  1:347324/134‑1 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGATCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                   >  2:347324/1‑134 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          >  1:956212/1‑95 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCtt                                                                          <  2:956212/95‑1 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  2:1006848/117‑1 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  1:1006848/1‑117 (MQ=255)
                                                                                                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAtt                              >  2:632818/1‑139 (MQ=255)
                                                                                                   cGGGGCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                            >  2:1255692/1‑139 (MQ=255)
                                                                                                    ggggCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                           <  2:993391/139‑1 (MQ=255)
                                                                                                    ggggCTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                           <  1:406483/139‑1 (MQ=255)
                                                                                                        cTGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                       >  1:782502/1‑139 (MQ=255)
                                                                                                         tGCACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                      >  2:353317/1‑139 (MQ=255)
                                                                                                           cACTTGGACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    >  1:988024/1‑139 (MQ=255)
                                                                                                           cACTTGGACGAGGAAAGGCTTAAAGATAAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    >  2:802899/1‑139 (MQ=255)
                                                                                                                ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCgaggag                                                        <  1:1311534/98‑1 (MQ=255)
                                                                                                                ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCgaggag                                                        >  2:1311534/1‑98 (MQ=255)
                                                                                                                ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    >  1:1240733/1‑102 (MQ=255)
                                                                                                                ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                    <  2:1240733/102‑1 (MQ=255)
                                                                                                                ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt               <  2:247990/139‑1 (MQ=255)
                                                                                                                ggACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt               <  1:1116416/139‑1 (MQ=255)
                                                                                                                ggACGAGGAAAGGCTGAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACACAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt               >  2:898107/1‑139 (MQ=255)
                                                                                                                 gACGAGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc              >  2:416386/1‑139 (MQ=255)
                                                                                                                     aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    >  2:762576/1‑129 (MQ=255)
                                                                                                                     aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    <  1:286526/129‑1 (MQ=255)
                                                                                                                     aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    <  1:762576/129‑1 (MQ=255)
                                                                                                                     aGGAAAGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                    >  2:286526/1‑129 (MQ=255)
                                                                                                                        aaaGGCTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgg       <  1:416386/139‑1 (MQ=255)
                                                                                                                             cTTAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac  <  2:955879/139‑1 (MQ=255)
                                                                                                                              ttAAAGATCAAGCCTTAACGAACTAAGACCCCCGCa                                                                                                        >  1:693554/1‑36 (MQ=255)
                                                                                                                              ttAAAGATCAAGCCTTAACGAACTAAGACCCCCGCa                                                                                                        <  2:58874/36‑1 (MQ=255)
                                                                                                                              ttAAAGATCAAGCCTTAACGAACTAAGACCCCCGCa                                                                                                        >  1:58874/1‑36 (MQ=255)
                                                                                                                                 |                                                                                                                                      
TTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAGTGATTAGCCTGGTCGTGATTAGCGTGGTGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAAC  >  NC_000913/3950291‑3950554

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