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

AGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCC  >  NC_000913/3815672‑3815941
                                                                                                                                          |                                                                                                                                    
aGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGc                                                                                                                                    <  1:641895/140‑1 (MQ=255)
             cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTAGc                                                                                                                      >  2:852711/1‑140 (MQ=255)
             cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTAGc                                                                                                                      >  2:1608623/1‑140 (MQ=255)
                  aTTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGACGCCTTaa                                                                                                                 >  2:14587/1‑140 (MQ=255)
                                    cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg                                                                                                <  2:66335/140‑1 (MQ=255)
                                     aaaCTTTACCCCCGGCTGGCTAAGCGAAAATACAATCAACTGGCGCTTATATGGTTACATGCCTCCGCTCCACAACTTACTTTTCAACAGA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgcc                                                                                              <  2:395956/140‑1 (MQ=255)
                                         ttcaccacCTGCTTGCTAAGCGCAAATTCAACACCCGGGCGCTGATTAGGTTACATGCCTTCCCTCCTCATCTTACTTGTCAACAGAC‑ACAACAAAGGCGACTCATCAGTCGACTTAAAAATCCGTTTGCCAGCgccgcc                                                                                           <  2:569220/135‑1 (MQ=255)
                                            aaCACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTC                                                                                        >  1:1121324‑M1/1‑140 (MQ=255)
                                            aaCACCTGATTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGATTCATGCATTCGCTCCTCAACTTACTTTACTAAAGAC‑AAAAAAAAGGCGACTCAACAGTCGCCTTAAAAATCAGATAGCCAGCGCCGCCTTC                                                                                        >  2:1554073‑M1/1‑140 (MQ=255)
                                                ccTGCTTGCTAAGCGCAAATTCAATACACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCCTCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCc                                                                                    <  1:1127830‑M1/140‑2 (MQ=255)
                                                  tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAATGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTc                                                                                  >  2:343555/1‑140 (MQ=255)
                                                  tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAATGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTc                                                                                  >  2:1620754/1‑140 (MQ=255)
                                                   gcttgctAAGCGCTAATTCAATAAACTGGCGCTGTTATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAATATAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGc                                                                                >  2:585848/1‑140 (MQ=255)
                                                   gcttgctAAGCGCAAATTCATTAAACTGGCGCTGATATGCTATCATGCCTTCGCTCCTCATCTTACTTTACTAAAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGACGc                                                                                >  2:123541/1‑140 (MQ=255)
                                                   gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                 >  2:1483420‑M1/1‑136 (MQ=255)
                                                   gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtc                                                                                 >  2:1482902‑M1/1‑136 (MQ=255)
                                                        ctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctg                                                                            <  1:1376723‑M1/140‑10 (MQ=255)
                                                           aGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcac                                                                         <  1:1143155‑M1/140‑13 (MQ=255)
                                                             cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:159078‑M1/140‑15 (MQ=255)
                                                             cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:1645078‑M1/140‑15 (MQ=255)
                                                             cgcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:1408314‑M1/140‑15 (MQ=255)
                                                             cgcAAATTCAAAAAACTCGCGCTGATATGGTTTCATGCATTCGCTCCTCAACTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCATAAAAAACACATTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:1488099‑M1/140‑15 (MQ=255)
                                                              gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcactt                                                                       <  2:103151‑M1/140‑15 (MQ=255)
                                                                  aTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcga                                                                  >  1:564504/1‑140 (MQ=255)
                                                                   ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt                                                                <  1:585848/140‑1 (MQ=255)
                                                                   ttCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGA‑‑AAAAAAAACGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattcgatt                                                                <  2:1563781/140‑1 (MQ=255)
                                                                       aTAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcccc                                                            >  2:1595934/1‑140 (MQ=255)
                                                                       aTAAACTGGCGCTGATATGGTTTCATGCCTACGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAAAGGATTCGATTcccc                                                            >  2:1638715/1‑140 (MQ=255)
                                                                                cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCag                                                   >  1:479744/1‑140 (MQ=255)
                                                                                cgcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTACCCTCGGGCCag                                                   >  2:748934/1‑140 (MQ=255)
                                                                                     atatGCTTTCATGCCTTCCATCCCCATCTTTATACACAACACAC‑AAAAAAAA‑GCGACACATCAGTCGCCATAAAAATCAGATTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCca                                              <  2:1064587/140‑1 (MQ=255)
                                                                                      tatGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcat                                              <  2:1605866‑M1/140‑40 (MQ=255)
                                                                                          gTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCACTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcact                                          <  2:475359‑M1/140‑44 (MQ=255)
                                                                                           tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                         >  2:1572218‑M1/1‑96 (MQ=255)
                                                                                           tttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcacta                                         >  2:1572291‑M1/1‑96 (MQ=255)
                                                                                             tCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacg                                       <  1:574911‑M1/140‑47 (MQ=255)
                                                                                               aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgtt                                     >  2:1472898‑M1/1‑92 (MQ=255)
                                                                                               aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                    >  1:209417/1‑140 (MQ=255)
                                                                                               aTGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑‑AAAAAAAGGCGACTCATCAGTCGCCTTAAAAAACAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAg                                    >  2:749273/1‑140 (MQ=255)
                                                                                                   cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                 <  1:1507043/140‑1 (MQ=255)
                                                                                                   cTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                 <  1:1170613/140‑1 (MQ=255)
                                                                                                    ttCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGc                                <  1:1484148/140‑1 (MQ=255)
                                                                                                            tcATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATgg                       >  1:265240/1‑140 (MQ=255)
                                                                                                             cATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                      >  2:290329/1‑140 (MQ=255)
                                                                                                             cATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                      >  2:682222/1‑140 (MQ=255)
                                                                                                             cATCTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                      >  2:682094/1‑140 (MQ=255)
                                                                                                               tCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcg                     <  2:1121324‑M1/140‑65 (MQ=255)
                                                                                                                cTTACTTTTCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAACATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                   <  2:1341647/140‑1 (MQ=255)
                                                                                                                cTTACTTTTCTACACAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGa                   <  2:495853/140‑1 (MQ=255)
                                                                                                                     ttttCTACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagt              >  2:1216528‑M1/1‑69 (MQ=255)
                                                                                                                      tttCTACACA‑‑AAAAAAAAAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCt             <  1:414995/140‑1 (MQ=255)
                                                                                                                          tACAGAC‑AATAAAAA‑GCGACTCATCAGTCCCCTTAAAAATCAGTTTGCCAGCTCCGCCATCTGCTTCGCTAAAATTGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGAAGCTCTTCATGGGTGAACGGCTCGcc         >  2:832334/1‑140 (MQ=255)
                                                                                                                          tACAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGcc         <  2:328994/140‑1 (MQ=255)
                                                                                                                          aaCAGAC‑AAAAAAAA‑GCGACTCATCAGTCGCCGTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGACCGGCTCGcc         >  2:832325/2‑140 (MQ=255)
                                                                                                                           aCAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttca         <  2:123212‑M1/140‑77 (MQ=255)
                                                                                                                             cgaC‑AAAAAAAAGGCGTCTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtctgtctctgcggcagagtcttcaac       <  2:62892‑M1/139‑79 (MQ=255)
                                                                                                                             aGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaac       <  2:63036‑M1/140‑79 (MQ=255)
                                                                                                                              gAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacg      >  1:1544928‑M1/1‑61 (MQ=255)
                                                                                                                                c‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCcgtcccctgcacttcaatgatgcgcccgtcttcggtcatcactacgttcatgtcggtctctgcggcagagtcttcaacgta    >  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)
                                                                                                                                          |                                                                                                                                    
AGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGAC‑AAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCC  >  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: ‑