Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,815,810 Δ1 bp intergenic (‑42/+24) 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.96.7% 265.1 / 2.8 60intergenic (‑41/+25)pyrE/rphorotate phosphoribosyltransferase/ribonuclease PH (defective);enzyme; Degradation of RNA; RNase PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (37/21);  minor base A (2/0);  total (39/21)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACT  >  NC_000913/3815663‑3815953
                                                                                                                                                  |                                                                                                                                                
tGAAGAAATAGGGGCTTTTGCGCCCGGAGTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAgcga                                                                                                                                               <  2:51806/150‑4 (MQ=255)
   agaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACtc                                                                                                                                           <  2:314328/150‑1 (MQ=255)
    gaaATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACtcat                                                                                                                                         >  1:680183/1‑151 (MQ=255)
      aaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACtcatta                                                                                                                                       <  1:167615/151‑3 (MQ=255)
      aaTAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACtcatca                                                                                                                                       <  1:167592/151‑1 (MQ=255)
              cTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCtt                                                                                                                               >  1:44745/1‑151 (MQ=255)
               ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCtt                                                                                                                               >  2:823459/1‑150 (MQ=255)
               ttttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCtt                                                                                                                               >  1:940071/1‑150 (MQ=255)
                tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTa                                                                                                                              <  1:83222/150‑1 (MQ=255)
                tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTa                                                                                                                              <  1:250134/150‑1 (MQ=255)
                tttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTa                                                                                                                              <  2:254948/150‑1 (MQ=255)
                 ttGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTaa                                                                                                                             <  1:382558/150‑1 (MQ=255)
                     gcCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAAt                                                                                                                         >  1:772175/1‑150 (MQ=255)
                      cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATc                                                                                                                        >  2:822315/1‑150 (MQ=255)
                      cccGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCa                                                                                                                       <  1:448521/151‑1 (MQ=255)
                         ggATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGtt                                                                                                                    >  2:831917/1‑151 (MQ=255)
                          gATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGttt                                                                                                                   <  1:281232/151‑1 (MQ=255)
                              tCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGcc                                                                                                                >  2:552618/1‑150 (MQ=255)
                               cAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCa                                                                                                               <  2:453998/150‑1 (MQ=255)
                                aGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgc                                                                                                             <  2:32067/151‑1 (MQ=255)
                                aGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgc                                                                                                             >  2:335922/1‑151 (MQ=255)
                                 gCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAgcg                                                                                                            >  1:144676/1‑151 (MQ=255)
                                    tAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgccg                                                                                                         >  1:3360/1‑151 (MQ=255)
                                    tAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGcgcc                                                                                                          >  2:970015/1‑150 (MQ=255)
                                        cTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCtt                                                                                                     >  1:276654/1‑151 (MQ=255)
                                            ccAAACTTTAACCCCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACGCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGc                                                                                                 <  1:933497/151‑1 (MQ=255)
                                             cAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCg                                                                                                >  2:489334/1‑151 (MQ=255)
                                               aaCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCACTGATATTGTTTCATGCCTTCGCTCCTCATCTTACGTGTCTACAGACACAAAAAA‑GCGACTCATCAGTCGCCTTAAAAAGCAGTGTGCCAGCGCCGCCTTCTTCGTc                                                                                              <  2:394663/151‑1 (MQ=255)
                                                aCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCg                                                                                             >  1:694129/1‑151 (MQ=255)
                                                aCTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GAAACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCg                                                                                             >  2:535584/1‑151 (MQ=255)
                                                    tAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTAc                                                                                         >  1:680718/1‑151 (MQ=255)
                                                      acacCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACa                                                                                        >  1:435992/1‑150 (MQ=255)
                                                           tgcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgat                                                                                  >  1:710621/1‑151 (MQ=255)
                                                            gcttgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTACAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgatt                                                                                 >  2:708908/1‑151 (MQ=255)
                                                               tgctAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGgattc                                                                                <  2:922112/149‑1 (MQ=255)
                                                                   aaGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTcc                                                                          >  1:993321/1‑151 (MQ=255)
                                                                       gcAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTc                                                                      <  2:762548/151‑1 (MQ=255)
                                                                         aaaTTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCgg                                                                    <  2:497991/151‑1 (MQ=255)
                                                                             tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTCCTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGTCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCggg                                                                   <  2:702984/148‑1 (MQ=255)
                                                                             tCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGcc                                                                 >  2:790247/1‑150 (MQ=255)
                                                                              cAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACCAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGcc                                                                 >  2:894770/1‑149 (MQ=255)
                                                                                       ggCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaacaa                                                      <  2:232192/151‑1 (MQ=255)
                                                                                          gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaaca                                                       >  2:913711/1‑147 (MQ=255)
                                                                                          gcTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGGCCCCTTAAAAATCAGTTTTCCAACGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAACcaa                                                         >  1:634034/1‑145 (MQ=255)
                                                                                          gcTGATATGGTTTCATGCCTTCGCACCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCcaaca                                                       >  2:514254/1‑147 (MQ=255)
                                                                                            tGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                  >  2:814419/1‑150 (MQ=255)
                                                                                             gATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATg                                                  >  2:60970/1‑149 (MQ=255)
                                                                                                 tGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGACGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAg                                             >  1:810833/1‑150 (MQ=255)
                                                                                                     ttCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCtt                                        >  1:946806/1‑151 (MQ=255)
                                                                                                           ccTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATggg                                  >  1:446090/1‑151 (MQ=255)
                                                                                                                  tcctcATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACgg                           >  2:413798/1‑151 (MQ=255)
                                                                                                                  tcctcATCTTACGTTTCTACAGACAAAAAAAA‑GCGACTCATCAGGCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTGCAATGGATTCGATTCCCCTCGGTCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTg                                <  2:39252/146‑1 (MQ=255)
                                                                                                                   ccttagcTTACTTCTGTACAGCCAAAAAAAA‑ACGACTCATCAGATGCCTGAAAAATCAGTTTGCCAGCGCCGACTTCTGCGTCGCTACAATGGATGCGATTCCCCACGGGCCAGAGCCAGCAGTAGGAGTAGCTCTTCAGGGGAGAGCGGc                          <  2:664086/145‑1 (MQ=255)
                                                                                                                   cctcATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGc                          >  1:382631/1‑151 (MQ=255)
                                                                                                                   cctcATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCACTTCATGGGTGATCGGc                          >  2:881339/1‑151 (MQ=255)
                                                                                                                     tcATCTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTc                        >  2:8501/1‑151 (MQ=255)
                                                                                                                         cTTACTTTTCTACAGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCt                    >  1:423919/1‑151 (MQ=255)
                                                                                                                                      aGACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTcccc        >  2:85338/1‑150 (MQ=255)
                                                                                                                                       gACAAAAAAAA‑GCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTg      <  2:522610/151‑1 (MQ=255)
                                                                                                                                         acaaaaaaaaGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGc     <  2:44741/149‑1 (MQ=255)
                                                                                                                                            aaaaaaAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACt  >  1:144088/1‑151 (MQ=255)
                                                                                                                                            aaaaaaAGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACt  >  2:293666/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TGAAGAAATAGGGGCTTTTGCGCCCGGATTTCAGCGTAAACTCGCCAAACTTTAACACCTGCTTGCTAAGCGCAAATTCAATAAACTGGCGCTGATATGGTTTCATGCCTTCGCTCCTCATCTTACTTTTCTACAGACAAAAAAAAGGCGACTCATCAGTCGCCTTAAAAATCAGTTTGCCAGCGCCGCCTTCTGCGTCGCTACAATGGATTCGATTCCCCTCGGGCCAGAGCCAACAAGATGAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACT  >  NC_000913/3815663‑3815953

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑