Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,599,236 C→A 8.8% E175D (GAG→GAT dapA ← dihydrodipicolinate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,599,2360CA 8.8% 269.1 / 13.5 98E175D (GAG→GATdapAdihydrodipicolinate synthase
Reads supporting (aligned to +/- strand):  ref base C (43/46);  new base A (6/2);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

Only 100 of 103 total aligned reads displayed.

CCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAG  >  NC_000913/2599089‑2599384
                                                                                                                                                   |                                                                                                                                                    
cccTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGCCCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAATATCATCTGAAACCAGctg                                                                                                                                                    >  2:1266624/1‑149 (MQ=255)
cccTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGCACAAAATCATCTGAAACCAGctc                                                                                                                                                    >  2:1062418/1‑150 (MQ=255)
cccTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCAGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATAGCATGAAGTCCAGCCCGCTCGCATCATCGCCGCTCACCAGAACAAAATCATCTGAAACCAGatc                                                                                                                                                    >  2:1698807/1‑150 (MQ=255)
 ccTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGctct                                                                                                                                                   >  1:907286/1‑150 (MQ=255)
  cTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCtt                                                                                                                                                  <  1:816782/150‑1 (MQ=255)
     ctgctgCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGa                                                                                                                                               <  2:352394/150‑1 (MQ=255)
      tgctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCGTTGAt                                                                                                                                              >  2:584232/1‑150 (MQ=255)
       gctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATc                                                                                                                                             >  2:582508/1‑150 (MQ=255)
       gctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATc                                                                                                                                             >  2:610551/1‑150 (MQ=255)
        ctgcCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCt                                                                                                                                            >  2:1023954/1‑150 (MQ=255)
            cAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGtt                                                                                                                                        <  2:1624558/150‑1 (MQ=255)
            cAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGtt                                                                                                                                        <  1:1551508/150‑1 (MQ=255)
             aGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGTAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCATAACAAAATCATCTGAAACCAGCTCTTTGATCTGGttt                                                                                                                                       >  2:852742/1‑150 (MQ=255)
             aGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGttt                                                                                                                                       <  1:610551/150‑1 (MQ=255)
                 tGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTacac                                                                                                                                   >  1:607703/1‑150 (MQ=255)
                 tGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTacac                                                                                                                                   <  1:642339/150‑1 (MQ=255)
                  gCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGGAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGATCAAAGAGCTGGtttcagat                                                                                                                                  <  1:702411/150‑6 (MQ=255)
                   cacaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgc                                                                                                                                 <  1:995153/150‑1 (MQ=255)
                   cacaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACAAGCTCTTTGATCTGGTTTacactc                                                                                                                                 >  2:1216209/1‑148 (MQ=255)
                    acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATCACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                <  2:816647/150‑1 (MQ=255)
                    acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                <  2:607703/150‑1 (MQ=255)
                    acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                >  1:1180375/1‑150 (MQ=255)
                    acaTCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACAcgcg                                                                                                                                <  2:371152/150‑1 (MQ=255)
                          gggCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCCGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCCGCTCTGTGATCTGGTTATCACGCGTTAAGt                                                                                                                          >  2:1558171/1‑150 (MQ=255)
                           ggCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGtt                                                                                                                         <  1:657854/150‑1 (MQ=255)
                            gCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGAATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGACCAAAATCATCTGAAACCAGCTCTTTGATCAGGTGTTCACGCGTTAAGTTc                                                                                                                        >  1:260830/1‑150 (MQ=255)
                             ccATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCCTGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGATcc                                                                                                                       <  2:1040254/150‑1 (MQ=255)
                                tatCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCACTg                                                                                                                    >  2:596676/1‑150 (MQ=255)
                                    aCGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGc                                                                                                                >  2:634588/1‑150 (MQ=255)
                                     cgcgCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGcc                                                                                                               >  1:464081/1‑150 (MQ=255)
                                       cgcTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCctc                                                                                                             <  1:642022/150‑1 (MQ=255)
                                         cTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCtt                                                                                                           <  1:511454/150‑1 (MQ=255)
                                          tGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCttt                                                                                                          <  2:470564/150‑1 (MQ=255)
                                            cGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGCACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGa                                                                                                        >  1:710846/1‑150 (MQ=255)
                                            cGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGa                                                                                                        <  1:504204/150‑1 (MQ=255)
                                              aCGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGAtt                                                                                                      <  2:1291987/150‑1 (MQ=255)
                                                gTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTcc                                                                                                    <  2:2776/150‑1 (MQ=255)
                                                gTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTcc                                                                                                    >  2:1346327/1‑150 (MQ=255)
                                                gTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGGTAAGTTCCCTGTTGCCTCTTTGATTcc                                                                                                    >  2:66461/1‑150 (MQ=255)
                                                  tAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGa                                                                                                  <  1:596676/150‑1 (MQ=255)
                                                   aGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCATCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGat                                                                                                 <  2:302335/150‑1 (MQ=255)
                                                       gTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGataata                                                                                             <  1:966254/150‑1 (MQ=255)
                                                           cGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttt                                                                                         >  1:1324318/1‑150 (MQ=255)
                                                            ggAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAttttt                                                                                        <  2:1296255/150‑1 (MQ=255)
                                                            ggAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAttttt                                                                                        <  2:464081/150‑1 (MQ=255)
                                                             taaaTAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCATCTCTTTGATCTGGTTTACACGCGTTTAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttttt                                                                                       <  1:695908/149‑1 (MQ=255)
                                                             gAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCTCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttttt                                                                                       <  1:94624/150‑1 (MQ=255)
                                                             gAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATAtttttt                                                                                       <  2:150331/150‑1 (MQ=255)
                                                                 tAACCCCATGACCGCCCAATTGCATGAAGTCAAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCACTGTTGCCTCTTTGATTCCGATAATATTTTTTTCtt                                                                                   >  1:134622/1‑150 (MQ=255)
                                                                    ccccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCg                                                                                <  1:192138/150‑1 (MQ=255)
                                                                    ccccATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCg                                                                                <  1:625036/150‑1 (MQ=255)
                                                                     cccATGACCGCCCAATTGCATGAAGTGCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGc                                                                               <  1:1097958/150‑1 (MQ=255)
                                                                        aTGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAg                                                                            <  1:1031215/150‑1 (MQ=255)
                                                                               cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGccc                                                                     >  2:659329/1‑150 (MQ=255)
                                                                               cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGccc                                                                     <  2:452225/150‑1 (MQ=255)
                                                                               cccAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGccc                                                                     >  2:426775/1‑150 (MQ=255)
                                                                                 cAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCAc                                                                   <  1:1023954/150‑1 (MQ=255)
                                                                                      gCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                                                              >  1:1232802/1‑150 (MQ=255)
                                                                                      gCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAAACATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAACTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGttt                                                              >  2:462226/1‑150 (MQ=255)
                                                                                        aTGAATTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAATATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTcc                                                            <  2:1505457/150‑1 (MQ=255)
                                                                                           aaGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                                         <  2:893006/150‑1 (MQ=255)
                                                                                            aGTCCAGCGCCCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTAGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTATCGCCAGACGGCCCACCGTTTCCGGGa                                                        >  2:891030/1‑150 (MQ=255)
                                                                                             gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                                         <  2:1290517/148‑1 (MQ=255)
                                                                                             gTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCggg                                                         >  1:1290517/1‑148 (MQ=255)
                                                                                               ccAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCa                                                     >  1:1191380/1‑150 (MQ=255)
                                                                                                       cTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGAGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAg                                             >  2:733545/1‑150 (MQ=255)
                                                                                                         cGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCGTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCCCCGTTTCCGGGCGCAGATCGCAGcc                                           >  2:1254680/1‑150 (MQ=255)
                                                                                                           catcatCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAg                                         <  1:582508/150‑1 (MQ=255)
                                                                                                               atcGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTAGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCAAGACGGCACACCGTTTCAGGGAGCAGATAGCAGCCAGTACg                                     >  1:1233386/1‑150 (MQ=255)
                                                                                                                tcGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACgg                                    <  1:66461/150‑1 (MQ=255)
                                                                                                                 cgccgcTCAGCAGCACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACggg                                   >  1:1339815/1‑150 (MQ=255)
                                                                                                                     gcTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGATAAGTTACCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCTCCAGACGGCCCACAGTTTCAGGGAGCAGCTCCAACCAAGTACGGGACgg                               >  2:336819/1‑150 (MQ=255)
                                                                                                                     gcTCACCAGAACAAAATCATCTGAAACCAGCCCTTTGATCTGGTGTACACGCGTGAAGTACCCTGGTGCCTCTTTGATGCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTGCCGGGAGCAGATCGCAGCCAGTACGGGACgg                               <  1:141110/150‑1 (MQ=255)
                                                                                                                        cagcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGcac                            >  2:437773/1‑150 (MQ=255)
                                                                                                                        cagcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGAGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGcac                            >  2:473182/1‑150 (MQ=255)
                                                                                                                          gcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACAt                          <  2:727335/150‑1 (MQ=255)
                                                                                                                          gcagAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACAt                          >  1:34709/1‑150 (MQ=255)
                                                                                                                             gAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCCGCCAGTACGGGACGGCACATtat                       >  1:1552642/1‑150 (MQ=255)
                                                                                                                                    atcatcTGAAACCAGCTCTTTCATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAAt                <  1:1120623/150‑1 (MQ=255)
                                                                                                                                     tgatcTGAATCCAGTTCTTTGATGTGGTTTACTGGGGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCTGCCAGTACGGGACGGCACATTATACAGAttt               <  1:1698807/148‑1 (MQ=255)
                                                                                                                                     tcatcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAAtt               <  1:426775/150‑1 (MQ=255)
                                                                                                                                     cgcgtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTTCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGCACGGCACATTATACAGAAtt               >  2:50906/8‑150 (MQ=255)
                                                                                                                                      gcgtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGCCGGCACATTATACAGAAttt              >  1:1631496/7‑150 (MQ=255)
                                                                                                                                       cgtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTg             >  2:1497525/6‑150 (MQ=255)
                                                                                                                                       cgtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTg             >  2:457540/6‑150 (MQ=255)
                                                                                                                                       atcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTg             >  2:216246/1‑150 (MQ=255)
                                                                                                                                        tcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGc            >  1:1415632/1‑150 (MQ=255)
                                                                                                                                        tcTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGAGAATATTTTTTACTTTCGCCCGACGGCCCCCCGTGTCCGGGAGCAGATCGCGGCCAGAACGGGACGGCACATTAGACAGACGTTGc            >  2:722679/1‑150 (MQ=255)
                                                                                                                                        gtgtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGc            >  2:238713/5‑150 (MQ=255)
                                                                                                                                         cTGAAACCCGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCGCACCGTTTCCGGGAGCAGATCGCCGCCAGTACGCGACGGCCCATTATAGAGAAGTTGCg           >  2:861933/1‑150 (MQ=255)
                                                                                                                                         cTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCg           <  1:854588/150‑1 (MQ=255)
                                                                                                                                         cTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCAAGTACGGGACGGCACATTATACAGAATTTGCg           >  1:1179605/1‑150 (MQ=255)
                                                                                                                                         cTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGGTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAAGTTGCg           >  2:1021482/1‑150 (MQ=255)
                                                                                                                                          gtaaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGTCACATTATACAGAATTTGCgg          <  1:1497525/148‑1 (MQ=255)
                                                                                                                                            ataCCATCTCTTTGATCTGTTTTACACGCGTTAAGTTCACTGTTGCCTCTTTTCTTCAGAAAATATTTTTAACTTTAGCCAGACTTCCCACAGTTTCCGGTAGCAGATCGCAGACAGTACGGGACGGCACATTATACAGAATTTGCGGCa        <  2:1018318/148‑1 (MQ=255)
                                                                                                                                            aaaCCAGATCAAAGAGCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGAAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCa        >  1:1700333/1‑150 (MQ=255)
                                                                                                                                             aaCCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAg       <  1:854338/150‑1 (MQ=255)
                                                                                                                                              accttctctTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTTCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCTCCGTTTCCGCGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAgg      <  1:447152/145‑1 (MQ=255)
                                                                                                                                                cAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTc    <  1:1062418/150‑1 (MQ=255)
                                                                                                                                                  gCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAg  <  2:624441/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCCTTCTGCTGCCAGTTTGCACATCTGGGCCATATCACGCGCTGCGACGTTAGCCGTAACGGAAATAACCCCATGACCGCCCAATTGCATGAAGTCCAGCGCGCTCGCATCATCGCCGCTCAGCAGAACAAAATCATCTGAAACCAGCTCTTTGATCTGGTTTACACGCGTTAAGTTCCCTGTTGCCTCTTTGATTCCGATAATATTTTTTACTTTCGCCAGACGGCCCACCGTTTCCGGGAGCAGATCGCAGCCAGTACGGGACGGCACATTATACAGAATTTGCGGCAGGTCAG  >  NC_000913/2599089‑2599384

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: