Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,509,571 C→A 81.6% pseudogene (93/315 nt) insO → pseudogene, IS911 transposase B;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,509,5710CA81.6% 263.8 / 68.6 114pseudogene (93/315 nt)insOpseudogene, IS911 transposase B;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related
Reads supporting (aligned to +/- strand):  ref base C (13/8);  new base A (49/44);  total (62/52)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.43e-01

Only 100 of 118 total aligned reads displayed.

AGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCA  >  NC_000913/4509425‑4509711
                                                                                                                                                  |                                                                                                                                            
aGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAActgatg                                                                                                                                            >  4:398688/1‑149 (MQ=255)
aGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAActgatg                                                                                                                                            <  1:476697/149‑1 (MQ=255)
         aaTGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAg                                                                                                                                   >  1:806838/1‑149 (MQ=255)
          aTGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGc                                                                                                                                  <  1:16772/149‑1 (MQ=255)
              cATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCaa                                                                                                                              >  3:33602/1‑149 (MQ=255)
              cATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCaa                                                                                                                              >  3:258844/1‑149 (MQ=255)
               aTGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAAt                                                                                                                             >  1:344517/1‑149 (MQ=255)
                tGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATg                                                                                                                            <  2:873605/149‑1 (MQ=255)
                tGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATg                                                                                                                            >  4:75800/1‑149 (MQ=255)
                tGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAAAAGCCAAATg                                                                                                                            >  1:852578/1‑149 (MQ=255)
                 gggAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATgg                                                                                                                           >  3:27459/1‑149 (MQ=255)
                  ggAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGa                                                                                                                          >  3:171840/1‑149 (MQ=255)
                   gAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCCGCGTGGAAActgc                                                                                                                                              <  4:882709/128‑1 (MQ=255)
                        ccGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGctt                                                                                                                    >  3:537817/1‑149 (MQ=255)
                        ccGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGctt                                                                                                                    >  4:281643/1‑149 (MQ=255)
                        ccGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGctt                                                                                                                    >  3:748844/1‑149 (MQ=255)
                              tAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAg                                                                                                              <  4:20017/149‑1 (MQ=255)
                               aaGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAgg                                                                                                             >  1:144785/1‑149 (MQ=255)
                                  cccGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGt                                                                                                          <  3:674899/149‑1 (MQ=255)
                                   ccGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTTCTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTc                                                                                                         <  4:241441/149‑1 (MQ=255)
                                   ccGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTc                                                                                                         <  4:27459/149‑1 (MQ=255)
                                    cGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCttcttc                                                                                                                >  2:538565/1‑141 (MQ=255)
                                    cGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCttcttc                                                                                                                <  1:538565/141‑1 (MQ=255)
                                    cGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGACCAATGGAGCGCTTCTTCAGGAGTCt                                                                                                        >  4:737219/1‑149 (MQ=255)
                                         gggTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTgaaga                                                                                                   >  4:204418/1‑149 (MQ=255)
                                          ggTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTgaagaa                                                                                                  >  1:350724/1‑149 (MQ=255)
                                           gTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAAc                                                                                                 <  1:674420/149‑1 (MQ=255)
                                           gTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAAc                                                                                                 <  2:550173/149‑1 (MQ=255)
                                             gATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAg                                                                                                              <  2:696929/134‑1 (MQ=255)
                                             gATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAg                                                                                                              >  1:696929/1‑134 (MQ=255)
                                             gATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGa                                                                                               >  3:883905/1‑149 (MQ=255)
                                              aTGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGaa                                                                                              >  4:362506/1‑149 (MQ=255)
                                                 ttCCAAGCGATCAAGTCAGACATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATgg                                                                                           >  2:114894/1‑149 (MQ=255)
                                                 ttCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATgg                                                                                           <  2:806838/149‑1 (MQ=255)
                                                 ttCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATgg                                                                                           >  1:151986/1‑149 (MQ=255)
                                                 ttCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATgg                                                                                           <  1:392474/149‑1 (MQ=255)
                                                 ttCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATgg                                                                                           <  4:883905/149‑1 (MQ=255)
                                                  tCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATggg                                                                                          >  2:830984/1‑149 (MQ=255)
                                                   ccAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGt                                                                                         <  2:409894/149‑1 (MQ=255)
                                                     aaGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCaa                                                                                                                              >  4:24956/1‑110 (MQ=255)
                                                     aaGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCaa                                                                                                                              <  3:24956/110‑1 (MQ=255)
                                                     aaGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGc                                                                                       >  4:828483/1‑149 (MQ=255)
                                                       gCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCa                                                                                     <  4:345112/149‑1 (MQ=255)
                                                         gATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGc                                                                                   >  1:352160/1‑149 (MQ=255)
                                                           tCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGa                                                                                 >  1:311650/1‑149 (MQ=255)
                                                              aGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACgg                                                                              >  4:783312/1‑149 (MQ=255)
                                                              aGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACgg                                                                              >  1:760778/1‑149 (MQ=255)
                                                               ggCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACggg                                                                             >  3:412815/1‑149 (MQ=255)
                                                                gCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGc                                                                            >  2:504107/1‑149 (MQ=255)
                                                                 cAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCt                                                                           >  1:524863/1‑149 (MQ=255)
                                                                 cAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCt                                                                           <  4:714048/149‑1 (MQ=255)
                                                                  aGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTa                                                                          <  4:33602/149‑1 (MQ=255)
                                                                  aGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTa                                                                          >  3:892604/1‑149 (MQ=255)
                                                                      aTTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGct                                                                                                                     >  4:465064/1‑102 (MQ=255)
                                                                      aTTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGct                                                                                                                     <  3:465064/102‑1 (MQ=255)
                                                                      aTTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTa                                                                      >  3:876599/1‑149 (MQ=255)
                                                                      aTTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTa                                                                      >  3:628584/1‑149 (MQ=255)
                                                                       ttATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTaa                                                                     <  4:379224/149‑1 (MQ=255)
                                                                         ataCGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGc                                                                   <  4:628584/149‑1 (MQ=255)
                                                                          taCGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGcc                                                                                      <  3:76572/129‑1 (MQ=255)
                                                                          taCGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGcc                                                                                      >  4:76572/1‑129 (MQ=255)
                                                                             gAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCa                                                               <  3:204418/149‑1 (MQ=255)
                                                                               gcaggcagTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATgg                                                                                           <  4:289926/119‑1 (MQ=255)
                                                                               gcaggcagTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGc                                                             <  1:504107/149‑1 (MQ=255)
                                                                                 aggcagTTCCGGCAGTTACTGTGGCGATACGGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTACAGGAGTCTGGAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGa                                                           >  3:888714/1‑149 (MQ=255)
                                                                                 aggcagTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGa                                                           <  3:55935/149‑1 (MQ=255)
                                                                                      gTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAg                                                      >  3:138544/1‑149 (MQ=255)
                                                                                       ttCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGc                                                     >  1:82907/1‑149 (MQ=255)
                                                                                          cGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGCCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCa                                                  >  2:467958/1‑149 (MQ=255)
                                                                                           ggCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCAc                                                 >  1:565751/1‑149 (MQ=255)
                                                                                              aGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAATCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCa                                              >  1:905210/1‑149 (MQ=255)
                                                                                                ttACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAAt                                            >  1:147904/1‑149 (MQ=255)
                                                                                                ttACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAAt                                            >  2:791667/1‑149 (MQ=255)
                                                                                                      tgGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGa                                      <  2:144785/149‑1 (MQ=255)
                                                                                                      tgGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGa                                      <  3:42830/149‑1 (MQ=255)
                                                                                                      tgGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGa                                      <  4:537817/149‑1 (MQ=255)
                                                                                                      tgGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGa                                      <  2:311650/149‑1 (MQ=255)
                                                                                                            tACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTatat                                >  1:394368/1‑149 (MQ=255)
                                                                                                              ccGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCg                              >  4:852394/1‑149 (MQ=255)
                                                                                                               cGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGt                             >  4:86080/1‑149 (MQ=255)
                                                                                                                ggATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGtt                            <  3:651912/149‑1 (MQ=255)
                                                                                                                ggATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGtt                            >  1:365506/1‑149 (MQ=255)
                                                                                                                ggATCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGtt                            <  4:655018/149‑1 (MQ=255)
                                                                                                                 gATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTg                           >  3:54543/1‑149 (MQ=255)
                                                                                                                  aTCAGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTgg                          >  1:817372/1‑149 (MQ=255)
                                                                                                                     aGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATggg                                                                                          >  4:340060/1‑82 (MQ=255)
                                                                                                                     aGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATggg                                                                                          <  3:340060/82‑1 (MQ=255)
                                                                                                                     aGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGata                       <  3:75800/149‑1 (MQ=255)
                                                                                                                     aGGCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGa                         >  4:685439/1‑147 (MQ=255)
                                                                                                                       gCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAg                                                                    <  3:127372/102‑1 (MQ=255)
                                                                                                                       gCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAg                                                                    >  4:127372/1‑102 (MQ=255)
                                                                                                                       gCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCaa                                             >  1:733721/1‑125 (MQ=255)
                                                                                                                       gCAGAGTATGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATAGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCaa                                             <  2:733721/125‑1 (MQ=255)
                                                                                                                               tGAGTCGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCa             <  3:86080/149‑1 (MQ=255)
                                                                                                                                   tCGGCGGGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTaa         <  1:258449/149‑1 (MQ=255)
                                                                                                                                    cGGCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATg                                                                        <  3:299398/85‑1 (MQ=255)
                                                                                                                                      gCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGAc      <  3:240156/149‑1 (MQ=255)
                                                                                                                                      gCGTGGAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAATCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGAc      <  2:905210/149‑1 (MQ=255)
                                                                                                                                          ggAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCa  <  4:748844/149‑1 (MQ=255)
                                                                                                                                          ggAAACTGATGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCa  <  2:71532/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
AGCACTGGAAATGGCATGGGAAACCCGTGGTAAGCCCGTCGGGGTGATGTTCCAAGCGATCAAGGCAGTCATTATACGAGCAGGCAGTTCCGGCAGTTACTGTGGCGATACCGGATCAGGCAGAGTATGAGTCGGCGTGGAAACTGCTGGGATAACAGCCCAATGGAGCGCTTCTTCAGGAGTCTGAAGAACGAATGGGTGCCAGCGACGGGCTATGTAAGCTTCAGCGATGCAGCTCACGCAATAACGGACTATATCGTTGGATATTACAGCGCACTAAGACCGCA  >  NC_000913/4509425‑4509711

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: