Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,599,236 C→A 9.2% E175D (GAG→GAT dapA ← 4‑hydroxy‑tetrahydrodipicolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,599,2360CA 9.2% 282.2 / 13.5 87E175D (GAG→GATdapA4‑hydroxy‑tetrahydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base C (33/46);  new base A (2/6);  total (35/52)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

TCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAG  >  NC_000913/2599093‑2599384
                                                                                                                                               |                                                                                                                                                    
tCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGc                                                                                                                                                      <  2:192524/144‑1 (MQ=255)
tCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGc                                                                                                                                                      >  1:192524/1‑144 (MQ=255)
tCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTg                                                                                                                                                >  1:149755/1‑150 (MQ=255)
tCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTg                                                                                                                                                <  1:424288/150‑1 (MQ=255)
tCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTg                                                                                                                                                <  2:1043709/150‑1 (MQ=255)
 ctgctgCCAGTGTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGa                                                                                                                                               <  2:123711/150‑1 (MQ=255)
  tgctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGAt                                                                                                                                              >  2:762343/1‑150 (MQ=255)
   gctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATc                                                                                                                                             <  2:612264/150‑1 (MQ=255)
   gctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATc                                                                                                                                             <  1:668748/150‑1 (MQ=255)
     tgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTg                                                                                                                                           <  2:943854/150‑1 (MQ=255)
      gcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTgg                                                                                                                                          >  1:1133019/1‑150 (MQ=255)
            ttGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCAAAGAGCTGGtttcag                                                                                                                                    <  2:1228159/150‑4 (MQ=255)
            ttGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCAAAGAGCTGGtttcag                                                                                                                                    <  1:625055/150‑4 (MQ=255)
             tGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCAAAGAGCTGGtttcaga                                                                                                                                   <  1:900655/150‑5 (MQ=255)
             tGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCAAAGAGCTGGtttcaga                                                                                                                                   <  1:651763/150‑5 (MQ=255)
              gCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcg                                                                                                                                  >  2:106905/1‑150 (MQ=255)
              gCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCAAAGAGCTGGtttcagat                                                                                                                                  <  2:1109727/150‑6 (MQ=255)
               cacaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgc                                                                                                                                 >  1:157004/1‑150 (MQ=255)
                acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                <  2:107959/150‑1 (MQ=255)
                 caTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGt                                                                                                                               <  1:1295899/150‑1 (MQ=255)
                  ttCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACTGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGtt                                                                                                                              <  2:421651/149‑1 (MQ=255)
                     tGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAg                                                                                                                           <  2:1145579/150‑1 (MQ=255)
                          cATATCACCCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTTCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTccc                                                                                                                      <  2:751479/150‑1 (MQ=255)
                           atatCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCt                                                                                                                     <  1:131797/150‑1 (MQ=255)
                               cACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTg                                                                                                                 >  1:1535/1‑150 (MQ=255)
                               cACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTg                                                                                                                 <  2:50692/150‑1 (MQ=255)
                                aCGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGc                                                                                                                >  1:836847/1‑150 (MQ=255)
                                 cgcgCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGcc                                                                                                               <  1:370799/150‑1 (MQ=255)
                                       gCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTGCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTg                                                                                                         <  2:148626/150‑1 (MQ=255)
                                       gCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTg                                                                                                         <  2:224538/150‑1 (MQ=255)
                                               aGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGat                                                                                                 <  1:1372789/150‑1 (MQ=255)
                                                gCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGata                                                                                                <  1:762343/150‑1 (MQ=255)
                                                 ccGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGataa                                                                                               <  1:178439/150‑1 (MQ=255)
                                                  cGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGataat                                                                                              >  2:1350487/1‑150 (MQ=255)
                                                  cGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGataat                                                                                              <  1:1410531/150‑1 (MQ=255)
                                                      aCGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAttt                                                                                          >  2:895295/1‑150 (MQ=255)
                                                       cGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttt                                                                                         >  1:1225400/1‑150 (MQ=255)
                                                         gAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttttt                                                                                       <  2:994814/150‑1 (MQ=255)
                                                           aaTAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTAc                                                                                     <  1:1360642/150‑1 (MQ=255)
                                                           aaTAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTAc                                                                                     <  2:807788/150‑1 (MQ=255)
                                                            aTAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACt                                                                                    >  1:353045/1‑150 (MQ=255)
                                                                        ccGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCCGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGCCgg                                                                        >  1:614790/1‑150 (MQ=255)
                                                                        ccGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACgg                                                                        >  2:990517/1‑150 (MQ=255)
                                                                                  tcATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                                                              >  2:1232829/2‑150 (MQ=255)
                                                                                  tcATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                                                              >  2:1226974/2‑150 (MQ=255)
                                                                                  gCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACCCGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                                                              >  1:1307240/1‑150 (MQ=255)
                                                                                  gCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                                                              >  2:1544493/1‑150 (MQ=255)
                                                                                   cATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTc                                                             >  2:1003957/1‑150 (MQ=255)
                                                                                    aTGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTcc                                                            >  1:1512642/1‑150 (MQ=255)
                                                                                       aaGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                                         <  2:69400/150‑1 (MQ=255)
                                                                                       aaGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGGTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                                         <  1:1076451/150‑1 (MQ=255)
                                                                                         gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAg                                                       <  2:353045/150‑1 (MQ=255)
                                                                                          tCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGc                                                      <  2:1225400/150‑1 (MQ=255)
                                                                                           ccAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCa                                                     <  1:101278/150‑1 (MQ=255)
                                                                                           ccAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCa                                                     >  1:823954/1‑150 (MQ=255)
                                                                                           ccAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCa                                                     <  1:329379/150‑1 (MQ=255)
                                                                                             aGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGa                                                   >  1:1100829/1‑150 (MQ=255)
                                                                                              gctcgcTCGCATCAGCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGGTTACACGCGTTAAGTTCCCTGTTGCGGCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGAt                                                  <  2:1537268/147‑1 (MQ=255)
                                                                                              gcgcgcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGAt                                                  <  2:1195222/150‑1 (MQ=255)
                                                                                              gcgcgcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGAt                                                  <  1:367693/150‑1 (MQ=255)
                                                                                               cgcgcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATc                                                 <  1:828290/150‑1 (MQ=255)
                                                                                                gcgcTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCg                                                <  2:1512642/150‑1 (MQ=255)
                                                                                                   cTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAg                                             >  2:1350846/1‑150 (MQ=255)
                                                                                                    tCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGc                                            >  1:800156/1‑150 (MQ=255)
                                                                                                     cGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGcc                                           <  1:1167631/150‑1 (MQ=255)
                                                                                                     cGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGcc                                           <  2:149755/150‑1 (MQ=255)
                                                                                                      gCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCa                                          <  2:1133019/150‑1 (MQ=255)
                                                                                                          catcGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTAc                                      <  1:1003957/150‑1 (MQ=255)
                                                                                                          catcGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTAc                                      >  1:459909/1‑150 (MQ=255)
                                                                                                          catcGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTAc                                      >  1:465390/1‑150 (MQ=255)
                                                                                                                   tCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGca                             >  1:930065/1‑150 (MQ=255)
                                                                                                                   tCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGca                             <  1:1350846/150‑1 (MQ=255)
                                                                                                                      gcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGt                                        >  2:1290078/1‑136 (MQ=255)
                                                                                                                      gcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGt                                        <  1:1290078/136‑1 (MQ=255)
                                                                                                                      gcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACAt                          <  2:1535/150‑1 (MQ=255)
                                                                                                                       cagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACAtt                         >  1:445155/1‑150 (MQ=255)
                                                                                                                          aaCAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTCTTGCCTCTTTGATTCCCATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATtata                      >  1:1078233/1‑150 (MQ=255)
                                                                                                                             aaaaTCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAg                   <  2:157004/150‑1 (MQ=255)
                                                                                                                              aaaTCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGa                  >  1:71957/1‑150 (MQ=255)
                                                                                                                               aaTCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGaa                 <  1:1130995/150‑1 (MQ=255)
                                                                                                                                atcatcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAAt                <  1:895295/150‑1 (MQ=255)
                                                                                                                                    gtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGc            >  2:1471723/5‑150 (MQ=255)
                                                                                                                                       taaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGc         <  1:97819/149‑1 (MQ=255)
                                                                                                                                        aaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCa        >  2:899197/1‑150 (MQ=255)
                                                                                                                                          aCCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAgg      >  2:1282974/1‑150 (MQ=255)
                                                                                                                                          aCCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAgg      >  2:540538/1‑150 (MQ=255)
                                                                                                                                            cAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTc    >  1:407224/1‑150 (MQ=255)
                                                                                                                                            cAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTc    <  1:540538/150‑1 (MQ=255)
                                                                                                                                              gCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAg  <  1:1544493/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAG  >  NC_000913/2599093‑2599384

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