Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 737,765 A→G pseudogene (941/1137 nt) ybfL → pseudogene, DDE domain transposase family;putative factor; Not classified; putative receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913737,7650AG100.0% 254.5 / NA 71pseudogene (941/1137 nt)ybfLpseudogene, DDE domain transposase family;putative factor; Not classified; putative receptor protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (34/37);  total (34/37)

TTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCA  >  NC_000913/737621‑737907
                                                                                                                                                |                                                                                                                                              
ttCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGt                                                                                                                                          >  2:794877/1‑151 (MQ=255)
   ggTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaat                                                                                                                                       >  2:141996/1‑151 (MQ=255)
     tCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaat                                                                                                                                       >  1:324966/1‑149 (MQ=255)
      cAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGACAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaat                                                                                                                                       >  2:213884/1‑148 (MQ=255)
        ataataGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTGGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTaatgaatg                                                                                                                                  <  2:551205/151‑1 (MQ=255)
          aataGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAAtgaatgaa                                                                                                                                <  1:357571/151‑1 (MQ=38)
            taGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAatgaatga                                                                                                                                 >  1:639527/1‑148 (MQ=38)
                 gAACAAACGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACt                                                                                                                         >  1:856415/1‑151 (MQ=255)
                   aCAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGc                                                                                                                       >  2:835300/1‑151 (MQ=38)
                    cAAAAGAAAGAGCCATAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTgg                                                                                                                                           >  2:659515/1‑130 (MQ=255)
                     aaaaGAAAGAGCCATAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGt                                                                                                                                          <  1:659473/130‑1 (MQ=255)
                      aaagaaagAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCaaa                                                                                                                    <  2:722429/151‑1 (MQ=38)
                             gagCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATaagaa                                                                                                             <  1:702197/151‑1 (MQ=38)
                                ccAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATaagaag                                                                                                            >  1:735919/1‑149 (MQ=25)
                                 cAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAgg                                                                                                         <  2:727045/151‑1 (MQ=25)
                                    aaaTGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaa                                                                                                       >  2:84933/1‑150 (MQ=18)
                                    aaaTGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGaa                                                                                                       >  2:607890/1‑150 (MQ=18)
                                      aTGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAAcg                                                                                                    <  1:775304/151‑1 (MQ=25)
                                      aTGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAAcg                                                                                                    <  2:705639/151‑1 (MQ=25)
                                           ggTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCa                                                                                               <  1:498263/151‑1 (MQ=38)
                                              cAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGaa                                                                                            <  1:234925/151‑1 (MQ=38)
                                               aGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAat                                                                                           <  2:555171/151‑1 (MQ=38)
                                                   attatATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAattatt                                                                                       >  2:365545/1‑151 (MQ=38)
                                                        atCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAg                                                                                  <  1:794874/151‑1 (MQ=38)
                                                          cAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAggg                                                                                <  2:525741/151‑1 (MQ=38)
                                                           aGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAggg                                                                                >  1:228722/1‑150 (MQ=38)
                                                           aGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGa                                                                               <  2:807705/151‑1 (MQ=38)
                                                                   tgATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGcac                                                                       >  2:845421/1‑151 (MQ=38)
                                                                    gATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGcac                                                                       >  1:130863/1‑150 (MQ=38)
                                                                    gATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGca                                                                        <  2:699130/149‑1 (MQ=38)
                                                                    gATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGca                                                                        <  1:758810/149‑1 (MQ=38)
                                                                    gATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGcac                                                                       <  1:388244/149‑1 (MQ=255)
                                                                        tAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGc                                                                  >  2:84084/1‑151 (MQ=255)
                                                                         aaCCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGCAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCt                                                                 >  1:349804/1‑151 (MQ=255)
                                                                           ccGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTAt                                                               >  1:325810/1‑151 (MQ=255)
                                                                            cGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTAtt                                                              >  2:823203/1‑151 (MQ=255)
                                                                              cAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAgg                                                                                 <  1:12103/130‑1 (MQ=255)
                                                                              cAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAgg                                                                                 >  2:12124/1‑130 (MQ=255)
                                                                                gAAAGGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAata                                                          <  2:213748/151‑1 (MQ=255)
                                                                                 aaaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAatat                                                         <  1:542280/151‑1 (MQ=255)
                                                                                  aaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATAtt                                                        <  2:577390/151‑1 (MQ=39)
                                                                                  aaaGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATAtt                                                        >  2:434893/1‑151 (MQ=39)
                                                                                    aGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATAtttt                                                      <  2:151746/151‑1 (MQ=38)
                                                                                       tCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTaa                                                    >  1:650436/1‑150 (MQ=255)
                                                                                        cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACg                                                  <  1:512800/151‑1 (MQ=255)
                                                                                          ccACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGaa                                                <  2:297792/151‑1 (MQ=255)
                                                                                           cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAAt                                               >  1:317339/1‑151 (MQ=255)
                                                                                           cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAAt                                               <  1:118788/151‑1 (MQ=255)
                                                                                           cacaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTATGGGAGTGGAAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAAt                                               >  1:708226/1‑151 (MQ=255)
                                                                                            acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATg                                              <  1:387609/151‑1 (MQ=255)
                                                                                            acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATg                                              >  1:636602/1‑151 (MQ=255)
                                                                                            acaGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATg                                              <  2:96813/151‑1 (MQ=255)
                                                                                             caGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATg                                              >  1:500513/1‑150 (MQ=255)
                                                                                                cAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATaa                                          <  2:344445/151‑1 (MQ=255)
                                                                                                  aTCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAgg                                        <  2:612353/151‑1 (MQ=255)
                                                                                                    ccGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGt                                       <  1:157452/150‑1 (MQ=255)
                                                                                                    ccGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTa                                      <  1:775225/151‑1 (MQ=255)
                                                                                                            aCTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGc                              >  2:241089/1‑151 (MQ=38)
                                                                                                                 cACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTTGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGt                         <  2:394758/151‑1 (MQ=255)
                                                                                                                 cACGTGGAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGt                         <  1:46878/151‑1 (MQ=38)
                                                                                                                    gTGGAGAATAAGCTTCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAAAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCGCATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTaa                      <  1:532299/151‑1 (MQ=255)
                                                                                                                      ggAGAATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTa                       >  2:843608/1‑148 (MQ=38)
                                                                                                                         gaATAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGt                 >  1:128451/1‑151 (MQ=38)
                                                                                                                          aaTAAGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGt                 >  2:511637/1‑150 (MQ=38)
                                                                                                                             aaGCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAg              <  1:583994/150‑1 (MQ=38)
                                                                                                                               gCTGCACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAACATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATg           >  1:484139/1‑151 (MQ=255)
                                                                                                                                   cACTGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGa        >  2:486561/1‑150 (MQ=38)
                                                                                                                                      tGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaaa      >  2:74019/1‑149 (MQ=38)
                                                                                                                                      tGGCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGaaa      >  2:179143/1‑149 (MQ=38)
                                                                                                                                        gCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagca  <  2:785106/151‑1 (MQ=38)
                                                                                                                                        gCGTCTGGGCGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAagca  <  1:316354/151‑1 (MQ=38)
                                                                                                                                                |                                                                                                                                              
TTCGGTCAATAATAGCAGAACAAAAGAAAGAGCCAGAAATGACGGTCAGATATTATATCAGTTCTGCTGATTTAACCGCAGAAAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTAATGAATGAAGACGACTGCAAAATAAGAAGAGGAAACGCCGCAGAATTATTTTCAGGGATACGGCACATCGCTATTAATATTTTAACGAATGATAAGGTATTCAAGGCAGGGTTAAGACGTAAGATGCGAAAAGCA  >  NC_000913/737621‑737907

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