Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,465,487 T→C 5.4% pseudogene (96/2513 nt) ydbA → pseudogene, autotransporter homolog; interrupted by IS2 and IS30

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,465,4870TC 5.4% 232.1 / 3.0 75pseudogene (96/2513 nt)ydbApseudogene, autotransporter homolog; interrupted by IS2 and IS30
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base C (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

ATTGGTCGGGTTAGTGCATCCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCA  >  NC_000913/1465339‑1465632
                                                                                                                                                    |                                                                                                                                                 
aTTGGTCGGGTTAGTGCATCCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGt                                                                                                                                                   <  1:89677/149‑1 (MQ=255)
    gTCGGGTTAGTGCATCCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTgaa                                                                                                                                                >  1:106897/1‑148 (MQ=255)
       gggTTAGTGCATCCTGCTTTCTTTATATTCGCGAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTgaaaag                                                                                                                                             <  1:46903/148‑1 (MQ=255)
          ttAGTGCATCCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTgaaaagaaa                                                                                                                                          >  7:362634/1‑148 (MQ=255)
            aGTGCATCCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGaaaaa                                                                                                                                        <  8:273088/148‑1 (MQ=255)
             gTGCATCCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATa                                                                                                                                      >  7:100193/1‑149 (MQ=255)
               gCATCCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATAc                                                                                                                                     >  7:326763/1‑148 (MQ=255)
                      gCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGa                                                                                                                             <  5:203012/149‑1 (MQ=255)
                       ctttctttATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGGCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGAc                                                                                                                            <  6:222813/149‑1 (MQ=255)
                            tttATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTc                                                                                                                       >  8:8330/1‑149 (MQ=255)
                            tttATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTc                                                                                                                       >  6:370931/1‑149 (MQ=255)
                            tttATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTc                                                                                                                       >  4:194220/1‑149 (MQ=255)
                              tataTGCGCCAGTCGGTTTTATTATGCAAAGGAAAACTCTATTTTCGGCCTGTATTGCATTAGCTCTGAGTGGGCAGGGTTGGGCGTCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATCCCAATGTGACTTGt                                                                                                                        <  2:124302/146‑1 (MQ=255)
                                      ccAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAATGAAAAATACCAATGTGACTTGTCCGGCAGAcc                                                                                                              >  5:91896/1‑148 (MQ=255)
                                      ccAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGAccc                                                                                                             <  3:119034/149‑1 (MQ=255)
                                       cAGAAGGATTTTTTATGCAAATGTAAAATCTTTTGTCGGCCTGTATTGCATTAGCTCTGTGTGGGCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCTGCAGAccc                                                                                                             <  6:91896/148‑1 (MQ=255)
                                         gAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAg                                                                                                           >  5:340198/1‑148 (MQ=255)
                                         gAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGCGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAg                                                                                                           >  3:23463/1‑148 (MQ=255)
                                          aaGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTgaaaagaaa                                                                                                                                          <  5:362682/116‑1 (MQ=255)
                                          aaGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTgaaaagaaa                                                                                                                                          >  6:362682/1‑116 (MQ=255)
                                                  attatGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCa                                                                                                  <  2:349130/148‑1 (MQ=255)
                                                  attatGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAg                                                                                                 >  6:135312/1‑149 (MQ=255)
                                                   ttatGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGCGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGt                                                                                                <  4:23463/149‑1 (MQ=255)
                                                         aaaGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGCAAAATACCAATGGGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCgg                                                                                            >  6:290591/1‑147 (MQ=255)
                                                         aaaGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGAc                                                                                                                            >  1:9042/1‑115 (MQ=255)
                                                         aaaGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGAc                                                                                                                            <  2:9042/115‑1 (MQ=255)
                                                         aaaGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGCAAACTCAGTCCgg                                                                                            >  7:80532/1‑147 (MQ=255)
                                                         aaaGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCgg                                                                                            >  7:148748/1‑147 (MQ=255)
                                                         aaaGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCgg                                                                                            >  3:142900/1‑147 (MQ=255)
                                                              aaaaCTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCt                                                                                     <  7:8330/149‑1 (MQ=255)
                                                                   tctATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCa                                                                                                                  <  7:54851/115‑1 (MQ=255)
                                                                   tctATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATAACAATGTGACTTGTCCGGCa                                                                                                                  >  8:54851/1‑115 (MQ=255)
                                                                     tATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAAc                                                                                <  5:370931/147‑1 (MQ=255)
                                                                     tATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGc                                                                              >  5:231640/1‑149 (MQ=255)
                                                                     tATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGCGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGc                                                                              >  2:31910/1‑149 (MQ=255)
                                                                       ttGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCtt                                                                            <  3:175655/149‑1 (MQ=255)
                                                                       ttGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCt                                                                             >  1:43476/1‑148 (MQ=255)
                                                                       ttGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGCGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCt                                                                             <  1:31910/148‑1 (MQ=255)
                                                                         gTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGAc                                                                                                                            >  4:251963/1‑99 (MQ=255)
                                                                         gTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGAc                                                                                                                            <  3:251963/99‑1 (MQ=255)
                                                                         gTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTa                                                                           <  6:340198/148‑1 (MQ=255)
                                                                          tCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTAc                                                                          <  8:100193/148‑1 (MQ=255)
                                                                             gCCTGTATTGCATTAGCTCTGAGTGGGCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCctc                                                                      <  4:35839/149‑1 (MQ=255)
                                                                              ccTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCctct                                                                     <  8:362634/149‑1 (MQ=255)
                                                                                   aTTGCATTAGCTCTGAGTGGTCAGGGTTTGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCGTAAACGCTTACCCTCTGAAt                                                                 <  1:384943/148‑1 (MQ=255)
                                                                                      gCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTa                                                                           >  3:13445/1‑135 (MQ=255)
                                                                                      gCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTa                                                                           <  4:13445/135‑1 (MQ=255)
                                                                                         ttAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAAt                                                                                                                                       >  5:202450/1‑72 (MQ=255)
                                                                                         ttAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAAt                                                                                                                                       <  6:202450/72‑1 (MQ=255)
                                                                                           aGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCt                                                         <  7:167744/148‑1 (MQ=255)
                                                                                            gCTCTGAGTTGTCAGGGTTGTGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCtt                                                        <  5:42820/148‑1 (MQ=255)
                                                                                             ctctGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAAt                                                                                                                                 <  3:145342/74‑1 (MQ=255)
                                                                                             ctctGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAAt                                                                                                                                 >  4:145342/1‑74 (MQ=255)
                                                                                              tctGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAg                                                     >  2:53118/1‑149 (MQ=255)
                                                                                              tctGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTATCCg                                                                                                                     <  4:71738/85‑1 (MQ=255)
                                                                                              tctGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTATCCg                                                                                                                     >  3:71738/1‑85 (MQ=255)
                                                                                                tGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGa                                                                   >  7:128110/1‑133 (MQ=255)
                                                                                                tGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGa                                                                   <  8:128110/133‑1 (MQ=255)
                                                                                                  aGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTc                                                   <  8:326763/147‑1 (MQ=255)
                                                                                                   gTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAg                                                                                                 <  1:276723/100‑1 (MQ=255)
                                                                                                   gTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAg                                                                                                 >  2:276723/1‑100 (MQ=255)
                                                                                                             ttGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGCAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGa                                      <  8:80532/149‑1 (MQ=255)
                                                                                                             ttGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGa                                      <  2:106897/149‑1 (MQ=255)
                                                                                                             ttGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGa                                      <  8:148748/149‑1 (MQ=255)
                                                                                                             ttGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGGGAAAAGAAAACTACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGa                                      <  2:53712/149‑1 (MQ=255)
                                                                                                               gggCTGCAGATATCAGAGAGGTAGAAACCAACACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGTAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGAt                                     <  5:135312/148‑1 (MQ=255)
                                                                                                                    gcAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATg                               <  3:194220/149‑1 (MQ=255)
                                                                                                                             acaGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCcac                      >  4:172536/1‑149 (MQ=255)
                                                                                                                                     aGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTa                                                                           >  6:138256/1‑88 (MQ=255)
                                                                                                                                     aGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTa                                                                           <  5:138256/88‑1 (MQ=255)
                                                                                                                                             ccacAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGtt       >  2:144285/1‑148 (MQ=255)
                                                                                                                                              cacaGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTaa     >  1:116793/1‑149 (MQ=255)
                                                                                                                                                caGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATc   <  1:53118/149‑1 (MQ=255)
                                                                                                                                                 aGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCa  >  4:306306/1‑149 (MQ=255)
                                                                                                                                                   gTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCcaca                     <  7:166895/128‑1 (MQ=255)
                                                                                                                                                   gTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCcaca                     >  8:166895/1‑128 (MQ=255)
                                                                                                                                                    tGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCct                                                                       <  4:131748/77‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ATTGGTCGGGTTAGTGCATCCGGCTTTCTTTATATTCGCCAGAAGGATTTATTATGCAAAGGAAAACTCTATTGTCGGCCTGTATTGCATTAGCTCTGAGTGGTCAGGGTTGGGCGGCAGATATCACAGAGGTAGAAACCACCACAGGTGAAAAGAAAAATACCAATGTGACTTGTCCGGCAGACCCAGGAAAACTCAGTCCGGAAGAGCTTAAACGCTTACCCTCTGAATGCTCTCCTTTAGTCGAACAAAACCTGATGCCATGGCTTTCCACAGGCGCTGCTGCGTTAATCA  >  NC_000913/1465339‑1465632

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