Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,809 Δ1 bp 32.2% intergenic (‑41/+25) pyrE ← / ← rph orotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,8090G.32.2% 103.3 / 84.0 59intergenic (‑41/+25)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (13/22);  new base . (15/4);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 4.48e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-01

AGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCC  >  NC_000913/3815672‑3815941
                                                                                                                                          |                                                                                                                                    
aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGc                                                                                                                                    <  1:641895/140‑1 (MQ=255)
             cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTAGc                                                                                                                      >  2:852711/1‑140 (MQ=255)
             cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTAGc                                                                                                                      >  2:1608623/1‑140 (MQ=255)
                  aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGACGCCTTaa                                                                                                                 >  2:14587/1‑140 (MQ=255)
                                    cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg                                                                                                <  2:66335/140‑1 (MQ=255)
                                     aaaCTTTACCCCCGGCTGGCTAAGCGAAAATACAATCAACTGGCGCTTATATGGTTACATGCCTCCGCTCCACAACTTACTTTTCAACAGAAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgcc                                                                                              <  2:395956/140‑1 (MQ=255)
                                         ttcaccacCTGCTTGCTAAGCGCAAATTCAACACCCGGGCGCTGATTAGGTTACATGCCTTCCCTCCTCATCTTACTTGTCAACAGACACAACAAAGGCGACTCATCAGTCGACTTAAAAATCCGTTTGCCAGCgccgcc                                                                                           <  2:569220/135‑1 (MQ=255)
                                            aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTC                                                                                        >  1:1121324‑M1/1‑140 (MQ=255)
                                            aaCACCTGATTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGATTCATGCATTCGCTCCTCAACTTACTTTACTAAAGACAAAAAAAAGGCGACTCAACAGTCGCCTTAAAAATCAGATAGCCAGCGCCGCCTTC                                                                                        >  2:1554073‑M1/1‑140 (MQ=255)
                                                ccTGCTTGCTAAGCGCAAATTCAATACACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCCTCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                    <  1:1127830‑M1/140‑2 (MQ=255)
                                                  tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAATGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTc                                                                                  >  2:343555/1‑140 (MQ=255)
                                                  tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAATGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTc                                                                                  >  2:1620754/1‑140 (MQ=255)
                                                   gcttgctAAGCGCTAATTCAATAAACTGGCGCTGTTATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAATATAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGc                                                                                >  2:585848/1‑140 (MQ=255)
                                                   gcttgctAAGCGCAAATTCATTAAACTGGCGCTGATATGCTATCATGCCTTCGCTCCTCATCTTACTTTACTAAAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGACGc                                                                                >  2:123541/1‑140 (MQ=255)
                                                   gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                 >  2:1483420‑M1/1‑136 (MQ=255)
                                                   gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                 >  2:1482902‑M1/1‑136 (MQ=255)
                                                        ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                            <  1:1376723‑M1/140‑10 (MQ=255)
                                                           aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                                         <  1:1143155‑M1/140‑13 (MQ=255)
                                                             cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:159078‑M1/140‑15 (MQ=255)
                                                             cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:1645078‑M1/140‑15 (MQ=255)
                                                             cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:1408314‑M1/140‑15 (MQ=255)
                                                             cgcAAATTCAAAAAACTCGCGCTGATATGGTTTCATGCATTCGCTCCTCAACTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCATAAAAAACACATTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:1488099‑M1/140‑15 (MQ=255)
                                                              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:103151‑M1/140‑15 (MQ=255)
                                                                  aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                  >  1:564504/1‑140 (MQ=255)
                                                                   ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt                                                                <  1:585848/140‑1 (MQ=255)
                                                                   ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAAAAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt                                                                <  2:1563781/140‑1 (MQ=255)
                                                                       aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcccc                                                            >  2:1595934/1‑140 (MQ=255)
                                                                       aTAAACTGGCGCTGATATGGTTTCATGCCTACGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAAGGATTCGATTcccc                                                            >  2:1638715/1‑140 (MQ=255)
                                                                                cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCag                                                   >  1:479744/1‑140 (MQ=255)
                                                                                cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTACCCTCGGGCCag                                                   >  2:748934/1‑140 (MQ=255)
                                                                                     atatGCTTTCATGCCTTCCATCCCCATCTTTATACACAACACACAAAAAAAAGCGACACATCAGTCGCCATAAAAATCAGATTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCca                                              <  2:1064587/140‑1 (MQ=255)
                                                                                      tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                              <  2:1605866‑M1/140‑40 (MQ=255)
                                                                                          gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCACTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                                          <  2:475359‑M1/140‑44 (MQ=255)
                                                                                           tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                         >  2:1572218‑M1/1‑96 (MQ=255)
                                                                                           tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                         >  2:1572291‑M1/1‑96 (MQ=255)
                                                                                             tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                       <  1:574911‑M1/140‑47 (MQ=255)
                                                                                               aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                     >  2:1472898‑M1/1‑92 (MQ=255)
                                                                                               aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                    >  1:209417/1‑140 (MQ=255)
                                                                                               aTGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAGGCGACTCATCAGTCGCCTTAAAAAACAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                    >  2:749273/1‑140 (MQ=255)
                                                                                                   cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                 <  1:1507043/140‑1 (MQ=255)
                                                                                                   cTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                 <  1:1170613/140‑1 (MQ=255)
                                                                                                    ttCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGc                                <  1:1484148/140‑1 (MQ=255)
                                                                                                            tcATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                       >  1:265240/1‑140 (MQ=255)
                                                                                                             cATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                      >  2:290329/1‑140 (MQ=255)
                                                                                                             cATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                      >  2:682222/1‑140 (MQ=255)
                                                                                                             cATCTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                      >  2:682094/1‑140 (MQ=255)
                                                                                                               tCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcg                     <  2:1121324‑M1/140‑65 (MQ=255)
                                                                                                                cTTACTTTTCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAACATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                   <  2:1341647/140‑1 (MQ=255)
                                                                                                                cTTACTTTTCTACACACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                   <  2:495853/140‑1 (MQ=255)
                                                                                                                     ttttCTACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagt              >  2:1216528‑M1/1‑69 (MQ=255)
                                                                                                                      tttCTACACAAAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt             <  1:414995/140‑1 (MQ=255)
                                                                                                                          tACAGACAATAAAAAGCGACTCATCAGTCCCCTTAAAAATCAGTTTGCCAGCTCCGCCATCTGCTTCGCTAAAATTGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGAAGCTCTTCATGGGTGAACGGCTCGcc         >  2:832334/1‑140 (MQ=255)
                                                                                                                          tACAGACAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc         <  2:328994/140‑1 (MQ=255)
                                                                                                                          aaCAGACAAAAAAAAGCGACTCATCAGTCGCCGTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGcc         >  2:832325/2‑140 (MQ=255)
                                                                                                                           aCAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttca         <  2:123212‑M1/140‑77 (MQ=255)
                                                                                                                             cgaCAAAAAAAAGGCGTCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtctgtctctgcggcagagtcttcaac       <  2:62892‑M1/139‑79 (MQ=255)
                                                                                                                             aGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac       <  2:63036‑M1/140‑79 (MQ=255)
                                                                                                                              gACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg      >  1:1544928‑M1/1‑61 (MQ=255)
                                                                                                                                cAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgta    >  1:1560565‑M1/1‑59 (MQ=255)
                                                                                                                                 acaaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg    <  2:1340670/138‑1 (MQ=255)
                                                                                                                                 acaaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg    >  1:1159994/3‑140 (MQ=255)
                                                                                                                                 acaaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg    <  2:429871/138‑1 (MQ=255)
                                                                                                                                 acaaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTGTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCCTGGGTGAACGGCTCGCCTTCTg    <  2:789034/138‑1 (MQ=255)
                                                                                                                                 acaaaaaaaaGCGACTCATCAGTCGCCTAAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCGACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTg    <  2:1341132/138‑1 (MQ=255)
                                                                                                                                   aaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGcc  <  2:629925/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
AGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCC  >  NC_000913/3815672‑3815941

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: