Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,633,634 C→T E38K (GAG→AAG)  yjjX ← non‑canonical purine NTP phosphatase, ITPase/XTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,633,6340CT98.4% 188.8 / NA 63E38K (GAG→AAG) yjjXnon‑canonical purine NTP phosphatase, ITPase/XTPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/27);  minor base G (1/0);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCA  >  NC_000913/4633488‑4633779
                                                                                                                                                  |                                                                                                                                                 
aCTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGa                                                                                                                                                 <  1:806790/149‑1 (MQ=255)
 cTGAAAGTGCTGTCACCATCGATGCCAGCTTCAAACGCCACCCAAAACTCAGCCTCTCGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGAc                                                                                                                                                <  1:1264253/149‑1 (MQ=255)
   gAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACgg                                                                                                                                              >  1:866926/1‑149 (MQ=255)
    aaaGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACAGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGc                                                                                                                                             >  1:59595/1‑149 (MQ=255)
           tGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACgatg                                                                                                                                      >  1:570334/1‑146 (MQ=255)
           tGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGAt                                                                                                                                      <  1:272952/149‑1 (MQ=255)
           tGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGAt                                                                                                                                      <  2:394610/149‑1 (MQ=255)
            gTCACCATCGATGCCAGCTTCAATCGCCCCCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACCCCGCTCTTGACGGCAAAGGATg                                                                                                                                     >  1:59923/1‑149 (MQ=255)
               aCCGTCGATGCCAGCTTCAATCGCCACCAAAAGCTCACCCGCTGGAACTACAGGGCGGGGATTGGCTACCCGATTTCGTGCGACAGGGCGCGGTTCCTAACAGCCAAAGGGCTGTTCCGCTACTCCGCTCTTGACGGCAACGGATGCaa                                                                                                                                  <  1:843661/149‑1 (MQ=255)
                   tCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATg                                                                                                                              <  1:1174802/149‑1 (MQ=255)
                        gCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCACCGGCTGCAATATGGCAgg                                                                                                                         >  2:240154/1‑149 (MQ=255)
                          cAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCCCTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAAGATGGCAGGAt                                                                                                                       >  1:1129056/1‑149 (MQ=255)
                          cAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCCCTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGAAATATGGCAGGaa                                                                                                                       >  2:66656/1‑148 (MQ=255)
                              ttCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCtt                                                                                                                   >  1:600495/1‑149 (MQ=255)
                              ttCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGAGGCAATATGGCAGGATCCtt                                                                                                                   >  1:424507/1‑149 (MQ=255)
                                 aaTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGAAATATGGCAGGATCCTTCGc                                                                                                                >  2:900550/1‑149 (MQ=255)
                                 aaTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTGGACGGCAACGGATGCAATATGGCAGGATCCTTCGc                                                                                                                >  1:895801/1‑149 (MQ=255)
                                  aTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGcc                                                                                                               <  2:211873/149‑1 (MQ=255)
                                    cGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCATTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGAGGCAATATGGCAGGATCCTTCGCCGa                                                                                                             >  1:31276/1‑149 (MQ=255)
                                     gCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCAGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGaa                                                                                                            >  1:1188733/1‑149 (MQ=255)
                                          ccAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGAtct                                                                                                       >  1:58052/1‑149 (MQ=255)
                                           cAAAAATCAGCCTCTGGAAGTAAACATCGGGCATTGGCTACCCGTTTTCGTGCGCCAGAGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGGACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGAtctc                                                                                                      <  2:1159530/149‑1 (MQ=255)
                                            aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCg                                                                                                     >  2:439850/1‑149 (MQ=255)
                                            aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCg                                                                                                     >  1:797247/1‑149 (MQ=255)
                                                     ccTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGcc                                                                                            >  1:1104051/1‑149 (MQ=255)
                                                      cgcttgAATTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCt                                                                                           <  2:1200885/144‑1 (MQ=255)
                                                        ctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGc                                                                                         <  1:439850/149‑1 (MQ=255)
                                                           gAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGa                                                                                      >  2:343131/1‑149 (MQ=255)
                                                             aGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAAt                                                                                    <  2:690446/149‑1 (MQ=255)
                                                              gTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATg                                                                                   <  1:744775/149‑1 (MQ=255)
                                                                aaaCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGc                                                                                 <  2:578607/149‑1 (MQ=255)
                                                                  aCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGTTCTCGTGAAATGCCTGCAGAATGGCCt                                                                               <  1:7954/149‑1 (MQ=255)
                                                                    ggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGa                                                                             >  2:220403/1‑149 (MQ=255)
                                                                      cggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAAt                                                                           <  1:382238/149‑1 (MQ=255)
                                                                          cATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGCAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATGTTa                                                                       >  2:810232/1‑149 (MQ=255)
                                                                          cATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTa                                                                       <  1:1102219/149‑1 (MQ=255)
                                                                            ttGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGc                                                                     <  1:326595/149‑1 (MQ=255)
                                                                                    ccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCTGGATtgg                                                             >  1:893328/1‑149 (MQ=255)
                                                                                     cGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTTCAGAATGGCCTGAATTTTAGCGGGATtggt                                                            >  1:575578/1‑149 (MQ=255)
                                                                                       gaaacgcgcgcCGGCGCACGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACCCCGCTCTTGACGGCAACGGATGCAATCTGGCAGGCTCCTTCGCGGAAGAGCTCGTGAACTGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggggg                                                          >  1:1222428/8‑149 (MQ=255)
                                                                                        tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCGGCAGAATGGCCTGAATTTTAGCGGGATTggtggt                                                         >  1:672830/1‑149 (MQ=255)
                                                                                                    gcgcgcGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTAGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCGTGAATTTTAGCGGGATTGGTGGTCGCACAGACAAc                                             >  2:1262603/1‑149 (MQ=255)
                                                                                                    gcgcgcGGTTCCTCACTGCCACAGGGCTGTTCGGGTACACCGCACTTGACGGCAACTGATGCAATATGGCAGGATCCTTCGCCGAATATCTCGTGAACTGCCTGCAGAATGGCCTGAATTTTAGCGGGCTTGGTGGTCGCACAGACAAc                                             <  2:859751/149‑1 (MQ=255)
                                                                                                        gcgTTTCCTCACTGCCCAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATACTTCGCCGAAGATCTCGTGAACTGCCTGCAGAATGGCCTGAATTTTAGCGGTATTGGTGGTCGCACAGACAACTTgg                                         <  2:371242/149‑1 (MQ=255)
                                                                                                          gTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAGTTTTAGCGGGATTGGTGGTCGCACAGCCAACTTGGTg                                       >  2:261339/1‑149 (MQ=255)
                                                                                                          gTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTg                                       >  2:787802/1‑149 (MQ=255)
                                                                                                          gTTTCCTCACTGCAAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGAAGCAATATGGCAGGATCCTTAGCCGCAGATCTCATGAAATGCATGAATAATTGCCTTAATTTTAGCTGGTATGGTGGTCTCACAGACAACttgttg                                       >  1:697192/1‑149 (MQ=255)
                                                                                                              ccTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATa                                   <  1:918947/149‑1 (MQ=255)
                                                                                                                tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAAt                                 <  1:987953/149‑1 (MQ=255)
                                                                                                                tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCGACGCCGAAGATCTAGTGACATGCCTTCAGAATGGCCTGAAATTTAGCGGGATTGGAGGTCGAACAGACAACTTGGTGCATAAt                                 >  2:1206005/1‑149 (MQ=255)
                                                                                                                 cACTGCCAATGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGATATGCCTGCAGAATGGCCTGAATTTTTGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                <  2:59595/149‑1 (MQ=255)
                                                                                                                 cACTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                <  1:208214/149‑1 (MQ=255)
                                                                                                                  aCTGCCAAAGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCa                               >  2:673158/1‑149 (MQ=255)
                                                                                                                        aaaGGGCTGTTCCGGTACACCGCTCTTGACGACAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTa                         >  2:80286/1‑149 (MQ=255)
                                                                                                                         aaGGGCTGTTCCGGTACACCGCTCTAGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGACTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTAc                        <  2:666472/149‑1 (MQ=255)
                                                                                                                          aGGGCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCATCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACt                       <  2:600495/149‑1 (MQ=255)
                                                                                                                           gggCTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTc                      <  2:58052/149‑1 (MQ=255)
                                                                                                                              cTGTTCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGa                   >  1:853366/1‑149 (MQ=255)
                                                                                                                                  tCCGGTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAAt               <  2:388213/149‑1 (MQ=255)
                                                                                                                                     ggTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTaa            <  2:19232/149‑1 (MQ=255)
                                                                                                                                      gTACACCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCTGAATGGCCTGAATTTTAGCGTGATTGGTGGTCGCACAGACAAATTGGTTCATAATCAGCATTACTCTGAAAATTaaa           >  2:550304/1‑148 (MQ=255)
                                                                                                                                          acCGCTCTTGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTa       >  1:194294/1‑149 (MQ=255)
                                                                                                                                          acCGCTCTTGACGGCAACGGATGCAATATCGCACGATCCTTCGCCGAAGATCTCGTGAAATGCCTCCAGCATGGCCTGAATTTTAGGGGGATTGGTGGACGCACAGCCAAGTTGGTGCATAATCTGCATTAGTCAGAAAATGAACGttc       >  1:900471/1‑148 (MQ=255)
                                                                                                                                               tcttGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGAAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagca  <  2:1188733/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACTGAAAGTGCTGTCACCATCGATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCA  >  NC_000913/4633488‑4633779

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: ‑