Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,573,528 A→G pseudogene (1115/1503 nt) yjiT → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,573,5280AG98.8% 295.1 / NA 84pseudogene (1115/1503 nt)yjiTpseudogene
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (31/52);  minor base C (1/0);  total (32/52)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCAATGTTGCTGGTT  >  NC_000913/4573385‑4573671
                                                                                                                                                  |                                                                                                                                                 
aGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGTGGTTTTTTGGGGGATTAAAATTTATGGCAgggga                                                                                                                                            >  1:668439/1‑147 (MQ=255)
  ccTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaaa                                                                                                                                          <  2:39458/150‑1 (MQ=255)
      tGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                                                                     <  2:212268/151‑1 (MQ=255)
      tGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                                                                     <  1:1104194/151‑1 (MQ=255)
      tGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                                                                     <  2:1034457/151‑1 (MQ=255)
      tGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                                                                     <  1:883856/151‑1 (MQ=255)
      tGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATATCGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                                                                     <  2:658747/151‑1 (MQ=255)
             ctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                                                                              <  1:606932/151‑1 (MQ=255)
             ctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCCGGGTGGGGGGGTGGTTTTTTGGGGGATTTAAATTTTTGGGAGGGGAATTGTCATTgggg                                                                                                                              >  1:690228/1‑149 (MQ=255)
               ggAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGTAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGAAGATTGTTGGGAGCATTTAAATTTGTGGGAGGGGAATTTCCATTtgtggg                                                                                                                            >  2:1018662/1‑149 (MQ=255)
               ggAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGATTGTCGTTAGTGGCTCGGGGCGGGGGGGGGGATGTTTGGGGGTTTTAACTTGATGGGAGGGGAAATGCCAATgggtgg                                                                                                                            >  1:70426/1‑147 (MQ=255)
                gAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCCGGGTGGAGGGGGGGTTGTTTGGGGGATTAAAATTTGTGGGAGGGGAATTGCCAtttgggggg                                                                                                                           >  2:741470/1‑144 (MQ=255)
                 aaaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAgg                                                                                                                          <  2:197720/150‑1 (MQ=255)
                           ggCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCCGGTTTGTCGTTAGTGGCTCTGGCTGGAGGGGGGATTTTTTGGAGCATTTAAACTTGTGGCCGGGGAAATTTCCTTGGGGGGGGGGCCGTGaa                                                                                                                 >  1:928450/1‑150 (MQ=255)
                           ggCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGGTGATTGTTGGGGGGATTAAACTTTGTGGCGGGGGAAATGCCATTggggggggggcctttaa                                                                                                                 >  1:320975/1‑135 (MQ=255)
                            gCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTa                                                                                                                  <  1:960618/147‑1 (MQ=255)
                                cAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaaa                                                                                                                                          <  2:621814/120‑1 (MQ=255)
                                   ttCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCg                                                                                                        <  2:323621/150‑1 (MQ=255)
                                         tACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                                                                                 <  2:58593/151‑1 (MQ=255)
                                          aCGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCACGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCa                                                                                                >  2:854620/1‑151 (MQ=255)
                                            gTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGGTGGGGGGGGGGTTTTTTGGGGGATTAAACCTTATGGCAGGGGAAATGGCCTTGGGGAGGGGGCGGTTAACTTTCTCGTTGATGgaga                                                                                              >  2:539048/1‑151 (MQ=255)
                                              aaaaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATc                                                                                            <  1:199693/151‑1 (MQ=255)
                                              aaaaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATc                                                                                            <  1:210454/151‑1 (MQ=255)
                                                   cGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATgg                                                                                       <  1:165938/151‑1 (MQ=255)
                                                   cGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATgg                                                                                       <  1:1154840/151‑1 (MQ=255)
                                                     aggTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                                                                      <  1:1060569/148‑1 (MQ=255)
                                                     aaGTGAGATTTGGCAGGCAAAATCCATCGGCGGGTTTGTCGGTAGTGGCTCAGGCTTGGGGGGTGGTTGTTGGGAGCCATTAAAATGATGGCTTTGAAAATTCCATTGGGGAGGGGGCGGTCAACTTTATCCTTGGTGTAGAACaaaatc                                                                                      >  1:479290/1‑146 (MQ=255)
                                                     aaGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                                                                      <  1:19920/150‑1 (MQ=255)
                                                     aaGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCt                                                                                     <  2:207135/151‑1 (MQ=255)
                                                       gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGTGGGGGGGTGGTTGTTTGGGGCATTTAACTTGGTGGGCGGGGAATTTTCCTTTGGGGGGGGGCGTTTACCTTCTTCTTTGTTGTAGATCAAAGGgtgt                                                                                   >  1:1044874/1‑151 (MQ=255)
                                                          agaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGGGCTGGGGGGGGGGTTGTTTGGGGGATTTAAATTTTTGGCGGGGAAAATTTCCTTTGGGGGGGGCCGTTAAACTTCATCCTTTGTGGAGATCAATGGGTGTTat                                                                                >  1:881990/1‑150 (MQ=255)
                                                             tttGGCAGGCAAAATCCAGCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                                                                             <  2:275713/151‑1 (MQ=255)
                                                              ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAggg                                                                            <  1:754586/151‑1 (MQ=255)
                                                                 gcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGAc                                                                          <  2:169028/150‑1 (MQ=255)
                                                                 gcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGAc                                                                          <  2:92915/150‑1 (MQ=255)
                                                                           tCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc                                                               <  2:440320/151‑1 (MQ=255)
                                                                           tCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc                                                               <  2:466931/151‑1 (MQ=255)
                                                                           tCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc                                                               <  2:620539/151‑1 (MQ=255)
                                                                           tCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc                                                               <  2:1002600/151‑1 (MQ=255)
                                                                           tCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc                                                               <  1:1204339/151‑1 (MQ=255)
                                                                           tCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc                                                               <  1:774703/151‑1 (MQ=255)
                                                                             cATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg                                                             <  2:982478/151‑1 (MQ=255)
                                                                             cATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg                                                             <  2:321567/151‑1 (MQ=255)
                                                                               tCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtg                                                           <  1:741786/151‑1 (MQ=255)
                                                                                       tttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagc                                                   <  2:910932/151‑1 (MQ=255)
                                                                                       tttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagc                                                   <  2:881655/151‑1 (MQ=255)
                                                                                       tttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagc                                                   <  2:838463/151‑1 (MQ=255)
                                                                                       tttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAag                                                    <  1:525414/150‑1 (MQ=255)
                                                                                          gTCGTTAGTGGCTCGTGCTGGGGGGGTGATTGTTTGGGGCATTTAACTTGATGGCCGTGAAATTGTCATTGGTGAGGGGGCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGTGcagcag                                                 >  2:1009314/1‑150 (MQ=255)
                                                                                            cGTTAGTGGCTCGTGCTGGAGGGGTGATTGTTGGGAGGATTAAAATTGATGGCAGTGAAAATGCCATTGGTGGGGGGCCGTTAAACTTCATCGTTGATGCAGATCAATGGCTGTTTAAGGGACAGGGTTCTTTCAGTGGGCGAAGCAGCGa                                              >  2:987278/1‑151 (MQ=255)
                                                                                               tAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGt                                           <  2:378278/151‑1 (MQ=255)
                                                                                                  tGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTc                                         <  1:285044/150‑1 (MQ=255)
                                                                                                  tGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTc                                         <  1:642762/150‑1 (MQ=255)
                                                                                                   ggCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGGCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGtt                                          >  2:899395/1‑148 (MQ=255)
                                                                                                        gTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATtgt                                  >  2:470040/1‑151 (MQ=255)
                                                                                                         tGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGGGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATtgtg                                 >  1:811685/1‑151 (MQ=255)
                                                                                                          gCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATtgtg                                 >  1:227623/1‑150 (MQ=255)
                                                                                                             ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTc                              >  2:639108/1‑150 (MQ=255)
                                                                                                             ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTc                              >  2:261481/1‑150 (MQ=255)
                                                                                                                 ggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCt                           <  1:854954/149‑1 (MQ=255)
                                                                                                                  gatgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCgg                        >  1:207875/1‑151 (MQ=255)
                                                                                                                    tgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGa                      <  1:257725/151‑1 (MQ=255)
                                                                                                                    tgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGa                      <  1:261542/151‑1 (MQ=255)
                                                                                                                    tgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGa                      >  1:421410/1‑151 (MQ=255)
                                                                                                                    tgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGa                      <  2:114531/151‑1 (MQ=255)
                                                                                                                      atTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGata                    >  2:303743/1‑151 (MQ=255)
                                                                                                                         gttgGGAGCATTAAACTTGATGGCAGTGACATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGataata                 >  2:998207/1‑151 (MQ=255)
                                                                                                                         gttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGataata                 >  2:617228/1‑151 (MQ=255)
                                                                                                                         gttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGataata                 >  1:925251/1‑151 (MQ=255)
                                                                                                                          ttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAg                <  1:30299/151‑1 (MQ=255)
                                                                                                                          ttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAg                <  2:456198/151‑1 (MQ=255)
                                                                                                                          ttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAg                <  1:756896/151‑1 (MQ=255)
                                                                                                                           tgtgAGCATTAACCTGGATGGCAGTGAACTTGCCATTGGCGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGGGTGGGAACCAGCGATTTTCTGAATGGGCTCCATCTGGATAGTAGc               <  2:404053/148‑1 (MQ=255)
                                                                                                                           tgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAg                >  2:648939/1‑150 (MQ=255)
                                                                                                                           tgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGc               <  1:893758/151‑1 (MQ=255)
                                                                                                                              gAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCaa             <  1:649817/150‑1 (MQ=255)
                                                                                                                                gCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCAATGt          >  2:289549/1‑151 (MQ=255)
                                                                                                                                gCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCAATGt          >  2:1240915/1‑151 (MQ=255)
                                                                                                                                gCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCAATGt          >  1:776536/1‑151 (MQ=255)
                                                                                                                                  cATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACGGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCAATg           <  1:899189/149‑1 (MQ=255)
                                                                                                                                  cATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCAATGt          <  1:1243738/150‑1 (MQ=255)
                                                                                                                                   aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCAATGTTg        <  1:593577/151‑1 (MQ=255)
                                                                                                                                       aaaCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGTCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTTTGCTCCCTCGGGATAATAGCAATGTTg        >  2:274814/1‑148 (MQ=255)
                                                                                                                                          cTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCAATGTTGCTGGtt  >  1:1205177/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGGATAATAGCAATGTTGCTGGTT  >  NC_000913/4573385‑4573671

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