Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,220,920 C→T 5.6% pseudogene (1320/2574 nt) ycgH → pseudogene; putative ATP‑binding component of a transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,220,9200CT 5.6% 186.2 / 2.2 71pseudogene (1320/2574 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base T (2/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GAAGCTTTCAGAAAATACAATCTGGAATATGAAAGATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGG  >  NC_000913/1220773‑1221067
                                                                                                                                                   |                                                                                                                                                   
gAAGCTTTCAGAAAATACAATCTGGAATATGAAAGATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTTTGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCg                                                                                                                                                    >  1:1027906/1‑149 (MQ=255)
   gCTTTCAGAAAATACAATCTGGAATATGAAAGATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACACAAGGAAAAACATTAACCGAAAAAGGTAATTATGTCGGta                                                                                                                                                 >  2:600195/1‑149 (MQ=255)
         aGAAAATACAATCTGGAATATGAAAGATGACTCCGTTGTTACCCATCTGACTAATTCATACAGTATTATAAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAATGTAATTATGTCGGAAATAATg                                                                                                                                           >  2:832092/1‑149 (MQ=255)
               aaCAATCTGGAATATGAAAGATGACTCCGTTGTTACCCATCTGACAAATTCAGACAGTATTAACAATCTATCGTCTGATGACGGTCAACCATTTACCCAACGACAAACAAACACCGTAAAACGTAATTATGTCGGTAATAATGGTCAGc                                                                                                                                     <  2:393707/148‑1 (MQ=255)
                aCAATCTGGAATATGAAAGATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAAAATTAAACGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCt                                                                                                                                    >  1:177368/1‑149 (MQ=255)
                aCAATCTGGAATATGAAAGATGACTCCGTTGGTACCCATCTGACTTATTCAGACAGTCTTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGATAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCt                                                                                                                                    >  2:945731/1‑149 (MQ=255)
                 cAATCTGGAATATGAACGATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTAGTATCATGATGGTCAAACATTTCCCCAAGGAAAAAGATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGGCAGCtt                                                                                                                                   <  2:695313/149‑1 (MQ=255)
                  aaTCTGGAATATGAAAGATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTa                                                                                                                                  <  1:495680/149‑1 (MQ=255)
                        gAATATGAAAGATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATcc                                                                                                                            <  1:86392/149‑1 (MQ=255)
                         aatatgaaatatgaCTCAGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATGTTTATTATGATGATGTTCAAACATTTAACCAAGTAAAAACATTAAAAGTAAAAGTTAATTTTGTCGGTAATTTTGGTCAGCTTAATATCCg                                                                                                                           <  2:80204/149‑1 (MQ=255)
                              gAAAGATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTTGGTAATAATGGTCAGCTTAATATCCGCACCg                                                                                                                      >  2:370074/1‑149 (MQ=255)
                                  gATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTAtt                                                                                                                  <  1:587814/149‑1 (MQ=255)
                                    tGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAg                                                                                                                <  2:351580/149‑1 (MQ=255)
                                    tGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAg                                                                                                                <  1:342608/149‑1 (MQ=255)
                                     taCTCCGTTCTTACCCATCTGACTAATTCAGACAGTTTTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAgg                                                                                                               <  1:492312/148‑1 (MQ=255)
                                        tCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGtga                                                                                                            <  1:107444/149‑1 (MQ=255)
                                          ctttgttACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATACTGGTCAGCTTAATATCCGCACCGTATTAGGggatg                                                                                                          <  1:704696/147‑1 (MQ=255)
                                          cGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGTTAATTATGTCGGTAAAAATGTTCAGCTTAATATCCGCACCGTATTAGGtgttt                                                                                                          >  2:675227/1‑146 (MQ=255)
                                          cGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGtgatg                                                                                                          >  1:475/1‑149 (MQ=255)
                                            ttgttACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAAGATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTgatgat                                                                                                        <  1:64675/149‑1 (MQ=255)
                                                  cccATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCg                                                                                                  >  2:247385/1‑149 (MQ=255)
                                                   ccATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCgg                                                                                                 >  1:756945/1‑149 (MQ=255)
                                                   ccATATGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCgg                                                                                                 <  1:234628/149‑1 (MQ=255)
                                                      tCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTa                                                                                              <  2:944925/149‑1 (MQ=255)
                                                         gACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCGACgg                                                                                           >  2:465681/1‑149 (MQ=255)
                                                          aCTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGa                                                                                          <  1:237456/149‑1 (MQ=255)
                                                          aCTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGa                                                                                          >  2:379180/1‑149 (MQ=255)
                                                           cTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGAc                                                                                         >  1:846706/1‑149 (MQ=255)
                                                            tAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACa                                                                                        >  2:631711/1‑149 (MQ=255)
                                                               ttCAGACAGTATTATCTATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGAc                                                                                     <  1:864881/149‑1 (MQ=255)
                                                               ttCAGACAGTATTATCAATCTTTCGTATGCTGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTCAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGGATTAGGGGATGATAAATCGGCTACGGACAGAc                                                                                     <  1:913483/149‑1 (MQ=255)
                                                               ttCAGACAGTATTATCAATCTTTCGTATGATTATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCAGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGAc                                                                                     >  2:104175/1‑149 (MQ=255)
                                                               ttCAGACAGTATTATCAATCTTTCGTATGATGACGGTCAAACATTTACCCACGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTACTATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGAc                                                                                     <  1:917336/149‑1 (MQ=255)
                                                                  aGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTATTTATGTCGGTAATATTGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATGTGCTACGGACAGagtaa                                                                                  >  1:341171/1‑145 (MQ=255)
                                                                  aGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTa                                                                                  >  1:347616/1‑149 (MQ=255)
                                                                    aCAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGGTACGGAAGGACTTAtt                                                                                >  1:701187/1‑149 (MQ=255)
                                                                     cAGTATTATCAATCTTTCGTATGATGATGGTCTAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTAttg                                                                               <  1:357457/149‑1 (MQ=255)
                                                                         attatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGa                                                                           <  2:475/149‑1 (MQ=255)
                                                                           tatCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCTTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCAGCACCGTATTAGGTGATGATAAATCGGCTACGGACATACTTATTGTTGAgg                                                                         >  1:77113/1‑149 (MQ=255)
                                                                            atCAATCTTTCGTATGATGATTGTCAAACATTTACCCAATGAAAAACATTCACCGTATAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCAGGGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAggg                                                                        <  2:312334/149‑1 (MQ=255)
                                                                            atCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAggg                                                                        <  2:29586/149‑1 (MQ=255)
                                                                             tCAATTTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCAGCACCGTATTATGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGt                                                                       <  2:445307/149‑1 (MQ=255)
                                                                             tCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGAAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCCGCTACGGCCAGACTTATTGTTGAGGGt                                                                       >  1:591386/1‑149 (MQ=255)
                                                                               aaTCTTTCGTATGATGATTGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCTGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTTTTGTTGAGGGTaa                                                                     <  2:302229/149‑1 (MQ=255)
                                                                               aaTCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTTGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTaa                                                                     >  2:565460/1‑149 (MQ=255)
                                                                                aTCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAAt                                                                    >  2:75628/1‑149 (MQ=255)
                                                                                        tATGATGATGGTCAAACATTTACCCACGGAAAAACATTTACCGTAAAAGGTAATTATGTCGGTAAACATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGGTTCGGACAGACTTATTGTTGAGGGTAATACTTCggg                                                            <  2:483455/149‑1 (MQ=255)
                                                                                         atgatgatgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTTTGTCGGTATTAAGGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGAATTATTGTTGAGGGTAATACTTCGGGt                                                           >  1:50616/1‑149 (MQ=255)
                                                                                         atgatgatgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTaa                                                                     >  1:326516/1‑139 (MQ=255)
                                                                                         atgatgatgGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTaa                                                                     <  2:326516/139‑1 (MQ=255)
                                                                                                  gTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCTGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTAcc                                                  <  2:124864/149‑1 (MQ=255)
                                                                                                  gTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTAcc                                                  <  1:85834/149‑1 (MQ=255)
                                                                                                      aaCATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGACGGTAATACTTCGGGTTAAACTACGGTCt                                              >  2:368215/1‑149 (MQ=255)
                                                                                                            tccccAAGGAAAAACATTAACCGTAAACGGCAATTATGCCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGAGGATAACTCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGa                                        <  2:701187/147‑1 (MQ=255)
                                                                                                             aCCCAAGGAAAAACATTACCCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTAATAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTACGGGTTCAACTACCGTCTATGAGaa                                       >  2:85726/1‑149 (MQ=255)
                                                                                                               ccAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGaaaa                                     >  1:445181/1‑149 (MQ=255)
                                                                                                                  aGGAAAAACAATAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATg                                  <  2:978344/149‑1 (MQ=255)
                                                                                                                    gAAAAACATTAACCGTAAAAGGTAATTAAGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCt                                <  1:136939/149‑1 (MQ=255)
                                                                                                                     aaaaaCATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTg                               <  1:75628/149‑1 (MQ=255)
                                                                                                                     aaaaaCATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTg                               <  2:846706/149‑1 (MQ=255)
                                                                                                                        aaCATTAACCGTAAAAGGTAATTATGTTGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTggag                            <  1:370074/149‑1 (MQ=255)
                                                                                                                           aTTAACCGTAAAACGTAATTATGTCGGTAATAATGGTCAGCTGAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAAAACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGCAGGaa                         >  2:1023550/1‑149 (MQ=255)
                                                                                                                                     aaaGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAgcggcgcggc               <  1:832092/149‑1 (MQ=255)
                                                                                                                                     aaaGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGATTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGcc              >  1:75364/1‑149 (MQ=255)
                                                                                                                                       aGGTAATTATGTTGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCa             <  1:565460/149‑1 (MQ=255)
                                                                                                                                       aGGAAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGGTCAACTACCGTCTATGTGAAAACTGCTGGAGGAAGCGGCGCGGCCa             <  1:430519/149‑1 (MQ=255)
                                                                                                                                          tAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCACCTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGc          <  2:756945/149‑1 (MQ=255)
                                                                                                                                               aTGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGCGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAAc     >  1:888606/1‑149 (MQ=255)
                                                                                                                                                 gTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATGGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACgg   <  2:652426/149‑1 (MQ=255)
                                                                                                                                                  tCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGTTAATAATTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGTAGGAATCGTCGCTGCCAAGCTAAAAggg  >  1:644859/1‑149 (MQ=255)
                                                                                                                                                  tCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAAAggg  >  1:648183/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAAGCTTTCAGAAAATACAATCTGGAATATGAAAGATGACTCCGTTGTTACCCATCTGACTAATTCAGACAGTATTATCAATCTTTCGTATGATGATGGTCAAACATTTACCCAAGGAAAAACATTAACCGTAAAAGGTAATTATGTCGGTAATAATGGTCAGCTTAATATCCGCACCGTATTAGGTGATGATAAATCGGCTACGGACAGACTTATTGTTGAGGGTAATACTTCGGGTTCAACTACCGTCTATGTGAAAAATGCTGGAGGAAGCGGCGCGGCCACGCTAAACGGG  >  NC_000913/1220773‑1221067

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