Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,633,657 C→T 69.6% C30Y (TGC→TAC)  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,6570CT69.6% 420.6 / 264.3 294C30Y (TGC→TAC) yjjXnon‑canonical purine NTP phosphatase, ITPase/XTPase
Reads supporting (aligned to +/- strand):  ref base C (46/43);  new base T (91/113);  total (138/156)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

Only 100 of 294 total aligned reads displayed.

GATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACA  >  NC_000913/4633509‑4633804
                                                                                                                                                    |                                                                                                                                                   
gATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCa                                                                                                                                                    <  2:379068/150‑1 (MQ=255)
           cAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCACCGGATGCAATATGGCAGGATCCTTCGc                                                                                                                                         >  1:2674639/1‑150 (MQ=255)
           cAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGc                                                                                                                                         >  1:548354/1‑150 (MQ=255)
            aaTCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGcc                                                                                                                                        >  2:1134181/1‑150 (MQ=255)
               cGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGaa                                                                                                                                     <  2:1573339/150‑1 (MQ=255)
               cGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGaa                                                                                                                                     <  2:2447091/150‑1 (MQ=255)
               cGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGaa                                                                                                                                     <  2:3934514/150‑1 (MQ=255)
               cGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGaa                                                                                                                                     <  2:4262287/150‑1 (MQ=255)
               cGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGaa                                                                                                                                     <  2:2444462/150‑1 (MQ=255)
               cGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGaa                                                                                                                                     <  2:2361198/150‑1 (MQ=255)
                 ccACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGa                                                                                                                                   >  1:3449005/1‑150 (MQ=255)
                  cACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAt                                                                                                                                  <  2:2519703/150‑1 (MQ=255)
                   aCCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGAtc                                                                                                                                 >  1:3399854/1‑150 (MQ=255)
                      cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCg                                                                                                                              >  1:1710545/1‑150 (MQ=255)
                      cAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCg                                                                                                                              >  1:3924027/1‑150 (MQ=255)
                       aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGt                                                                                                                             >  2:59111/1‑150 (MQ=255)
                       aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGt                                                                                                                             <  1:2778264/150‑1 (MQ=255)
                       aaaaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGt                                                                                                                             <  1:2264717/150‑1 (MQ=255)
                         aaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGa                                                                                                                           <  2:4136563/150‑1 (MQ=255)
                         aaaTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGa                                                                                                                           <  2:2784242/150‑1 (MQ=255)
                           aTCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGaaa                                                                                                                         <  1:4207176/150‑1 (MQ=255)
                            tCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAAt                                                                                                                        <  1:1266646/150‑1 (MQ=255)
                              aGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGc                                                                                                                      >  1:2107811/1‑150 (MQ=255)
                                 ctctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTg                                                                                                                   <  1:3722011/150‑1 (MQ=255)
                                  tctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGc                                                                                                                  >  2:4059536/1‑150 (MQ=255)
                                  tctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGc                                                                                                                  >  1:1214311/1‑150 (MQ=255)
                                   ctGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCa                                                                                                                 >  2:2636182/1‑150 (MQ=255)
                                    tGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAg                                                                                                                >  1:1028477/1‑150 (MQ=255)
                                     ggAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGa                                                                                                               <  1:1925854/150‑1 (MQ=255)
                                      gAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGaa                                                                                                              >  2:683492/1‑150 (MQ=255)
                                      gAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGaa                                                                                                              <  2:1990909/150‑1 (MQ=255)
                                      gAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACAGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGaa                                                                                                              <  2:329575/150‑1 (MQ=255)
                                          tAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGc                                                                                                          >  1:1484389/1‑150 (MQ=255)
                                            aaCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCt                                                                                                        >  2:1293098/1‑150 (MQ=255)
                                            aaCGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCt                                                                                                        >  1:2520260/1‑150 (MQ=255)
                                            aaCGGCGGGCAGTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCCAAGGGCTGTTCCGGTACACCGCTCGCGACGGCAACGGATGCAATCTGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCt                                                                                                        >  2:290228/1‑150 (MQ=255)
                                              cggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGa                                                                                                      >  1:3891441/1‑150 (MQ=255)
                                               ggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGaa                                                                                                     <  2:4282296/150‑1 (MQ=255)
                                               ggcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGaa                                                                                                     >  2:3762729/1‑150 (MQ=255)
                                                gcggGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAAt                                                                                                    <  1:2942389/150‑1 (MQ=255)
                                                    gCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTa                                                                                                >  2:779667/1‑150 (MQ=255)
                                                     cATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAg                                                                                               >  1:1589023/1‑150 (MQ=255)
                                                      aTTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGc                                                                                              <  1:611922/150‑1 (MQ=255)
                                                           cTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGAt                                                                                         >  2:2933998/1‑150 (MQ=255)
                                                           cTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGAt                                                                                         >  1:1459871/1‑150 (MQ=255)
                                                           cTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGAt                                                                                         >  2:3512739/1‑150 (MQ=255)
                                                             aCCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtg                                                                                       >  1:2368720/1‑150 (MQ=255)
                                                              cccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtgg                                                                                      >  2:221101/1‑150 (MQ=255)
                                                              cccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtgg                                                                                      >  2:757655/1‑150 (MQ=255)
                                                               ccGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATtggt                                                                                     <  1:772431/150‑1 (MQ=255)
                                                                   tttCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTc                                                                                 <  1:2594475/150‑1 (MQ=255)
                                                                       gTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcac                                                                             <  1:1667852/150‑1 (MQ=255)
                                                                       gTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcac                                                                             >  2:2057204/1‑150 (MQ=255)
                                                                       gTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcac                                                                             >  2:988137/1‑150 (MQ=255)
                                                                        tGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGcaca                                                                            <  1:206377/150‑1 (MQ=255)
                                                                          cgcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                                                          >  1:416837/1‑150 (MQ=255)
                                                                          cgcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                                                          <  2:504548/150‑1 (MQ=255)
                                                                          cgcCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCCCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGa                                                                          <  1:859787/150‑1 (MQ=255)
                                                                                cgcgcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACtt                                                                    <  2:15811/150‑1 (MQ=255)
                                                                                   gcgTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGt                                                                 <  1:3152291/150‑1 (MQ=255)
                                                                                     gTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGc                                                               >  2:3615872/1‑150 (MQ=255)
                                                                                       ttCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCAt                                                             >  1:4964/1‑150 (MQ=255)
                                                                                         ccTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATaa                                                           <  2:623885/150‑1 (MQ=255)
                                                                                          cTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAAt                                                          >  1:1540502/1‑150 (MQ=255)
                                                                                          cTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAAt                                                          <  1:87207/150‑1 (MQ=255)
                                                                                           tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                                         >  2:2667540/1‑150 (MQ=255)
                                                                                           tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCCGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                                         <  1:1022436/150‑1 (MQ=255)
                                                                                           tCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATc                                                         >  2:3567967/1‑150 (MQ=255)
                                                                                             aCTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAg                                                       >  2:388405/1‑150 (MQ=255)
                                                                                             aCTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAg                                                       >  2:1333876/1‑150 (MQ=255)
                                                                                              cTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGc                                                      >  2:1885665/1‑150 (MQ=255)
                                                                                               tGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCa                                                     >  1:1496680/1‑150 (MQ=255)
                                                                                                 ccAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTATAGCGGGAATGTTGGTCGCACAGACAACATGGAGAATAATCAGCAtt                                                   >  2:1024617/1‑150 (MQ=255)
                                                                                                 ccAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCAtt                                                   <  2:3224267/150‑1 (MQ=255)
                                                                                                      gggCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCa                                              <  2:1487778/150‑1 (MQ=255)
                                                                                                      gggCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCa                                              <  1:2031070/150‑1 (MQ=255)
                                                                                                        gCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGa                                            <  1:681825/150‑1 (MQ=255)
                                                                                                           gTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGaaaa                                         <  2:3674414/150‑1 (MQ=255)
                                                                                                           gTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGaaaa                                         >  1:2334619/1‑150 (MQ=255)
                                                                                                               cGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTaa                                     <  1:3602096/150‑1 (MQ=255)
                                                                                                               cGGTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTaa                                     <  1:334475/150‑1 (MQ=255)
                                                                                                                gttACACCGCTCTCGAAGTCACCGGATGCAATATTTTAGGATCATTAGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAAc                                    <  2:2313964/148‑1 (MQ=255)
                                                                                                                ggTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAAc                                    >  1:3789001/1‑150 (MQ=255)
                                                                                                                ggTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAAc                                    <  1:3023906/150‑1 (MQ=255)
                                                                                                                ggTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAAc                                    <  1:2818702/150‑1 (MQ=255)
                                                                                                                ggTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAAc                                    <  1:1134357/150‑1 (MQ=255)
                                                                                                                 gTACACCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACg                                   >  2:1264537/1‑150 (MQ=255)
                                                                                                                     acCGCTCTCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAc                               >  2:3002935/1‑150 (MQ=255)
                                                                                                                         ctctCGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagca                           >  2:53024/1‑150 (MQ=255)
                                                                                                                          tctcGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagcag                          >  1:3694904/1‑150 (MQ=255)
                                                                                                                          tctcGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTAcagcag                          >  2:214574/1‑150 (MQ=255)
                                                                                                                           ctcGACGGCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGt                         <  2:3388007/150‑1 (MQ=255)
                                                                                                                               aCGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGtataa                     <  2:1206982/150‑2 (MQ=255)
                                                                                                                                  gCAACGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGa                  <  1:3986131/150‑1 (MQ=255)
                                                                                                                                    aaCGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaa                >  1:4098238/1‑150 (MQ=255)
                                                                                                                                     aCGGATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaa               >  2:1781238/1‑150 (MQ=255)
                                                                                                                                       ggATGCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaa             >  1:1765983/1‑150 (MQ=255)
                                                                                                                                        gATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGaaaaaaa            <  1:1145012/150‑1 (MQ=255)
                                                                                                                                           gCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCa         <  1:2217412/150‑1 (MQ=255)
                                                                                                                                           gCAATATGGTAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCa         >  1:4269300/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GATGCCAGCTTCAATCGCCACCCAAAAATCAGCCTCTGGAAGTAAACGGCGGGCATTGGCTACCCGATTTCGTGCGCCAGCGCGCGTTTCCTCACTGCCAAAGGGCTGTTCCGGTACACCGCTCTCGACGGCAACGGATGCAATATGGCAGGATCCTTCGCCGAAGATCTCGTGAAATGCCTGCAGAATGGCCTGAATTTTAGCGGGATTGGTGGTCGCACAGACAACTTGGTGCATAATCAGCATTACTCAGAAAATTAACGTTACAGCAGTATACGGAAAAAAAGCATGTTACA  >  NC_000913/4633509‑4633804

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